ue  FILE  COPN  AD-A177  649 


APOSR.TK*  87-0  190 


Approved  for  public  release* 

ion  unlimited. 


4^ 


w  AT,  er.TTM-TyTC  RESEARCH  (AfSC) 

m  ^  ‘T.  - :  iosuo  l3 

•'.f  ‘  .  v,  - -i  r^.fi.rwed  •'no  13 

This  i  •  ■  ,lr  1  ! n  IA.V  AF/l  19>l2‘ 

t?„ . . ilimitcd. 

Bist”  t  • _ ; 

TT  T  7  ^ 


IKATTin 


«  :--o 


FCElef .  Iconic- 

■  |  m  |  September  1986 

*T 


rmation  Division 


LIDS-TH- 1599 


ON  ROBUST  CONTROL  DESIGNS  FOR 
INFINITE  DIMENSIONAL  SYSTEMS 
by 

Wing  Hong  Ricky  Lee 


DTIC 


i 

1 

» 

Laboratory  for  Information  and  Decision  Systems  J 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY,  CAMBRIDGE,  MASSACHUSEnS  02139  | 


f 

! 


I  DISTRIBUTION  statement  a~1 

I  Approved  fox  public  rrloam  S7  2  27  0  0  1 

'  DUtributk.n  TJnl>>  ' 

’ 

wicwtji  ininTmMMMilMiamBsai  T.f.ivr.ixr-rjr.  attlifcii  *fJ  Mi  f  g  |'|  f|  f iJiXaXlXlXiJlljXLCLtlL^tLCL^Ai 


Unclassified 

SECURITY  CLASSIFICATION  Of  THIS  PAPE 


’|)^^sWv  CLA”,F,CAT,ON 


m  6</ i 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


J*.  SECURITY  CLASSIFICATION  AUTHORITY 


3b.  DECLASSIFICATION/DOWNG  RACING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 


3.  DISTRIBUTIONS VAI LABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AF08R.TR.  8  7  -  0  i  90 


Sa.  NAME  OF  PERFORMING  ORGANIZATION  Bb.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

(It  nppllcabla) 

Massachusetts  Inst,  of  Tech.  AFOSR/NM 


6c.  ADDRESS  (City.  State  and  ZIP  Coda) 

Room  35-311 
Cambridge,  MA  02139 


Sa.  NAME  OF  FUNOINQ/SPONSORING 
ORGANIZATION 

AFOSR 


*fit  Slatt  ond  ZIP  Co<u> 

Bolling  AFB  DC  20332-6448 


7b.  ADDRESS  (City.  Slata  and  ZIP  Coda) 

Bldg  410 

Bolling  AFB  DC  20332-6448 


8b.  OFFICE  SYMBOL  S.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

NM  AF0SR-82-0258 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 

61102F 


M\ 


4im 


’MWBttWtfcnj  Ricky  Lee 


13a.  TYPE  OF  REPORT 

Reprint 


13b.  TIME  COVERED 


BlUBgaiiil 


COSAT I  COOES 


18.  SUBJECT  TERMS  (Continue  on  reuerae  if  neceaaory  ond  identify  by  block  number ) 


\v 


SUB-  GR. 


|  This  thesis  deals  with  the  robustness  of  stability  of 

'I'e’eMTOirT  distributed.  linear- time- Invariant  (XTI)  feedback  control  aysteas. 

IS.  ABSTRACT  iConimua  on  ret  ^  mtn  gQal  „  t#  forM|,BM  a  praitlc.{  „thod  for  .valuatll* 

feedback  designs  based  on  the  actual  XTI  ay, tea  clmracterlstics. 

As  a  result,  a  design  procedure  can  be  developed  for  XTI  aysteas  to 
synthesize  feedback  controllers  that  are  guaranteed  to  be 
closed-loop  stable.  We  have  developed  a  robustness  characterization 
for  DLTI  systems,  and  have  shown  that  linear  quadratic  (LQ)  optimal 
control  systems  have  nice  robustness  properties  and  can  serve  as 
good  reference  designs  for  the  actual  Implementation  of  the  feedback 
controller.  We  have  studied  In  detail  linear  hereditary 
differential  systems  and  a  vibration  suppression  problem  for  a 
flexible  beam.  We  stress  the  study  of  lmplemen table  controllers, 
which  are  finite  dimensional.  In  contrast  to  optimal  controllers 
that  are  typically  lnflnlte-dlmenslonal.  However,  one  can  integrate 
uur  multivariable  robustness  results  with  the  LQ  optimal  control  to 
derive  a  finite-dimensional  suboptlmal  control  law  which  Is 
closed  loop  stable.  We  shoe  hew  this  can  be  dene  by  ualr« 
spatially-sampled  measurements  along  the  flexible  beam.  Also  we 
have  demonstrated  by  using  our  robustness  results  that  the  inherent 
damping  In  the  flexible  structure  plays  a  vital  role  In  determining 
whether  a  physically  realizable  closed-loop  stabilizing  controller, 
based  on  spatially-sampled  measurements,  can  be  achieved. 


30.  OISTRI BUTION/A  VAI  LABI  LIT  Y  OF  ABSTRACT 

UNCLASSIF I EO/UN LIMITED  □  SAME  AS  RPT.  □  OTIC  USERS  □ 


33s.  NAME  OF  RESPONSIBLE  INOIVIOUAL 

Major  Crowley 


DD  FORM  1473,  83  APR  EDITION  OF  1  JAN  73  IS  OBSOLETE. 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b.  TELEPHONE  NUMBER 
(Includa  Area  Code) 

(202)  767-5026 


22c.  OFFICE  SYMBOL 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


September  1986 


LIDS-TH-1599 


ON  ROBUST  CONTROL  DESIGNS  FOR 
INFINITE  DIMENSIONAL  SYSTEMS 

by 

Wing  Hong  Ricky  Lee 


This  report  is  based  on  the  unaltered  thesis  of  Wing  Hong  Ricky  Lee. 
submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy  at  the  Massachusetts  Institute  of  Technology, 
Laboratory  for  Information  and  Decision  Systems  with  partial  support 
provided  by  the  Air  Force  Office  of  Scientific  Research  Grant 
no.  AF0SR-82-0258.  and  NASA  Grant  no.  NAG  2-297. 


Laboratory  for  Information  and  Decision  Systems 
Massachusetts  Inrtitute  of  Technology 
Cambridge,  Massachusetts  02139 


Accesion  For 


NTIS  CRA&I 
DTIC  TAB 

Unanno.i  >ccd 
Justiticatio.i 


£ 


□ 

□ 


By . 

Dist.  ibutio../ 


Availability  Codes 


Dist 

A  / 


Avail  a-d/or 
Special 


ON  ROBUST  CONTROL  DESIGNS  FOR 
INFINITE  DIMENSIONAL  SYSTEMS 
by 

Wing  Hong  Ricky  Lee 


Submitted  to  the  Department  of  Electrical  Engineering 
and  Computer  Science  on  May  23,  1986  in  par ti cal 
fulfillment  of  the  requirements  for  the  Degree  of 
Doctor  of  Philosophy 

Abstract 

\  ' 

— This  thesis  deals  with  the  robustness  of  stability  of 
distributed,  linear-time-lnvarlant  (DLTI)  feedback  control  systems. 
The  main  goal  is  to  formulate  a  practical  method  for  evaluating 
feedback  designs  based  on  the  actual  DLTI  system  characteristics. 

As  a  result,  a  design  procedure  can  be  developed  for  DLTI  systems  to 
synthesize  feedback  controllers  that  are  guaranteed  to  be 
closed-loop  stable.  We  have  developed  a  robustness  characterization 
for  DLTI  systems,  and  have  shown  that  linear  quadratic  (LQ)  optimal 
control  systems  have  nice  robustness  properties  and  can  serve  as 
good  reference  designs  for  the  actual  implementation  of  the  feedback 
controller.  We  have  studied  in  detail  linear  hereditary 
differential  systems  and  a  vibration  suppression  problem  for  a 
flexible  beam.  We  stress  the  study  of  implementable  controllers, 
which  are  finite  dimensional,  in  contrast  to  optimal  controllers 
that  are  typically  infinite-dimensional.  However,  one  can  integrate 
our  multivariable  robustness  results  with  the  LQ  optimal  control  to 
derive  a  finite-dimensional  suboptimal  control  law  which  is 
closed-loop  stable.  We  show  how  this  can  be  done  by  using 
spatially-sampled  measurements  along  the  flexible  beam.  _  Also  we 
have  demonstrated  by  using  our  robustness  results  that1  the  inherent 
damping  in  the  flexible  structure  plays  a  vital  role  in  determining 
whether  a  physically  realizable  closed-loop  stabilizing  controller, 
based  on  spatially-sampled  measurements,  can  be  achieved. 


Thesis  Co-Supervisor: 

Title: 


Dr.  Alan  S.  Willsky 

Professor  of  Electrical  Engineering 


Thesis  Co-Supervisor:  Dr.  Gunter  Stein 

Title:  Adjunct  Professor  of  Electrical  Engineering 


ON  ROBUST  CONTROL  DESIGNS  FOR 
INFINITE  DIMENSIONAL  SYSTEMS 

by 

Wing  Hong  Ricky  Lee 

S.B.E.E.,  Massachusetts  Institute  of  Technology 

(1980) 

S.B.  Math.,  Massachusetts  Institute  of  Technology 

(1980) 

S.M. .  Massachusetts  Institute  of  Technology 

(1980) 

E.E. ,  Massachusetts  Institute  of  Technology 

(1980) 

SUBMITTED  IN  PARTIAL  FULLFILLMENT 
OF  THE  REQUIRMENTS  FOR  THE 
DEGREE  OF 

DOCTOR  OF  PHILOSOPHY 
at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
June.  1986 

c  Wing  Hong  Ricky  Lee,  1986 

The  author  hereby  grants  to  M.I.T.  permission  to  reproduce  and  to 
distribute  copies  of  this  thesis  document  in  whole  or  in  part. 


Signature  of  the  Author. 


Department  of  Electrical  Engineering  and 
Computer  Science,  May  23,  1986 


Accepted  by 


Arthur  C.  Smith 
Chairman.  Departmental  Graduate  Committee 


3 


ACKNOWLEDGEMENTS 

I  would  like  to  express  my  deepest  gratitude  to  Professors 
Gunter  Stein  and  Alan  Willsky.  The  completion  of  this  thesis  is 
made  possible  by  their  care,  Interests,  and  patience.  They  have 
helped  me  to  see  the  real  world  and  understand  the  goals  of  my 
thesis  research. 

I  would  like  to  thank  Prof.  Sanjoy  Mitter  for  serving  as  a 
reader  of  this  thesis.  He  has  made  key  remarks  and  directed  me  to 
the  appropriate  sources  of  information. 

I  am  grateful  to  Dr.  Nils  Sandell,  Jr.,  Prof.  Bernard  Levy. 
Prof.  Michael  Athans,  and  Dr.  Sol  Gully.  Dr.  Sandell  suggested  to 
me  the  flexible  beam  problem  and  the  Lyapunov  control  approach, 
while  I  was  an  employee  at  ALPHATECH,  INC.,  working  under  Dr.  Gully. 
Dr.  Sandell  and  Prof.  Athans  provided  much  of  the  insight  when  this 
thesis  work  was  being  formulated.  The  results  on  LQHD  systems  are 
based  on  earlier  work  I  did  with  Prof.  Levy. 

I  have  benefit  ted  greatly  from  the  many  technical  discussions 
with,  and  the  friendship  of,  a  number  of  LIDS  classmates.  In 
particular,  I  have  learned  a  great  deal  on  robustness  theory  from 
Dr.  Norman  Lehtomaki,  and  on  semigroup  theory  from 
Dr.  Marcel  Coderch.  During  my  last  year  at  MIT,  I  was  fortunate  to 

•  t 

have  Mr.  Cuneyt  Ozveren,  Mr.  Jem  Robin  Rohlicek  and  Mr.  Mike  Chin  as 
my  officemates.  I  could  always  count  on  them  for  technical  help  and 
support.  I  eum  very  grateful  to  all  of  them.  I  am  also  thankful  to 
Mr.  Petros  Kapasouris  for  helping  me  with  MacHoneyX. 


Ms.  Helene  George  and  Mr.  Arthur  Giordan!  deserve  special 
thanks.  Helene  somehow  managed  to  translate  my  totally  unorganized 
drafts  into  this  report,  and  she  did  an  excellent  job  in  typing. 

Mr.  Giordani's  artwork  has  always  been  superb,  even  though  I  rushed 
my  figures  to  him  at  the  last  minute. 

This  research  and  the  preparation  of  this  report  was  partially 
funded  by  the  Air  Force  Office  of  Scientific  Research  Grant 
#AF0SR-82-0258 ,  and  NASA  Grant  #NAG  2-297. 


% 


i 

* 

v 


5 


TABLE  OF  CONTENTS 

Page 


ABSTRACT  2 

ACKNOWLEDGEMENTS  3 

TABLE  OF  CONTENTS  5 

LIST  OF  FIGURES  8 

CHAPTER  1:  INTRODUCTION  10 

1.1  An  Overview  and  Motivation  10 

1.2  A  Summary  of  the  Results  13 

1.3  Organization  of  this  Thesis  16 

CHAPTER  2:  ROBUSTNESS  OF  LINEAR  HEREDITARY  DIFFERENTIAL 

SYSTEMS  19 

2.1  Introduction  19 

2.2  Linear  Quadratic  Hereditary  Differential  Systems 

and  the  Multivariable  Kalman  Frequency  Domain 
Inequality  20 

2.2.1  Linear  Quadratic  Hereditary  Differential 

System  21 

2.2.2  Kalman  Inequality  for  LQHD  Systems  - 

Statement  24 

2.2.3  Kalman  Inequality  for  LQHD  Systems  - 

Derivation  28 

2.3  Robustness  Characterization  of  Linear  Hereditary 

Differential  Feedback  Systems  33 

2.3.1  Nominal  and  Perturbed  LHDF  Systems  34 

2.3.2  Open-Loop  Unstable  Poles  of  LHDF  Systems  38 

2.3.3  Multivariable  Nyquist  Theorem  42 

2.3.4  Robustness  Theorems  for  Hereditary 

Differential  Systems  42 

2.4  Robustness  Properties  of  Linear  Quadratic 

Hereditary  Differential  Systems  49 


2.5  Discussion  on  LQHD  Control  Designs  56 

2.5.1  Examples  of  LQHD  Systems  56 

2.5.2  Extensions  of  LQHD  Control  Designs  60 

CHAPTER  3:  ROBUSTNESS  PROPERTIES  OF  DISTRIBUTED, 

LINEAR-TIME-INVARIANT  SYSTEMS  63 

3.1  Introduction  63 

3.2  Yakubovich  Frequency  Domain  Theorem  65 

3.2.1  Evolution  Formulation  in  which  A  and  B 

are  Unbounded  70 

3.2.2  Semigroup  representation  of  Linear 

Hereditary  Differential  Systems  74 

3.2.3  Hyperbolic  Systems  with  Inherent  Damping  77 

3.3  Robustness  Characterization  for  the  Case 

Dim  U  is  Finite  79 

3.3.1  Preliminaries  79 

3.3.2  Nyquist  Theorem  for  Distributed  Systems  81 

3.3.3  Robustness  Theorems  84 

3.3.4  Robustness  Properties  of  Distributed 

LQ  Systems  89 

3.4  Lyapunov  Control  Systems  96 

CHAPTER  4:  SYNTHESIS  OF  A  ROBUST  CONTROLLER  FOR  THE 

SUPPRESSION  OF  VIBRATION  OF  A  FLEXIBLE  BEAM  103 

4.1  Introduction  103 

4.2  A  Simply  Supported  Flexible  Beam  105 

4.3  Open-Loop  Analysis  108 

4.4  Optimal  Lyapunov  Control  112 


4.5  Implementation  of  a  Closed-Loop  Stable  Suboptimal 
Design 


123 


7 


4.6  Relationship  between  Sample  Density  and  Damping 

Ratio  129 

4.7  Closed-Loop  Stability  of  Systems  with  Constant 

Inherent  Damping  Coefficient  132 

4.7.1  Optimal  and  Suboptlmal  Control  Laws  133 

4.7.2  Suboptlmal  Implementation  136 

CHAPTER  5:  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH  145 

5.1  Conclusions  145 

5.2  Account  of  Results  In  Chapters  2-4  148 

5.3  Suggestions  for  Future  Research  150 

APPENDIX  A:  Derivation  of  Identities  1  -  5  for  LQHD 

Systems  154 

APPENDIX  B:  Robustness  Theory  for  Infinite  Dimensional 

Control  Space  160 

APPENDIX  C:  Modal  Expansion  of  Flexible  Beam  184 

APPENDIX  D:  Green's  Function  Characterization  of  Kg  186 

APPENDIX  E:  Optimal  Loop  Function  Analysis  188 

APPENDIX  F:  Derivation  of  Suboptlmal  Loop  Function  194 

APPENDIX  G:  Analysis  of  Implementation  Error  201 

REFERENCES 


208 


8 


LIST  OF  FIGURES 

Page 


2.1  LQHD  Feedback  System  25 

2.2  Nyquist  Dr  Contour  27 

2.3  Nominal  LHDF  System  35 

2.4  Perturbed  LHDF  System  37 

2.5  Multiplicatively  Perturbed  LQHD  System  50 

2.6  Perturbed  Feedback  System  Representation  for 

the  Derivation  of  LQHD  Robustenss  Properties  52 

2.7  Single-Input/Single-Output  LQHD  Example  57 

2.8  Return  Difference  of  Optimal  and  Suboptimal  Control 

for  Example  2  61 

3.1  LQ  Feedback  System  69 

3.2  Nominal  Feedback  System  82 

3.3  Perturbed  Feedback  System  85 

3.4  Multiplicatively  Perturbed  Feedback  System  91 

3.5  Perturbed  Feedback  System  Representation  for 

the  Derivation  of  LQ  Robustness  Properties  92 

4.1  Flexible  Beam  Example  106 

4.2  Optimal  Loop  Function  124 

4.3  Suboptimal  Control  Implementation.  Constant  Damping 

Ratio  Case  127 

4.4  Relationship  between  Damping  Ratio  and  Sample  Density  131 

4.5  Envelope  of  Peaks  for  Suboptimal  Loop  Function, 

Constant  Damping  Coefficient  Case  139 

4.6  Loop  Function  for  Colocated  Actuator-Sensor 

Pair,  Constant  Damping  Coefficient  Case  141 

4.7  Loop  Function  for  Colocated  Actuator-Sensor 

Pair,  Constant  Damping  Ratio  Case  142 


4.8  Suboptlmal  Control  Implementation,  Constant  Damping 
Coefficient  Case 

B.l  Illustration  of  Singular  Values  of  Linear  Bounded 
Operator  A 

B.2  Nominal  Unity  Feedback  System 

B.3  Unity  Feedback  System  Subject  to  Additive  Perturbation 
B.4  Unity  Feedback  System 


10 


CHAPTER  1 
INTRODUCTION 


1.1  An  Overview  and  Motivation 

This  thesis  deals  with  the  robustness  of  stability  of 
distributed,  linear- time-invariant  (DLTI)  feedback  control  systems. 
The  main  goal  is  to  formulate  a  practical  method  for  evaluating 
feedback  designs  based  on  the  actual  DLTI  system  characteristics. 

As  a  result,  a  design  procedure  can  be  developed  for  DLTI  systems  to 
synthesize  feedback  controllers  that  are  guaranteed  to  be 
closed-loop  stable. 

This  research  is  motivated  by  the  fact  that  although  infinite 
dimensional  systems  are  commonly  found  in  practice  and  it  is  widely 
acknowledged  [60]  that  robustness  considerations  are  of  utmost 
importance  in  the  implementation  of  feedback  compensators  for 
infinite  dimensional  systems,  a  formal  treatment  of  the  robustness 
issue  has  not  been  found  in  literature.  In  our  research,  we  have 
developed  a  robustness  characterization  for  DLTI  systems,  and  have 
shown  that  linear  quadratic  (LQ)  optimal  control  systems  have  nice 
robustness  properties  and  can  serve  as  good  reference  designs  for 
the  actual  implementation  of  the  feedback  controller.  The 
"goodness"  of  this  reference  is  measured  by  the  robustness  of  the 
optimal  system.  We  stress  in  this  thesis  that  the  optimal 
controllers  are  in  general  infinite  dimensional,  and  the  control 
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designer  must  achieve  a  design  that  is  finite  dimensional  and 
physically  realizable  [1].  We  use  a  flexible  beam  vibration 
suppression  example  to  highlight  this  discussion. 

Often  in  current  control  design  practice  for  infinite 
dimensional  systems,  a  reduced-order  model  (e.g.  [57],  [58],  [59]) 
is  first  generated  to  approximate  the  system,  and  the  control  design 
is  synthesized  using  this  reduced-order  model.  These  reduced-order 
design  methods,  however,  do  not  always  produce  feedback  control 
designs  that  are  robust,  i.e.,  remain  stable  in  the  presence  of 
uncertainties.  In  fact,  the  resulting  controllers  are  often 
unstable  when  evaluated  with  the  true  infinite  dimensional  plant. 

Another  shortcoming  of  reduced-order  design  methods  is  that 
they  are  all  problem-dependent.  There  is  no  simple  a  priori  rule  to 
determine  which  model  reduction  method  will  produce  the  best  result, 
or  what  degree  of  approximation  is  adequate.  In  particular,  one  is 
required  to  perform  the  model  reduction  on  the  open-loop  system 
(such  as  selecting  open-loop  poles),  but  the  validation  of  the 
procedure  has  to  be  checked  on  the  resulting  closed-loop  design 
(where  the  poles  have  changed).  Modeling  errors  resulting  from  the 
approximation,  such  as  parameter  variations,  nonlinear! ties, 
neglected  disturbances,  and  control  and  observation  spillovers  [2], 
[56]  can  all  cause  the  system  to  become  unstable  [3],  [4].  If  they 
do,  it  becomes  necessary  to  iterate  the  entire  approach.  A  fresh 
approximate  model  is  selected  and  the  complete  design  and  analysis 
process  is  repeated.  In  addition,  the  computational  requirements 
may  become  stringent  when  a  larger-order  model  is  used. 


It  is  hoped  that  by  using  the  optimal  control  as  a  reference 
design,  the  control  designer  can  avoid  the  aforementioned 
limitations  of  the  reduced-order  design  method.  That  is,  they  can 
start  with  an  infinite-dimensional  structural  model,  develop  an 
exact  optimal  solution,  which  is  robust,  and  use  the  robustness 
properties  of  the  optimal  solution  to  facilitate  the  actual  (finite 
dimensional)  Implementation. 

We  take  a  control  design  approach  that  has  found  a  great  deal 
of  success  in  dealing  with  finite  dimensional  systems  [5],  [6]. 

This  approach  integrates  the  multivariable  robustness  theory  with 
the  LQ  control  theory.  The  crucial  difference  between  the  finite 
dimensional  study  and  ours  is  that  in  general  one  cannot  use  the 
optimal  design  directly  in  the  infinite  dimensional  system.  Usually 
a  suboptimal  scheme  is  required  for  the  implementation.  This 
additional  design  consideration  underscores  the  importance  of  a 
robustness  study  for  infinite  dimensional  systems.  It  motivated  us 
to  investigate  the  significance  of  the  LQ  robustness  results  in  the 
synthesis  of  a  closed-loop  stable  and  implementable  suboptimal 
control  design. 

LQ  optimal  control  [54],  [55]  has  long  been  a  popular  method 
with  feedback  control  design  engineers  due  to  its  conceptual 
simplicity  and  ease  of  computation.  However  it  was  recognized  just 
over  ten  years  ago  that  optimal  control  is  merely  a  convenient  tool 
to  synthesize  a  design.  In  order  for  a  control  design  to  be  called 
truly  successful,  one  must  do  a  commendable  job  in  trading  off  the 
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system  design  performance,  robustness,  sensitivity  and  other  system 
characteristics.  These  trade-offs  are  easier  to  accomplish  in  the 
frequency  domain.  In  recent  years  we  have  witnessed  the  development 
of  many  practical  robustness  results  for  multivariable  finite 
dimensional  systems,  Including  the  guaranteed  stability  margins  of 
LQ  systems  [7],  [8],  [9],  the  introduction  of  the  singular  values  of 
the  return  difference  of  a  feedback  system  as  a  robustness  measure 
[10],  [5],  and  the  loop  transfer  recovery  technique  [3],  [6]. 
Consequently,  a  much  more  unified  treatment  for  multi-input-multi¬ 
output  control  design  methods  has  been  made  possible.  These  methods 
generally  fall  into  two  categories:  (1)  the  full-state  feedback 
type  [11],  [12],  [13];  and  (2)  the  observer  type,  i.e.  the 
linear-quadratic-guassian/loop-transfer-recovery  approach  [6],  [7], 
[14].  Although  this  theory  is  by  no  means  complete,  much  success  in 
actual  applications  and  feasibility  studies  have  been  reported  [15], 

[16].  [17].  [18]. 

We  believe  that  by  following  a  similar  development,  we  have 
obtained  a  framework  for  dealing  with  control  designs  for  infinite 
dimensional  systems  that  is  equally  effective  as  that  for  finite 
dimensional  systems. 

1.2  A  Summary  of  the  Results 

Two  important  examples  of  infinite  dimensional  systems  are 
systems  with  time  delays  and  large  flexible  structures.  They  are 
used  in  this  thesis  to  develop  the  concepts  behind  our  research.  In 


particular.  Instead  of  beginning  with  an  abstract  framework  of 
describing  infinite  dimensional  systems  using  semigroups,  we  start 
by  considering  a  general  class  of  time-delay  systems  called  the 
linear  hereditary  differential  system  [19].  We  show  that  by  using 
only  elementary  complex  function  theory,  a  very  general 
multivariable  robustness  characterization,  using  singular  values  of 
the  return  difference,  of  this  class  of  time-delay  systems  can  be 
obtained.  Based  on  this  characterization,  we  present  an  extension 
of  the  well  known  [7]  [8]  robustness  results  of  finite  dimensional 
LQ  control  systems. 

It  has  been  shown  in  one  form  or  another  [20],  [21],  [7]  that 
the  LQ  approach,  when  applied  to  finite  dimensional  systems: 

(1)  produces  a  closed-loop  system  that  satisfies  an  optimal 
frequency  domain  condition  which  is  commonly  referred  to  as  the 
Kalman  frequency  domain  inequality; 

o 

(2)  guarantees  that  the  closed-loop  control  system  has  at  least  60 
phase  margins  and  infinite  gain  margins  simultaneously  in  all 
control  input  channels  when  a  diagonal  control  weighting  matrix 
is  used. 

The  robustness  properties  in  (2)  are  direct  consequences  of  property 
(1).  In  this  thesis,  both  of  these  results  (1)  and  (2)  are  extended 
to  include  all  linear  hereditary  differential  systems.  The  results 
are  derived  directly  in  the  frequency  domain. 

A  natural  question  to  raise  here  is  "how  generally  can  one 
extend  the  results  (1)  and  (2)  above?"  We  provide  an  answer  to  this 


question  by  combining  the  Yakubovich  frequency  domain  theorem  [22], 
[23],  [26],  which  is  the  most  general  (infinite  dimensional)  version 
of  the  Kalman  frequency  domain  equality  known  to  date,  and  Desoer 
and  Wang’s  generalized  Nyqulst  thoerem  [24]  to  establish  that  for 
linear  systems  described  as  a  semigroup  over  a  Hilbert  space,  (1) 
and  (2)  above  are  true  if 

(a)  the  control  space  is  finite  dimensional;  and 

(b)  the  loop  transfer  function  belongs  to  a  class  of  transfer 
functions  which  have  finitely  many  poles  in  the  closed 
right-half  plane  (CRHP). 

Their  result  allows  us  to  generalize  the  singular  value 

robustness  characterization  to  DLTI  systems. 

In  this  thesis  we  also  study  the  infinite  dimensional  Lyapunov 

control  approach.  The  finite  dimensional  Lyapunov  control  was 

investigated  by  Wong  [12],  and  later  by  Lehiomakl  [5].  This  kind  of 

design  is  a  subclass  of  the  LQ  systems.  It  has  the  additional 

property  that  the  optimal  control  system  is  guaranteed  to  have  at 
o 

least  90  phase  margins  and  100%  gain  reduction  margins.  These 
superior  stability  margins  are  especially  useful  because  they 
provide  protection  against  Instability  due  to  actuator  failures. 

We  employ  a  flexible  beam,  simply  supported  at  both  ends,  to 
illustrate  the  use  of  the  Lyapunov  control  methodology  to  synthesize 
a  closed-loop  stable  control  design  for  the  vibration  suppression 
problem.  We  demonstrate  by  using  our  robustness  results  that  the 
inherent  damping  in  the  large  flexible  structure  plays  a  vital  role 
in  determining  whether  or  not  a  physically  realizable  closed-loop 


stable  suboptimal  (finite-dimensional)  controller  can  be  achieved  by 
spatial  sampling.  Specifically  we  show  that  we  can  design  such  a 
controller  if  the  modes  of  the  structure  has  a  constant  damping 
ratio.  The  number  of  spatial  samples  required  to  stabilize  the 
system  increases  as  this  damping  ratio  decreases.  If  instead,  the 
damping  coefficient  is  constant,  then  a  dynamical  compensator  will 
be  required  in  addition  to  Increasing  the  number  of  spatial  samples 
to  achieve  stability. 

Our  approach  focuses  on  the  loop  characteristics.  Quite  often 
researchers  (e.g.  [27])  in  the  infinite  dimensional  system  control 
area  equate  the  implementability  of  the  controller  with  whether  the 
controller  can  be  represented  as  a  compact  operator.  This  is  not 
adequate  because  in  order  to  realize  the  basis  for  that  operator,  a 
certain  basis-building  device  must  be  implemented.  This  device  may 
be  unrealizable.  Indeed  it  may  actually  be  noncausal.  For  example, 
in  large  flexible  structures  that  use  modal  feedback  control  [28] 
even  though  it  can  be  shown  that  a  compact  optimal  solution  (for  the 
modal  basis)  exists  and  the  optimal  system  is  closed-loop  robust,  a 
closed-loop  stable  suboptimal  implementation  might  not  be  possible 
because  ideal  band-pass  filters  are  not  available  in  practice. 

1.3  Organization  of  this  Thesis 

The  rest  of  this  thesis  is  organized  as  follows. 

Chapter  2  contains  all  results  on  the  linear  hereditary 
differential  systems.  We  present  a  general  robustness 
characterization,  using  singular  values  for  this  whole  class  of 


systems.  We  show  that  these  systems  have  a  finite  number  of 
open-loop  poles  In  the  CRHP.  When  restricted  to  the  LQ  control,  the 
optimal  system  satisfies  a  Kalman  frequency  domain  Inequality. 

Hence  a  set  of  guaranteed  stability  margins  can  be  readily  derived 
for  the  feedback  system.  All  results  in  this  chapter  contains  the 
known  [5]  finite  dimensional  system  results  as  a  special  case. 

Unlike  the  usual  abstract  treatment  for  linear  hereditary 
differential  systems  (e.g.  space  and  semigroup  treatment  in  [29], 
[30]).  we  are  able  to  obtain  all  results  in  this  chapter  using  only 
elementary  complex  function  theory. 

Chapter  3  contains  a  generalization  of  the  result  in  Chapter  2 
to  include  a  much  wider  class  of  DLTI  systems,  which  have  a 
finite-dimensional  control  space  and  a  finite  number  of  closed-loop 
poles  in  the  CRHP.  These  results  are  obtained  by  formalizing  the 
concepts  behind  the  Yakubovich  frequency  domain  theorem  and  the 
Desoer  and  Wang  generalized  Nyquist  theorem.  The  Lyapunov  control 
problem  along  with  the  associated  robustness  properties  are 
introduced.  As  an  illustration  we  study  the  Lyapunov  control 
problem  for  the  linear  hereditary  differential  system.  Also  we  show 
that  one  can  use  the  Yakubovich  theorem  to  rederive  the  Kalman 
inequality  in  Chapter  2. 

In  Chapter  4  we  present  an  extensive  discussion  on  the  use  of 
the  Lyapunov  control  to  suppress  the  vibration  of  a  simply 
supported,  uniform,  flexible  beam.  Since  the  optimal  solution 
requires  the  use  of  perfect  measurements,  all  along  the  beam. 


subop timal  controls  using  spatial  sampling  are  analyzed.  Also  we 
use  this  opportunity  to  study,  using  classical  control  arguments, 
the  importance  of  inherent  damping  in  a  flexible  structure.  It  is 
shown  that  if  the  inherent  damping  is  sufficiently  large,  then  the 
guaranteed  stability  margins  of  the  optimal  system  allow  us  to 
synthesize  a  closed-loop  stable  suboptimal  control. 

Chapter  5  contains  the  conclusions  and  some  suggestions  for 


future  research. 


CHAPTER  2 


ROBUSTNESS  PROPERTIES  OF  LINEAR 
HEREDITARY  DIFFERENTIAL  SYSTEMS 


2.1  Introduction 

In  this  chapter  we  study  the  robustness  properties  of  a  class 
of  time  delay  systems  whose  dynamics  are  described  by  a  linear 
hereditary  differential  equation  [19].  We  present  a 
characterization  of  robustness  for  this  class  of  systems  under 
feedback  control.  Also  we  demonstrate  that  an  optimal  frequency 
domain  condition  of  the  Kalman  inequality  type  [20],  [7]  can  be 
derived  for  linear  quadratic  control  with  infinite  horizon  [31]. 

This  result  provides  us  with  a  uniform  lower  bound  on  the  minimum 
singular  value  of  the  optimal  return  difference  matrix.  This  bound 
allows  us  to  establish  some  guaranteed  robustness  properties  for  the 
optimal  system  in  an  efficient  manner. 

In  their  paper,  Lehtomakl,  Sandell  and  Athans  [7]  considered 
the  LQ  controller  formulation  for  finite  dimensional  linear  time 
invariant  systems,  and  used  the  finite  dimensional  version  of  the 
multivariable  Kalman  inequality  to  study  the  robustness  of  the  LQ 
feedback  system.  All  results  presented  in  this  chapter  contain  the 
finite  dimensional  system  as  a  special  case. 

This  chapter  is  organized  as  follows.  In  Section  2.2  we  derive 


the  optimal  frequency  domain  condition  for  the  linear  quadratic 


hereditary  differential  (LQHD)  system.  The  optimal  condition  yields 

an  inequality  which  we  refer  to  as  the  Kalman  frequency  domain 

inequality  (KFDI).  In  Section  2.3  we  characterize  the  robustness  of 

linear  hereditary  differential  feedback  (LHDF)  systems.  This 

characterization  is  based  on  a  generalized  Nyquist  theorem  for  LHDF 

systems.  In  Section  2.4  we  combine  the  results  in  Sections  2.2  and 

2.3  to  show  that  LQHD  systems  have  good  robustness  properties.  In 

o 

particular,  these  systems  have  guaranteed  60  phase  margins, 
infinite  gain  margins  and  50%  reduction  gain  margins  simultaneously 
in  all  control  channels.  Section  2.5  contains  a  discussion  about 
the  robustness  results,  highlighted  by  two  numerical  examples. 

2.2  Linear  Quadratic  Hereditary  Differential  Systems  and  the 
Multivariable  Kalman  Frequency  Domain  Inequality 
We  consider  a  general  class  of  linear  hereditary  differential 
systems.  This  class  of  systems  includes  systems  with  time  delays, 
of  the  retarded  type  [32].  It  is  well  known  that  these  systems 
require  an  infinite  dimensional  state  space  realization,  i.e.  they 
have  infinitely  many  modes.  Multiple  discrete  time  delays  as  well 
as  distributed  time  delays  are  allowed  in  this  formulation. 

A  number  of  researchers  have  investigated  the  LQ  problem  for 
the  linear  hereditary  differential  system.  They  include 
Gibson  [30],  Delfour  and  Mitter  [19],  Delfour,  McCalla  and 
Mitter  [31],  Kwong  and  Willsky  [33],  [39]  and  Kwong  [34], 

Vinter  [35].  In  their  research,  the  optimal  control  solution  as 


well  as  its  closed-loop  stability  were  established.  However,  none 
of  these  researchers  have  addressed  the  robustness  problem. 


We  organize  this  section  as  follows.  In  subsection  2.2.1,  we 
review  the  hereditary  differential  system  and  the  solution  to  its  LQ 
problem  with  infinite  horizon.  In  subsection  2.2.2  we  state  the 
Kalman  frequency  domain  inequality  for  this  class  of  linear 
quadratic  hereditary  differential  systems.  We  present  the  proof  of 
this  result  in  subsection  2.2.3. 

2.2.1  Linear  Quadratic  Hereditary  Differential  System 

We  define  the  LQHD  system  by  following  [31].  Let  N  >  1  be  an 
integer  and  r  >  0,  and  also 


-t  =  eN  <  .  <  <  0 


Then  we  consider  the  following  autonomous  hereditary  differential 
system 


37  (t)  =  Ax(t)  + 


l 


xCt+Gj) 

0 


.  t+e^o 

,  otherwise 


+ 


x(t+e)  ,  t+0>o 

0  ,  otherwise 


(2.1) 
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x(a)  =  h(a),  a  e  [-t.O],  is  known, 


(2.2) 


where  x(t)£l^,  and  u(t)£Rm,  and  A,  A^,  Aq(0)  are  matrices  of 
dimension  nxn,  and  B  is  an  nxm  matrix;  Aq(0)  is  bounded  for 
0e[— r.O].  Also  we  assume  that  the  system  is  stabllizable,  i.e. 
there  exist  matrices  CQ  and  (^(0).  0  £  [-t,  0]  such  that  the  control 
law 


u(t)  =  CQx(t)  + 


Cj(0)x(t+0)d0 


is  closed-loop  stable. 

The  quadratic  cost  criterion  is  given  by 


J(u)  = 


[x'(t)Qx(t)+u'(t)Ru(t)]dt 


(2.3) 


where  the  weighting  matrices  Q  and  R  are  symmetric  and  positive 
definite.  The  system  of  equations  from  Eq.  (2.1)  to  (2.3)  is  called 
the  LQHD  system. 

The  optimal  control  for  the  LQHD  system  is  given  in  the 
feedback  compensation  form  [31]  by 


u(t)  =  -R-1B'K, 


0 


,x(t) 


-£ 


R"1B,K1(0)x(t+0)d0 


(2.4) 
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where  Kq,  Kj(0)  satisfy  the  following  set  of  equations, 


In  < 


A’Kq  +  KqBR  ‘B'Kq  +  Q  +  Kj(0)  +  Kj(0)  =  0  (2.5) 


K  ’  =  K0  >  0 


(2.6) 


N-l 

1 

i=l 


59^(6)  =  [A1-KoBR"1R']K1(0)  +  l  K^fife-Bj)  +  KQAo(0)  +  lyo.G); 

(2.7) 


Ki(-t)  =  KqAn 


where  6(0-0^  is  the  Dirac  delta  function  at  0=0^  and 


K2(O,0)  =  k|(-0-t)An  -  J8  K j (-0+a)BR_1B ' Kj (a)da 


N-l  r  A.K, (0+0. )  ,  -T<0+0. 

♦  5  J 


i=l 


,  otherwise 


N-l  |  K1(-0+01)Ai  .  0<01 


i=l 


,  otherwise 


j«  |  VfJKjff+O)  ,  f >-0-t 


-T 


,  otherwise 


} 


df 


(2.8) 


raj 


24 


+ 


Kj(f-e)A0(f) 

0 


f>e-T 

otherwise 


(2.9) 


The  optimal  LQHD  system  is  closed-loop  stable  [31].  Thus  it  makes 
sense  to  talk  about  the  robustness  of  LQHD  systems. 

Before  we  turn  to  give  a  statement  of  the  optimality  condition 
in  the  frequency  domain,  we  remark  that  by  following  [31]  we  have 
assumed  Q>0  in  the  system.  Partial  work  in  relaxing  this  assumption 
for  the  LQ  problem  to  the  observability  condition  is  given  in  [34]. 
The  results  in  this  thesis  can  be  adapted  to  such  cases  with  only 
minimal  changes. 


2.2.2  Kalman  Inequality  for  LQHD  Systems  -  Statement 

The  LQHL'  system  is  depicted  in  its  feedback  form  in  Figure  2.1. 
In  the  figure,  F(s)  and  H(s)  denote  the  open-loop  and  feedback 
compensator  transfer  function  matrices,  respectively.  From 
Eq.  (2.1),  one  obtains 


N 

F(s)  =  [si  -  A  -  l 
i=l 


*i  -  A0(s)]_1B 


(2.10) 


where 


A0(s)  =  J°  AQ(a)es<7da  . 


(2.11) 


Similarly,  from  Eq.  (2.4), 


H(s)  =  R^B’Kq  +  R’Vk^s) 


(2.12) 


where 


Kj(  s) 


K1(0)eaede  . 


(2.13) 


Therefore,  the  loop  transfer  function  G(s),  with  the  loop  broken  at 
the  control  channel,  i.e.,  the  point  X  in  Figure  2.1  is 

N  s0 

G(s)  =  H(s)F(s)  =R_1B'[K0  +  K1(s)][sl  -A-  J  e  \  -  A0(s)]_1B. 

i=l 

(2.14) 

Figure  2.2  illustrates  the  Nyquist  Contour  D  of  radius  r.  It 

r 

is  a  closed  contour  constructed  such  that  it  encloses  all  the 
open-loop  poles  of  F(s)  in  the  closed  right-half  plane.  Indentation 
of  radius  1/r  is  used  to  avoid  any  such  pole  on  the  imaginary  axis. 
We  shall  see  in  the  next  section  that  F(s)  has  only  finitely  many 
poles  within  the  Nyquist  contour. 

Now  we  can  state  the  Kalman  frequency  domain  inequality  (KFDI) 
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as  the  following  theorem. 

Theorem  2. 1 
Let 


0  =  {s  e  lim  D  :  Re(s)<0)  (2.15) 

r-»eo  r 

Then 

[I+G(s)]*  R[I+G(s)]  >  R.  for  all  set!  (2.16) 

where  I+G(s)  is  the  return  difference  transfer  function  matrix,  and 

K 

[I+G(s)]  denotes  its  complex  conjugate  transpose,  and  R  is  the 
control  weighting  matrix. 

Proof  of  Theorem  2.1  is  presented  in  the  next  subsection.  In 
Section  2.4  we  shall  discuss  some  consequences  of  this  theorem. 

2.2.3  Kalman  Inequality  for  LQHD  Systems  -  Derivation 

First  let  us  compute  the  left-hand-side  of  Eq.  (2.16)  using  Eq. 
(2.14). 

[  I+G(  s  )  ]*R[  I+G(  s  )  ]  =  R+B'CKq+K^s^A^sJB  (2.17) 

i* 

+  B-A"1  (s)[K0+Kj(s)]B 
|| 

+  B’A-1  (sJCKQ+KjfsJjBR^^'CKQ+K^sJ^^sJB 


where  we  have  let 


A(s) 


=  sI-A- 


r 

z  e  Ws) 

1=1 


(2.18) 


Comparing  this  with  Eq.  (2.16).  we  see  that  it  suffices  to  show 

AM(s)[K0+K1(s)HK0+Kj(s)]A(s) 

+  [Kq+K^sJIBR'^'CKq+K^s)]  >  0  (2.19) 

For  the  case  in  which  the  Nyquist  contour  does  not  coincide 
with  the  imaginary  axis,  it  is  somewhat  tedious  to  carry  out  all  the 
algebra  for  the  rest  of  this  proof  here.  Since  there  is  no  loss  in 
the  central  idea  of  the  proof,  for  the  sake  of  simplicity,  from 
hereon,  we  shall  assume  that  F(s)  has  no  open-loop  poles  on  the 
imaginary  axis,  i.e.  fi  =  imaginary  axis.  We  remark,  however,  that 
the  derivation  presented  here  can  be  adapted  to  the  general  case  in 
a  straightforward  fashion. 

We  establish  Eq.  (2.19)  by  deriving  the  following  identities: 

Identity  1 


sK^s)  =  K1(0)-K0ANeST 


es0  K^ejde 


(2.20) 


Identl tv  2 


r 


-T 


e’6  ae  Kl(8)d9  =  [A'-KoBR^B-^fs) 


N-l 


S0  . 


*  I  KoAie  1+KoAo(s> 


i=l 


+  e_STK1(-s)A^ 


-  *(s) 


+  r^sj+r^s) 


+  H^sJ+HgCs) 


where 


*(s)  =  J°  es9  J6  K j (-0+ct)Epv-1B 'Kj (a)dad0 


N-l 


>(s)  =  2  AiKi<e+ei>e 


S0 


d0 


1=1 


N-l 


r2{s)  =  I  J°  K^-O+OjJAje^d© 


1=1  '6i 


(2.21) 


(2.22) 


(2.23) 


(2.24) 
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Hjts)  =  fX  Ao(f)Kj(f+0)dfe 


S0 


d0 


(2.25) 


h2(s)  =  J°  J°  K0(f-e)A0(ndCe 


S0 


d0 


(2.26) 


Identity  3 


-1, 


*(s)  +  #*(-s)  =  K^-sJBR  *Kj( s) 


(2.27) 


Identity  4 


N-l 


—s0 . 


r1(s)+r^(-s)  =  l  Ajx^sje  1 


(2.28) 


1=1 


Identity  5 


— 1  ( s  )  +  H2(-S)  =  A^-sJK^s) 


(2.29) 


We  refer  the  readers  to  Appendix  A  for  the  proofs  of  Identities  1-5. 
Combining  Identities  1-5,  we  see  that 


st  .  _  — sr 


sKj(s)+(-s)Kj(-s)  =  KjfOJ+Kj  (O)A^e  -A^KQe 


-  { [A 1  -K0BR' [B  ’  ]Kj  ( s  )+Kj  (-s ) [A-BR_1B ‘ K-] } 


/ 

At 


>2? 

% 

s\s 
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"-1  se  .  -a6 

-  <  l  V  ♦  2  W  ) 

i=l  i=l 


-  [K0A0(s)+A0(-s)K0] 


+  Kj(-s)BR*’1B‘K1(s) 


N-l 


N-l 


-s0i 

% 

i=l  1=1 


s0. 


-  [  }  ^(sje  1  +  2  Ki(-s)Aie  ] 


-  [AqC-sJK^sJ+KjC-sJA^s)]  (2.30) 


After  some  further  manipulations,  Eq.  (2.30)  can  be  rewritten  as 


-In., 


A1(-s)[K0+K1(s)>[K0+K1(-s)]A(s)+[K0+K1(-s)]BR  *B‘ [K^K^s)] 


—  In '  i 


=  “A '  Kq-KqA+KqBR  B'Kq-Kj(0)-Kj(0) 


(2.31a) 


But  by  Eq.  (2.5), 


In  > 


-A’Kq-KqA+KqBR  ‘B'Kq  -  Kj(0)-K‘ (0)  =  Q  >  0  (2.31b) 


The  inequality  of  Eq.  (2.19)  is  established  by  combining  Eqs. 
(2.31a)  and  (2.31b).  Hence  Theorem  2.1  is  proven. 

By  following  the  above  proof  closely,  one  can  actually  give  a 
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stronger  statement  than  the  Kalman  Inequality.  Specifically,  by 
combining  Eqs.  (2.10),  (2.14),  (2.17),  (2.31a)  and  (2.31b),  one 
obtains 


(I+G(s))*R(I+G(s))  =  F*(s)QF(s)  +  R  (2.32) 

This  identity  is  called  the  Kalman  equality.  By  letting  R  =  pi, 
where  p  is  a  positive  number,  the  singular  values  of  the  return 
difference  can  be  readily  determined  without  solving  the  Riccati 
equation.  This  is  a  very  useful  fact.  In  Section  2.5,  we  use  this 
result  to  study  two  design  examples. 

2 

We  remark  that  by  considering  degenerate  operators  and  if 
spaces  [29],  [35],  the  LQHD  system  can  be  formulated  as  a  semigroup 
[30].  This  leads  to  an  alternate  proof  for  Theorem  2.1,  which  we 
present  in  Chapter  3  as  a  special  case  of  the  more  general  semigroup 
formulation  for  the  LQ  robustness  problem. 

2.3  Robustness  Characterization  of  Linear  Hereditary  Differential 
Feedback  Systems 

In  this  section  we  develop  a  robustness  characterization  of 
LHDF  systems  by: 


(1)  showing  that  the  open-loop  dynamics  of  the  linear  hereditary 
differential  system  have  finitely  many  poles  in  the  closed 
right-half  complex  plane  (CRHP),  Re  s  >  0; 


sr 
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(2)  stating  a  version  of  the  multivariable  Nyquist  theorem  for  LHDF 
systems; 

(3)  characterizing  the  robustness  measure  by  using  the  singular 
values  of  the  return  difference  and  error  matrices. 

Steps  (2)  and  (3)  are  a  generalization  of  the  robustness 
characterization  presented  by  Lehtomaki  [5]  for  finite  dimensional 
systems. 

The  robustness  characterization  is  stated  in  terms  of  the 
nominal  and  perturbed  LHDF  systems.  The  formulation  we  used  is  very 
general,  and  Includes  all  the  LQHD  system  that  we  consider  in 
Section  2.2. 

2.3.1  Nominal  and  Perturbed  LHDF  Systems 

The  nominal  LHDF  system  is  depicted  in  Figure  2.3  where  the 
loop  transfer  function  matrix  is  assumed  to  incorporate  both  the 
open-loop  plant  dynamics  and  compensation  employed. 

Specifically,  we  assume 

G(s)  =  H(s)F(s)  (2.33) 

where  F(s)  is  the  combined  pre-compensator  and  open-loop  plant 
transfer  function  and  H(s)  represents  the  post-compensator.  It  is 
further  assumed  that  H(s)  is  analytic  in  the  CRHP  except  at  a  finite 
number  of  isolated  points,  and  F(s)  has  the  following  form 
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F(s)  =  [si  -  A  -  J  e*01  At  -  J°  A^aje^da]”^  (2, 

j  4  "T 


34) 


1=1 


where  A,  A. ,  A  (a)  are  matrices  of  dimension  nxn,  and  B  is  an  nxm 
io 


matrix;  A n(o)  is  bounded  for  a  t  [-t,  0];  N  is  a  positive  integer 


and  -r  =  0M  <  ...  <6,  <0. 
N  1 


The  perturbed  feedback  system  is  depicted  in  Figure  2-4  where 
the  loop  transfer  matrix  is  assumed  to  incorporate  both  the 
open-loop  plant  dynamics  and  compensation  employed.  In  this  case, 
however,  either  or  both  of  the  plant  and  compensator  dynamics  may  be 
perturbed  from  their  nominal  values. 

Specifically,  we  assume 


G(s)  =  H(s)F(s) 


(2.35) 


where  F(s)  is  the  perturbed  pre-compensator  and  open-loop  plant 
transfer  function  and  H(s)  represents  the  post-compensator.  It  is 

/v 

further  assumed  that  H(s)  is  analytic  in  the  CRHP  except  at  a  finite 
number  of  isolated  points,  and  F(s)  has  the  form 


F(s)  =  [si  -  A  -  J  e  ^iAi  -  j°v  A0(o)eSadCT]-1I  (2.36) 
i=l  ~T 


where  A,  A^,  Aq(o)  are  matrices  of  dimension  nxn,  and  B  is  an  nxm 

'v  /v  /v 

matrix;  Aq(o)  is  bounded  for  a  e  [-r,  0];  N  is  a  positive  integer 
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and  -t  =  0~  < .  <0,  <0. 

N  1 

Note  that  both  G(s)  and  G(s)  are  mxm  matrices. 


2.3.2  Open-Loop  Unstable  Poles  of  LHDF  Systems 

In  this  subsection  we  present  a  study  on  the  open-loop  poles  of 
the  LHDF  system  in  the  CRHP.  It  is  shown  that  for  the  nomlal  system 
described  by  Eqs.  (2.33)  and  (2.34),  the  system  has  at  most  a  finite 
number  of  open-loop  poles  in  the  CRHP. 

Various  stability  tests  for  time  delay  systems  have  been 
published  in  literature.  For  examples,  Lee  and  Dianat  [36]  used  the 
Lyapunov  approach.  Thowsen  [37]  presented  a  Routh-Hurwitz  type  of 
criterion.  Tsypkin  [38]  discussed  a  graphical  method  that 
determines  the  critical  frequency  and  critical  delay  time.  None  of 
these  methods,  however,  deals  with  distributed  time  delays.  More 
importantly,  they  do  not  seem  to  provide  a  convenient  way  to  count 
the  number  of  unstable  poles.  These  two  problems  are  resolved  by 
using  an  algebra  d  studied  by  a  number  of  researchers,  e.g.  Desoer 
and  Vidyasagar  [40],  Callier  and  Desoer  [41]. 

Let  R  and  R+  denote  the  real  line  and  the  nonnegative  real 
line,  respectively,  and  denotes  the  set 


L 


1 


H.  f  |f(t)|dt  <  »  } 
Jo 


(2.38) 


The  convolution  algebra  si  consists  of  the  elements  of  the  form 
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r  0 


t  <  0 


f(t) 


Lfa(t)+  l 

1=0 


t  >  0 


(2.39) 


where 


(1)  fa(.)  e  Lj 

(2)  =  0  and  tj  >  0  for  i  =  1 ,  2. . . . . 

(3)  f  j  e  R;  and 
00 

(4)  l  |fj|  <  » 

1=0 

A  A 

Let  f  denote  the  Laplace  transform  of  f(*).  Let  si  be  the  set 

d  =  (f  |f  e  rf}  (2.40) 

A  /V 

Associate  with  si  the  pointwise  product.  Then  si  is  a  commutative 
algebra  with  the  following  property  (Callier  and  Desoer  [41]):  If 

A  A  A 

f  e  si  and  f  is  bounded  away  from  zero  at  infinity  in  the  CRHP,  then 

A  A 

f  has  a  finite  number  of  zeros  in  the  CRHP.  The  function  f  is  said 
to  be  bounded  away  from  zero  at  infinity  in  the  CRHP  if  and  only  if 
there  exist  T)  >  0,  p  >  0  such  that  for  all  |s|  >  p  in  the  CRHP, 

|f(s)  |  >  T]. 

Now  consider  G(s)  described  by  Eqs  (2.33)  and  (2.34).  Since 


det(I  +  H(s)F(s))  =  det(I  +  F(s)H(s)) 


(2.41) 
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it  Is  easy  to  see  that 


det(I  +  G(s))  = 


♦ct<5> 


(2.42) 


where 


N 

*0<!(s)  ~  det(sI  "  A  "  J 

1=1 


J0  A^oje^do)  (2.43) 


and 


4>c^(s)  =  det(sl  -  A  - 


i=l 


£,• 


AQ(a)eS£7da+BH(  s  )  ) 


(2.44) 


The  following  theorem  is  fundamental  to  the  derivation  of  a 
Nyquist  theorem  for  LHDF  systems,  because  it  allows  us  to  count  the 
number  of  encirclements. 


Theorem  2.2 

$0g(s)  as  defined  by  Eq.  (2.43)  has  finitely  many  zeros  in  the 

CRHP. 


t 


Proof  of  Theorem  2.2 

We  consider  two  cases. 

Case  1:  Ag(a)  =  0  for  all  a  t  [-r,  0].  The  proof  of  this  case  has 
been  briefly  described  by  Callier  and  Desoer  [41].  It  is  trivial  to 

S0, 


see  that 


1 


and 


all  belong  to  si,  as  are  their 


(s+1)  1  (s+1)  ““  (s+1) 

products.  Applying  Cramer’s  rule  shows  that  p(s)  is  a  linear 

s0i  0 

combination  of  1,  s  and  e  and  their  products.  Therefore 


(s+1) 


n 


£  si. 


f.y 


Case  2:  Aq(o)  /  0.  Note  that  |  aQ(o)eSada,  a ^(a)  bounded,  is 
equivalent  to  the  Laplace  transform  of  a  bounded  function  with 

/s 

compact  support.  Hence  it  is  an  element  of  si.  The  product 
f°  so  (°  so 

aQ(a)e  da  bg(a)e  da  is  also  in  si  because  it  corresponds  to 
the  Laplace  transform  of  the  convolution  of  two  functions  which  both 


°f 


have  a  compact  support.  Finally,  multiplication 
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1  s  ©  A  A 

with  •  777  ■  or  '  also  yields  a  product  that  is  in  s4. 

Extending  this  argument  as  in  the  previous  case,  we  see  that 


(a)eSCTda 


(s+1) 


n 


s6 


Now  observe  that  e 


*  and  r.  rao<a 


SO 

)e  da  are  uniformly  bounded 


for  all  s  e  CRHP.  Therefore  as  s  tends  to  ®,  ^0^(s)  is  of  the  order 


n 


s  .  Consequently 


(s+1) 


n 


tends  to  1  as  |s|  tends  to  00  in  the  CRHP. 


Thus  according  to  the  property  of  si. 


(s+1) 


n  .  hence  ^o5(s).  has  only 


*.inMrwirmrv  «rw  imx  *.  k*  mwi 
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a  finite  number  of  zeros  in  the  CRHP.  (End  of  Proof). 

2.3.3  Multivariable  Nyquist  Theorem 

Because  of  Theorem  2.2,  a  version  of  the  multivariable  Nyquist 
theorem  can  now  be  stated  in  the  following  form.  This  version 
follows  easily  from  the  standard  application  of  the  principle  of 
arguments  of  complex  variable  theory,  e.g.  [42]  [43]. 

Theorem  2.3  (Nyquist  Theorem  for  LHDF  Systems! 

Suppose  that  G(s)  is  defined  by  Eqs.  (2.33)  and  (2.34).  Then 
the  system  of  Figure  2.3  is  closed-loop  stable,  (in  the  sense  that 
<Pc{,  (s)  Eq.(2.44)  has  no  zero  in  the  CRHP,)  if  and  only  if  both  the 
following  conditions  are  satisfied: 

(1)  det  [I  +  G(s)]  /  0  ,  for  all  s  t  Q  ,  where  fi  is  defined  by 
Eq.  (2.15); 

(2)  det[I  +  G(s)] | stQ  encircles  the  origin  p  times  in  the 
counterclockwise  sense,  where  p  denotes  the  number  of  zeros  of 
$o g(s)  (Eq-  2.43)  in  the  CRHP,  counting  multiplicities. 

2.3.4.  Robustness  Theorems  for  Hereditary  Differential  Systems 
In  this  subsection,  we  develop  theorems  that  guarantee  the 
stability  of  the  perturbed  closed-loop  system  of  Figure  2.4.  The 
development  here  closely  parallels  that  in  [5].  Theorem  2.3  allows 
us  to  derive  a  simple  test  of  robustness  for  linear  hereditary 
differential  systems.  As  in  the  well  known  robustness  theorems  for 
finite  dimensional  systems  by  Doyle  [10],  Lehtomaki ,  et  al .  [7], 


Doyle  and  Stein  [3],  and  Lehtomaki  [5],  the  notion  of  singular 
values  of  matrices  are  used  to  establish  these  robustness  tests. 

Let  G(s)  satisfy  Eqs.  (2.33)  and  (2.34),  and  let  G(s)  satisfy 
Eqs.  (2.35)  and  (2.36).  Then  define  the  error  matrix  of  the 
perturbed  system  by 


(2.45) 

Also  we  know  that 


E(s)  =  G(s)  -  G(s) 


~  KM 

det  (I  +  G(s))  =  - - 

♦of <s> 


(2.46) 


where 


♦0f(s>  =  det 


N  ~ 

(si  -  A  -  Y  eS9iAi  "  J°~  A0(a)eSCTdc7)  (2.47) 
i  =  l  ~T 


N 

1  s0  (0 

♦cg(s)  =  det  (si  -  A  -  Y  e  1ai  ~  J  • ~  AQ(a)eSCrda+BH(s))  (2.48) 

1=1  -T 


We  can  now  state  the  robustness  theorem  for  additive  modeling 


error  as  follows. 
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The  perturbed  system  of  Figure  2.4  is  closed-loop  stable  (in 
the  sense  that  $cg(s)  has  no  zeros  in  the  CRHP)  if  the  following 
conditions  hold** 


(1)  a.  $o^(s)  and  $0g(s)  have  the  same  number  of  zeros  in  the  CRHP; 

b*  lf  ♦ 0e =  °*  then  =  0; 

c.  $cg(s)  has  no  zeros  in  the  CRHP. 

(2)  omln[I  +  G( s) 3  >  ^(s)]  for  all  s  e  0. 

In  the  above.  ♦o5(s) •  ^(s).  ^(s)-  are  deflned  by 

Eqs.  (2.43),  (2.44),  (2.47),  (2.48),  respectively.  The  additive 

error  matrix  E(s)  is  defined  by  Eq.  (2.45).  The  path  0  is  defined 

by  Eq.  (2.15).  The  notations  a  .  [•]  and  a  [•]  denote  the  minimum 

min-  max- 

and  maximum  singular  values  of  the  matrix. 

This  theorem  says  that  the  size  of  the  modeling  and 
implementation  error  that  a  feedback  system  can  tolerate  without 
becoming  unstable  is  given  by  the  quantity  ^min  [I+G(s)]t  when  no 
structure  of  the  error  is  assumed.  The  proof  of  this  theorem  is  of 
the  imbedding  type  [45]. 

Proof  of  Theorem  2.4 

Let  \  e  [0,1].  It  is  well  known  (e.g.  [5]  [44])  in  singular 
value  theory  that 


VnC1  +  G(S>  +  >  "mi^1  +  G<S>]  “ 


TE(s)]  (2.49) 
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Therefore  condition  (2)  guarantees  that  I  +  G(s)  +  XE(s)  is 
nonsingular  for  all  s  e  G. 

Now  suppose  that  as  X  is  varied  continuously  from  zero  to 
unity,  the  number  of  encirclements  of 

f (X, s)  =  det  (I  +  G(s)  +  XE(s))  (2.50) 

around  the  origin  changes.  Since  f(X,s)  is  continuous  in  (X,s)  e 
[0,  1]  x  (1,  its  locus  on  the  path  0  forms  a  closed  bounded  contour 
in  the  complex  plane  for  any  X  t  [0,  1].  Therefore  the  only  way  for 
a  change  in  the  number  of  encirclements  to  occur  is  for  the  locus  of 
f(X0«s)  to  pass  through  the  origin  for  some  XQ  e  [0,  1],  This  is 
equivalent  to  requiring 

det(I  +  G(sq)  +  XqE(Sq))  =  0  for  some  (Xq.Sq)  e  [0.  1]  x  0. 

(2.51) 

However  we  have  shown  that  I  +  G(s)  +  XE(s)  cannot  be  singular  by 
assuming  condition  (2).  Hence  we  conclude  that  det(I  +  G(s)), 
corresponding  to  X  =  1,  encircles  the  origin  the  same  number  of 
times  as  det(I  +  G(s)),  corresponding  to  X  =  0,  along  the  path 
s  e  G.  By  applying  Theorem  2.3  (Nyquist  Theorem),  we  see  that  the 
perturbed  system  is  closed-loop  stable.  (End  of  Proof) 


In  his  thesis,  Lehtomaki  [5]  presented  six  different  error 
models  for  describing  the  perturbed  system.  For  each  model,  a 
corresponding  robustness  theorem  similar  to  Theorem  2.4  was  stated 
and  proved.  Due  to  our  results  in  subsection  2.3.2,  Lehtomaki *s 
robustness  theorems  of  the  perturbed  system  can  be  carried  over  to 
the  LQHD  case  without  any  modification.  Theorem  2.4  is  an  example 
of  this  generalization.  Of  the  other  five  models,  we  shall  also 
describe  the  multiplicative  error  model  below.  Readers  who  are 
interested  to  see  details  on  the  rest  of  the  error  model  types  are 
referred  to  [5.  Ch.  3]  for  details. 

A  useful  way  to  describe  the  perturbed  system  with  respect  to 
the  nominal  system  is  to  define  a  multiplicative  factor  matrix  L(s) 
by 


G(s)  =  G(s)  L(s)  (2.52) 

This  representation  generalizes  the  gain  and  phase  margin 
description  of  the  single-input-single-output  case,  and  is 
particularly  useful  in  describing  the  robustness  properties  of  the 
LQHD  system  that  we  discuss  in  the  next  section.  The  corresponding 
robustness  theorem  for  multiplicative  modeling  error  is  the 
following. 

Theorem  2.5  fRobustness  Theorem  for  Multiplicative  Error): 

The  perturbed  system  of  Figure  2.4  is  closed-loop  stable  if  the 


following  conditions  hold: 


(1)  a.  $0g(s)  and  $cg(s)  have  the  same  number  of  zeros  in  the  CRHP; 
b>  if  =  °'  then  =  0: 

e.  $cg(s)  has  no  zeros  in  the  CRHP. 

(2)  L(s)  has  no  eigenvalue  in  R+,  for  all  s  e  fl. 

(3)  amin[I  +  G(s):i  >  amax  [L_1(s)  "  ^  for  a11  s  £  (2-53) 

»V 

In  the  above  $og(s),  $cg(s)  •  $Dg(s),  ^c^(s)  are  defined  by  Eqs. 
(2.43),  (2.44),  (2.47),  (2.48),  respectively.  The  multiplicative 
error  matrix  L(s)  is  defined  by  Eq.  (2.52).  The  path  Q  is  defined 
by  Eq.  (2.15). 

Proof  of  Theorem  2.5 

Let  f(X,s):  [0.1]  x  0  ■»  C  be  the  complex-valued  function 
defined  by 

f(A.s)  =  det(I  +  G(s)[I  +  ML_1(s)  -  I)]"1)  (2.54) 

The  inverses  exist  because  of  condition  (2).  It  is  easy  to  see  that 


(i) 

f(0,s) 

=  det(I 

+  G(s)); 

(2.55) 

(ii) 

f(l.s) 

=  det(I 

+  G(s)) ; 

(2.56) 

(iii)  f(X,s)  is  continuous  in  [0,  1]  x  0. 
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Now  suppose  that  as  X  is  varied  continuously  from  zero  to 
unity,  the  number  of  encirclements  of  f(X,s)  around  the  origin 
changes.  Since  f(X,s)  is  continuous  in  (X.s)  a  [0,  1]  x  0,  its 
locus  on  the  path  ft  forms  a  closed  bounded  contour  in  the  complex 
plane  for  any  X  a  [0,  1].  Therefore  the  only  way  for  a  change  in 
the  number  of  encirclements  to  occur  is  for  the  locus  of  fJX^.s)  to 
pass  through  the  origin  for  some  Xq  a  [0,  1].  This  is  equivalent  to 
requiring 


f(X().so)  =  0  for  some  (Xq.Sq)  a  [0,  1]  x  0.  (2.57) 


Since 


I+G(s)[I+X(L-1(s)-I)]_1  =  [I+G(s)+X(L“1(s)-I)][I+X(L“1(s)-I)]“1 

(2.58) 

Eq.  (2.57)  is  the  same  as  requiring  I  +  G(Sq)  +  XQ(L  *(Sq)  -  I)  be 
singular  for  some  (Xq,  Sq)  in  [0,  1]  x  [?.  However,  condition  (3) 
guarantees  that  for  all  |x|  <  1,  and  all  s  a  0 

W1  +  G(S)  +  ML_1(S)  -  I)]  >  0  (2.59) 

/v 

contradicting  Eq.  (2.57).  Therefore  we  conclude  that  det(I  +  G(s)), 
corresponding  to  X  =  1,  encirles  the  origin  the  same  number  of  times 
as  det(I  +  G(s)),  corresponding  to  X  =  0,  along  the  path  s  a  Q.  By 


applying  Theorem  2.3  (Nyquist  theorem),  we  see  that  the  perturbed 
system  is  closed-loop  stable.  (End  of  Proof) 

In  the  next  section  we  show  that  LQHD  systems  have  many  of  the 
robustness  properties  which  have  been  known  to  be  enjoyed  by  finite 
dimensional  LQ  system.  These  robustness  properties  are  consequences 
of  Theorem  2.5  and  the  KFDI  derived  in  Section  2.2. 

2.4  Robustness  Properties  of  Linear  Quadratic  Hereditary 

Differential  Systems 

The  KFDI  proven  in  Section  2.2  for  LQHD  systems  is  the  same  as 
the  KDFI  for  LQ  systems  with  finite  dimensional  state 
representations  [7],  [21],  As  in  the  finite  dimensional  case,  the 
KFDI  has  some  very  important  consequences  in  terms  of  robustness  for 
LQHD  systems.  The  purpose  of  this  section  is  to  discuss  these 
robustness  properties  of  the  LQHD  system. 

For  the  purposes  of  investigating  the  robustness  properties  of 
the  LQHD  system,  depicted  in  Figure  2.1,  it  is  convenient  to 
represent  the  model  uncertainties  as  a  multiplicative  perturbation 
factor,  as  shown  in  Figure  2.5.  In  this  configuration,  the  matrix 
L(s)  represents  the  model  uncertainties:  of  the  system  reflected  to 
the  point  X  in  Figure  2.1  where  the  loop  is  broken,  as  a 
multiplicative  perturbation  in  the  control  channel.  The  matrix  G(s) 
is  the  nominal  system  loop  transfer  function  matrix,  given  by  Eqs. 
(2.33),  (2.34).  Nominally  L(s)  is  the  identity  matrix  and  the  LQHD 


system  is  closed-loop  stable.  However  because  of  modeling  errors 

and  implementation  constraints,  the  matrix  L(s)  's  subject  to 

changes  from  its  nominal  value.  Consequently,  the  product  G(s)L(s) 

can  be  viewed  as  the  actual  system  loop  transfer  function  matrix. 

This  closed-loop  system  is  stable  if  and  only  if  the  closed-loop 

1/2 

system  of  Figure  2.6  is  stable.  In  the  figure,  R  is  a  square 
root  of  the  control  weighting  matrix,  in  the  sense  that 


R 


(2.60) 


With  this  formulation,  we  can  state  some  of  the  robustness 
properties  of  the  LQHD  system  as  the  following  theorem. 


Theorem  2.6  (Stability  Margin  of  LQHD  Systems) 

The  multiplicatively  perturbed  LQHD  system  in  Figure  2.5  is 

/v 

closed-loop  stable,  i.e.  the  polynomial  ^c^(s)  (Eq.  (2.48))  has  no 
CRHP  zeros,  provided  that  the  following  conditions  are  satisfied: 

(1)  G(s)  is  specified  by  Eqs.  (2.10)  -  (2.14),  as  the  closed-loop 
solution  to  the  LQHD  control  problem. 

(2)  (a)  ^0g(s)  nnd  $0g(s)  ^iave  tbe  same  number  of  CRHP  zeros 

(b)  if  =  0  then  =  0 

where  ^0g(s)*  ^0g(s)’  are  defined  in  Eqs.  (2.43),  (2.47), 
respectively. 

(3)  L(s)  has  no  eigenvalue  in  R+. 
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(4)  amax(R1/2L~1(s)R~1/2  -  I)  <  1  for  all  s  e  fl. 

Proof  of  Theorem  2.6 

It  is  well  known  that  $cg(s)  (Eq.  (2.44))  has  no  zero  in  the 
CRHP  [31].  Note  that 

det[I  +  (R1/2G(s)R_1/2)(R1/2L(s)R"1/2)]  =  det[I  +  G(s)L(s)] 

(2.61) 


Hence 


det[I  +  (R1/2G(s)R_1/2)(R1/2L(s)R‘1/2)]  =  tr^ 

+ce^ 


(2.62) 


Similarly, 


det[I  +  R1/2G(s)R~1/2]  - 


Vf) 


(2.63) 


Theorem  2.1  (KFDI)  implies  that 

1/9  -1/9 

a  [I  +  R  G(S)R  ]  >  1  for  all  sen  (2.64) 


Therefore  in  conjunction  with  Theorem  2.5,  condition  (3)  guarantees 
that  ^c^(s)  has  no  zeros  in  the  CRHP.  (End  of  Proof) 
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The  following  are  corollaries  to  Theorem  2.6.  Their  proofs  are 
omitted  because  they  are  completely  identical  to  the  proofs  of  their 
counterparts  for  finite  dimensional  systems.  For  examples 
interested  readers  are  referred  to  [5,  Ch.  5]. 

Corollary  2.7 

The  multiplicatively  perturbed  LQHD  feedback  system  is  stable 
if 


amax{L~1{s)_I)  <  Ccond(R)]‘1/2  .  s  e  R 


(2.65) 


where  cond(R)  denotes  the  condition  number  of  the  matrix  R, 


a  (R) 
maxv 


cond(R)  =  ■  — 

£WR' 


(2.66) 


defined  as  the  ratio  of  its  maximum  singular  value  to  its  minimum 
singular  value. 


Corollary  2.8 

If  R  >  0  is  diagonal,  then  simultaneously  in  each  feedback 
loop,  the  LQHD  system  has  the  following  guaranteed  stability 
margins, 

1.  [1/2,  «]  gain  margin 

o 

2.  +  60  phase  margin 


In  other  words,  let  L(s)  be  diagonal  and  let  each  of  these 

diagonal  elements  be  either  a  positive  number  larger  them  or 

<4  o  o 

complex  number  of  the  form  e  ,  -  60  <  $  <  60  .  Then  the  perturbed 

system  Is  closed-loop  stable. 


Corollary  2.9 


If  R  Is  block  diagonal  of  the  form 


Rj  0 

0  fc, 


(2.67) 


and  L(s)  consists  of  crossfeed  perturbation  of  the  form 


'l  X| 

P  I 


(2.68) 


Then  the  LQHD  system  is  stable  if 


,,  ,  .amin*R2).l/2 

"JWI  <  t - (RJ3  '  S“n 

max'-  1' 


(2.69) 


Similarly,  if  L(s)  is  of  the  form 


L(s)  = 


I  0 


X(s)  I 


(2.70) 


then  the  LQHD  system  is  stable  if 
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»-.«■»  <  t;HiS>1/2  • !  £  n  <2-71> 

max'  2 ’ 

In  particular,  if  R  =  pi,  where  p  is  a  positive  scalar,  then  the 
LQHD  system  is  stable  if 

amax{X(S))  <  1  •  s  n  (2.72) 

for  both  types  of  crossfeed  perturbations  (Eqs.  (2.68)  and  (2.70)). 

2.5  Discussion  on  LQHD  Control  Designs 

In  this  section,  we  look  at  two  scalar  control  examples  to 
highlight  the  previous  results  for  LQHD  systems,  and  discuss  the 
applicability  of  some  of  the  well-known  design  techniques  for 
finite-dimensional  systems  to  the  LQHD  case. 

2.5.1  Examples  of  LQHD  Systems 

The  first  example  is  a  special  case  which  shows  that  one  can 
often  obtain  the  same  robustness  measure  as  the  no  time  delay  case 
as  long  as  one  has  perfect  knowledge  of  the  time  lag.  The  second 
example  compares  a  suboptimal  design  with  the  optimal  one. 

Examp 1 e  1 : 

Consider  the  following  single-input/single-output  system  in 
which  there  is  a  time  delay  r  in  its  actuator  dynamics. 


5? 


u 

ae'sr 

c 

y 

s+a 

s+b 

actuator 

plant 

Figure  2.7  Single-Input/Single-Output  LQHD  Example 
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Let  the  quadratic  criterion  be 


J(U)  =  I  [y2(  t )  +  U2(t)]dt  (2.73) 

J0 

Then  the  Kalman  equality  (Eq.  (2.32))  can  be  used  to  compute  the 
magnitude  of  the  optimal  return  difference,  which  is  given  by 


ll  +  gWl2*!-  l(jw4(Ju,b)|2  (2.74) 

It  is  interesting  to  note  in  this  equation  that  the  quantity 
|l+g(j(j)|,  which  is  the  robustness  measure  formulated  in  Section 
2.3,  is  independent  of  t.  The  reason  for  this  is  that  if  the  LQ 
system  has  full  knowledge  of  the  time  delay  element,  then  as  a 
consequence  of  the  full  state  feedback,  the  LQ  system  will 
compensate  so  that  the  same  stability  margin  is  obtained  as  in  the 
t  =  0  case. 


Strictly  speaking,  this  corresponds  to  a  positive  semi-definite  Q 
matrix  which  is  not  considered  in  this  chapter.  However,  this 
limitation  is  removed  in  Chapter  3. 


This  example  is  drawn  from  Gibson  [30],  who  has  computed  the 
control  gains  for  the  example  numerically.  We  consider  the 
suboptimal  design  using  those  control  gains  here. 

Let  the  dynamics  of  the  system  be 

x(t)  =  x(t)  +  x(t-l)  +  u(t)  (2.75) 

where  x(»)  and  u(«)  are  scalar  functions.  Gibson  divides  the  unit 
interval  into  N  subintervals  in  order  to  approximate  the  state 
x(t+»)  over  the  unit  interval  as  a  piecewise  constant  function. 
Gibson  then  derived  an  N-th  order  approximate  model  for  the  system. 
He  solved  this  system  numerically  for  the  control  gains  k^,  kj(»)  in 

u(t)  =  -  kgx(t)  -  kj(0)x{ t+0)d9 

The  resulting  values  are  given  for  0  =  0.0,  -0.1,  -0.2 . -1.0,  and 

N  =  17  ,  29  ,  50  ,  74,  by  the  following  table. 
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Table  2.1 
29  50 


0 

ko 

2.8260 

«- 

2.8190  2.8148 

kj(0) 

2.8130 

0.0 

0.6684 

0.6547 

0.6469 

0.6435 

-0.1 

0.7726 

0.7169 

0.7179 

0.7273 

-0.2 

0.8467 

0.8239 

0.8209 

0.8258 

-0.3 

0.9961 

0.9508 

0.9434 

0.9607 

-0.4 

1.0822 

1 . 1020 

1.0895 

1.1023 

-0.5 

1.2811 

1.2822 

1.2633 

1.2694 

-0.6 

1.5200 

1.4963 

1.4693 

1.4965 

-0.7 

1.6606 

1.7501 

1.7125 

1.7315 

-0.8 

1.9802 

2.0499 

1.9987 

2.0480 

-0.9 

2.3648 

2.4033 

2.3347 

2.3748 

-1.0 

2.5852 

2.6730 

2.7284 

2.7541 

We  consider  the  suboptimal  control  of  the  following  form. 


u(t)  =  -  kQx(t) 


10 

J=1 


-  ill 
10 


x(t+0)d0 


where  belongs  to  one  of  three  cases 


a)  i}  =kl(-  Jfi); 


•JL 


Pi'iv 


IVS'J V\  ‘ 


ft) 


W1 


■ill  9mm 
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(2)  kj  -  [k;(-  +  kj(-  fo)]/2s 

(3)  k^kjC-Jj). 

The  return  difference  for  these  suboptlmal  controllers  are 
plotted  in  Figure  2.8,  for  the  N  =  17  approximation.  Also  shown  is 
the  optimal  return  difference  which  is  computed  using  the  Kalman 
equality  (Eq.  (2.32)). 

It  is  clear  from  these  plots  that  the  suboptlmal  controller  (2) 
approximates  the  optimal  system  very  closely.  Also  by  looking  at 
the  plots  for  N  =  50,  74,  23,  we  see  that  only  insignificant 
improvement  can  be  obtained  by  using  those  higher  order  models. 

2.5.2  Extensions  of  LQHD  Control  Designs 

In  this  thesis  we  have  not  considered  the  effect  of  observers 
in  the  controller  designer.  However,  the  most  important  control 
design  that  has  evolved  from  the  multivariable  robustness  study  for 
finite  dimensional  systems  is  probably  the  linear-quadratic- 
gaussian/loop-transfer-recovery  approach  [3],  [6].  This  approach 
combines  an  observer  (usually  the  Hainan  filter  [61])  design  for  the 
system  with  the  asymptotic  properties  of  the  LQ  regulator  [62]  to 
synthesize  controllers  which  have  desirable  loop  properties.  It 
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seems  to  us  that  this  approach  is  also  valid  for  the  LQHD  system, 
provided  that  the  system  is  minimum-phase  (a  requirement  of  the 
finite  dimensional  case  also).  A  good  filter  design  to  use  is  the 
one  for  LQHD  systems  by  Kwong  and  Willsky  [39]. 

Since  our  singular-value  formulation  for  robustness  is 
identical  to  the  one  for  finite  dimensional  systems,  all  the  finite 
dimensional  results  that  depend  only  on  the  manipulation  of  the 
singular  value  decomposition  are  also  valid  for  the  LQHD  system.  In 
particular,  the  design  techniques  which  exploit  the  structural 
information  of  a  system  are  appropriate  for  LQHD  system  designs 
also.  Examples  of  these  techniques  are  the  celebrated  design 
analysis  method  using  the  error  structure  by  Lehtomakl,  et  al.  [63], 
and  the  design  adjustment  technique  using  the  control  weighting 
matrix  by  Lee,  et  al .  [13],  and  the  localizing  robustness  analysis 
by  Lee,  et  al .  [64].  Without  being  repetitious  here,  we  refer  the 
interested  readers  to  the  last  chapter  of  this  thesis  for  an 
elaboration  of  these  techniques. 


CHAPTER  3 


ROBUSTNESS  PROPERTIES  OF  DISTRIBUTED. 
LINEAR-TIME-INVARIANT  SYSTEMS 

3.1  Introduction 

In  Chapter  2  we  have  developed  a  multivariable  robustness 
characterization  for  LHDF  systems,  using  singular  values  of  the 
return  difference.  Also  we  have  shown  that  the  LQHD  control  design 

(1)  produces  a  closed-loop  system  that  satisfies  an  optimal 
frequency  domain  condition  which  is  commonly  referred  to  as  the 
Kalman  frequency  domain  inequality  (KFDI); 

O 

(2)  guarantees  that  the  closed-loop  control  system  has  at  least  60 
phase  margins  and  infinite  gain  margins  simultaneously  in  all 
control  input  channels  when  a  diagonal  control  weighting  matrix 
is  used. 

For  systems  described  by  semigroups,  the  above  results  are  also 
valid.  Yakubovich  [22]  [23]  proved  the  result  (1)  for  linear 
systems  in  which  both  the  system  and  input  operators  are  linear 
bounded,  and  in  a  separate  paper  with  Likhtarnikov  [26],  for  systems 
described  by  linear  unbounded  operators. 

In  order  to  establish  the  result  (2),  we  need  to  extend  the 
singular-value  robustness  characterization  to  the  present  case.  It 
appears  that  the  paper  by  Desoer  and  Wang  [24]  contains  the  most 
general  Nyquist  Theorem  available  for  our  purpose.  Their  result  was 
derived  by  using  a  quotient  algebra  of  transfer  functions  [41].  We 
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use  this  algebra  to  describe  the  nominal  and  perturbed  system.  Also 
we  study  the  infinite  dimensional  Lyapunov  control  design  approach. 
The  finite  dimensional  Lyapunov  control  was  investigated  by  Wong 
[12].  This  kind  of  design  is  a  subclass  of  the  LQ  system.  It  has 
the  additional  property  that  the  optimal  control  system  is 

O 

guaranteed  to  have  at  least  90  phase  margins  and  100%  gain 
reduction  margins. 

The  rest  of  this  chapter  is  organized  as  follows.  In  Section 
3.2.  we  discuss  the  Yakubovich  Frequency  Domain  Theorem  (YFDT),  also 
we  revisit  the  LQHD  system  to  show  that  the  YFDT  can  be  used  to 
derive  the  KFDI  in  Chapter  2.  In  Section  3.3  we  develop  the 
robustness  characterization  for  distributed,  linear-time-invariant 
(DLTI)  systems,  described  by  linear  operators  with  a  finite 
dimensional  control  space.  The  result  is  then  used  to  derive  the 
robustness  properties  of  LQ  optimal  systems.  In  Section  3.4,  we 
discuss  the  Lyapunov  control  system. 

In  Section  3.3  we  emphasize  that  the  robustness 
characterization  is  only  valid  when  the  control  space  is  finite 
dimensional,  even  though  the  YFDT  applies  to  the  infinite 
dimensional  case  also.  This  shortcoming  is  due  to  the  absence  of 
any  Nyquist-type  of  theorem  for  infinite  dimensional  control  space. 
It  is,  however,  possible  to  formulate  a  robustness  characterization 
based  on  the  singular  values  of  linear  operators,  for  which  the 
control  space  is  allowed  to  be  infinite  dimensional.  We  present 
this  result  in  Appendix  B. 


3.2  Yakubovich  Frequency  Domain  Theorem 


Yakubovich  [22],  [23]  proved  that  under  certain  rather  general 
assumptions,  the  KFDI-type  of  criterion  (see  Chapter  2)  is  a 
necessary  and  sufficient  condition  for  the  existence  of  LQ  optimal 
control  for  infinite  dimensional  systems.  The  necessary  part  of  his 
theorem  provides  a  basis  for  us  to  study  the  robustness  properties 
of  LQ  systems  represented  as  a  semigroup. 

Let  X  and  U  be  two  Hilbert  spaces.  Each  of  these  spaces  has  an 
inner  product.  With  an  abuse  of  notation,  let  each  inner  product  be 
denoted  by  (•,•).  and  the  norm  associated  with  the  inner  product  be 
denoted  by  |  |*|  | . 

First  we  present  the  results  for  systems  described  by  bounded 
operators.  This  covers  all  finite  dimensional  systems  and  is 
sufficient  to  treat  a  number  of  systems  described  by  linear  compact 
and  integral  operators.  The  results  for  linear  unbounded  operators 
are  given  in  subsection  3.2.1. 

Let  Q  be  a  linear  bounded,  self-adjoint  operator  on  X  such  that 

(x,  Qx)  >  0  for  all  x  t  X  (3.1) 

Let  R  be  a  linear  bounded,  self-adjoint  operator  on  U  such  that 
for  come  6  >  0, 

(u,  Ru)  >  6  |  |u | | ^  for  all  u  e  U  (3.2) 


Saties 
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Definitions:  Let  A:  X  -»  X,  B:  U  -»  X  be  linear  operators  and 
I  =  identity  operator  in  X  (or  U).  The  spectrum  sp(A)  of  A  is  the 
set  of  all  complex  numbers  X  such  that  XI-A  does  not  have  a  linear 
bounded  inverse.  A  is  called  Hurwitz  if  sp(A)  is  disjoint  from  the 
CRHP.  The  pair  A,  B  is  said  to  be  stabllizable  if  there  is  a  linear 

It 

bounded  operator  H  such  that  A  +  BH  is  Hurwitz.  Let  A  ,  H  denote 
the  adjoints  [46]  of  A,  H,  respectively. 

Theorem  3.1  (Yakubovich  Frequency  Domain  Theorem  (YFOTl) 

Let  A,  B  be  a  stabllizable  pair  of  linear  bounded  operators  and 
suppose  that  no  point  in  sp(A)  is  on  the  imaginary  axis.  Then  a 
necessary  and  sufficient  condition  for  the  existence  of  a  linear 

It 

bounded  self-adjoint  operator  K-  X  X  and  an  operator  H  :  X  -»  U 
such  that 


KA  +  A*K  +  Q  -  HRH*  =  0  (3.3) 

KB  -  HR  -  0  (3.4) 

and  | |exp[(A  -  BH*)t]xQ||  is  squared  integrable  over  t, 
i.e.  e  ^((O,®).  X),  for  all  Xq  t  X  is 

((jwl  -  A)-1Bu,  Q(jul  -  A)-1Bu)  +  (u.  Ru) 

=  ((I  +  G(ju))u.  R(I  +  G(ju))u),  for  all  u  e  U  +  jU,  u  >  0 

(3.5) 


where 


G(j(j)  =  H*( JojI  -  A)-1B  .  (3.6) 

Eq.  (3.5)  implies  the  following  domain  inequality  that  corresponds 
to  the  KFDI  in  Chapter  2, 

((I  +  G( Jw))u,  R(I  +  G(j<j))u)  >  (u.  Ru) 

for  all  u  t  U  +  JU. 

In  subsection  3.2.1,  we  state  the  corresponding  theorem  for 
linear  unbounded  operators. 

The  relationship  of  Theorem  3.1  to  the  LQ  optimal  control  is  as 
follows.  Let  A  be  the  infinitesimal  generator  of  a  strongly 
continuous  semigroup  of  operators  T(t)  over  X,  i.e. 

lim  ||T(t)x-x| |  =  0  for  all  x  e  X  (3.7) 


The  optimal  control  problem  is  defined  as  [22]: 

Find  u(»)  £  ^([O.w);  U)  to  minimize  the  quadratic  cost  functional 


(x(t),  Qx(t))dt  +  P  (u(t),  Ru(t))dt  (3.8) 

J0 
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where  Q  and  R  are  defined  in  Eqs.  (3.1)  and  (3.2),  respectively, 
subject  to  the  dynamics 


=  Ax(t)  +  Bu(t)  ;  x(0)  given 


(3.9) 


*  1/2 

Fur therf more  we  assume  that  A  ,  Q  is  a  stabilizable  pair,  where 
1/2 

Q  is  the  square  root  of  Q: 

Q  =  Q1/2  Q1/2 

The  optimal  control  solution  is  characterized  by  Theorem  3.2. 


Theorem  3.2  (Yakubovich  [22]) 

(1)  The  solutions  K,  H "  to  Eqs.  (3.3)~(3.4)  with  the  property  of 
exp[(A-BH*)t]  Xq  e  ^((O,  00  ) ,  X)  are  unique. 

(2)  The  control  given  by 

u°(t)  =  -  H*x(t)  (3.10) 

is  the  optimal  control  solution  to  Eqs.  (3.8)-(3.9). 

(3)  The  optimal  system  is  closed-loop  exponentially  stable. 

The  optimal  control  can  be  represented  in  feedback  form  in  the 

A  A 

frequency  domain  by  Figure  3.1,  where  u(s)  and  x(s)  are  the  Laplace 
transforms  of  u(t)  and  x(t),  respectively.  By  breaking  the  loop  at 

A 

u(s),  one  sees  that  the  loop  transfer  function  is 

G(s)  =  H*(sl  -  A)-1B  (3.11) 


Correspondingly  the  operator-valued  function  I  +  G(s)  is  the  return 
difference. 

Theorems  3. 1-3.2  allow  one  to  extend  the  finite  dimensional  LQ 
robustness  theorems  [7]  [8]  [20]  to  the  infinite  dimensional  case. 
This  is  possible  when  the  dimension  of  the  control  space,  i.e.  U,  is 
finite.  Section  3.3  deals  with  this  case.  A  generalized  Nyquist 
theorem  by  Desoer  and  Wang  [24]  is  used  to  extend  the  finite 
dimensional  results. 

In  order  to  characterize  the  robustness  of  distributed  systems, 
it  is  essential  that  we  consider  the  case  in  which  the  operator  A  is 
linear  unbounded,  since  the  Infinitesimal  generator  of  a  semigroup 
is  generally  unbounded. 

3.2.1  Evolution  Formulation  in  which  A  and  B  are  Unbounded 

The  case  in  which  A  and  B  are  unbounded  operators  is  considered 
in  this  subsection.  Under  a  given  set  of  conditions,  Likhtarnikov 
and  Yakubovich  [26]  showed  that  the  YFDT  (Theorem  3.1)  is  also  true 
for  this  case.  We  discuss  their  extended  result  for  the  necessary 
part  of  Theorem  3.1  here.  That  result  allows  us  to  study  the 
robustness  characterization  of  LQ  optimal  systems  in  which  A  is 
unbounded . 


The  problem  considered  in  [26]  is  an  evolution  problem  of  the 
parabolic  type.  In  subsections  3.2.2  and  3.2.3,  we  discuss  how  the 
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LQHD  system  and  the  distributed  system  described  by  a  second-order 
partial  differential  equation  can  be  adapted  to  use  this  frequency 
domain  theorem,  even  though  their  natural  descriptions  are  not 
parabolic. 

The  result  in  [26]  assumes  that  the  parabolic  evolution 
equation  of  the  system  satisfies  a  coercive  and  a  regularity 
condition.  Its  assumptions  are  similar  to  those  in  Lions  [75,  Ch. 
3],  which  contains  a  detailed  study  on  the  optimal  control  problem. 
However.  [75]  did  not  consider  the  case  in  which  B  is  unbounded. 

Let  X  and  U  be  two  Hilbert  spaces  given  as  before.  Also  let  V 
be  a  Hilbert  space  such  that  V  is  a  dense  subspace  of  X.  Let  V  be 
the  dual  space  of  V.  Then 

V  C  X  C  V 

Let  (•.•)y  and  ||«||  denote  the  inner  product  and  its 
corresponding  norm  on  V,  and  when  necessary  for  the  distinction,  we 
will  use  (•.*)x.  I  M  lx  for  X. 

Let  S£(Y,Z)  denote  the  space  of  linear  continuous  operators  from 
Y  to  Z. 

Let  A,  B  be  linear  unbounded  operators: 

A  e  tf(V.V) 

B  e  i£(U,  V  ) 


Ms 
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such  that  the  operator  A  is  an  infinitesimal  generator  of  a  strongly 
continuous  semigroup  of  linear  bounded  operators  in  X,  which 
satisfies  the  following  coercive  condition:  there  exists  real 
numbers  X,  a  such  that  a  >  0  and 

(v,-Av)v  +  (-Av,v)y  +  X|  |vj  |*  >  a  |  |v| 

for  all  v  £  V  (3.12) 

Then  A  and  B  describe  an  evolution  system  as  follows, 

-  Ax(t)  +  Bu(t)  (3.13) 

x(0)  =  XQ.  XQ  £  X. 

where  x(t)  e  X,  u(t)  £  U  and  t  >  0,  and  u(*)  £  ^([O,®0).  U) . 

Other  conditions  such  as  the  regularity  condition  for 
sufficient  "smoothness"  which  is  required  for  the  synthesis  of  the 
optimal  control  are  described  in  [26].  As  remarked  earlier,  we  do 
not  elaborate  on  these  conditions  here  since  we  are  only  interested 
in  the  necessary  part  of  the  frequency  domain  theorem  in  [26]. 
Instead  we  simply  assume  that  the  optimal  control  solution  exists 
and  state  the  frequency  domain  condition  for  the  case  in  hand.  This 
is  stated  as  a  continuation  of  theorem  3.1. 
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Theorem  3.1A  (Unbounded  Operator  Casel: 

Suppose  that  no  point  in  sp(A)  is  on  the  imaginary  axis  and 
there  exists  a  solution  K  e  tf(V  ,X)  D  5£(X,V),  self-adjoint  in  X,  for 
the  equations 

KA  +  A*K  +  Q  -  HRH*  =  0 

KB  -  HR  =  0 

Then  Eq.  (3.5)  holds,  i.e. 

(( j(oI-A)_1Bu,  Q(juI-A)_1Bu)  +(u,Ru) 

=  ((I  +  G( j(j))u,  R(I  +  G(ju))u)  for  all  u  e  U  +  jU  (3.14a) 


where 

G(  ju)  =  H*(jul  -  A)-1B 

and  Q  e  S£(V,V')  and  R  e  S£(U,  U)  are  self-adjoint. 
If  in  addition, 

(Qx,x)v  >  0  for  all  x  e  V 

and  for  some  6  >  0 

(Ru,u)  >  6||u||^  for  all  u  e  U 


i 


((I  +  G(ju))u,  R(I  +  G(ju))u)  >  (u,Ru) 


for  all  u  e  U  +  jU  (3.14b) 

In  the  next  subsection  we  show  that  the  LQHD  system  can  be 
represented  in  the  form  of  Eq.  (3.13),  and  satisfies  Eq.  (3.12).  In 
subsection  3.2.3,  we  do  the  same  for  hyperbolic  systems  with 
sufficient  Inherent  damping. 

3.2.2  Semigroup  Representation  of  Linear  Hereditary  Differential 
Systems 

Consider  the  linear  hereditary  differential  system  that  we 
define  in  Eqs.  (2.1)  -  (2.2).  This  can  be  represented  in  semigroup 
form  [21]  [29]  as  follows. 

Let  the  space 

M2  =  Rn  x  ^((O.  -t),  Rn) 

Then  Eqs.  (2.1)  -  (2.2)  can  be  described  equivalently  as 
^|=  Az(  t)  +  Bu(t) 
z(t)  =  h;  h  e  D(A) 

2 

where  A  is  a  closed  linear  operator  on  M  ,  with  a  dense  domain  in 
M2,  defined  by 
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(Ah) (0)  = 


N  rO 

Ah(0)  +  l  Aih(0i)  +  £  Ao(0)h(0)d0; 


0  =  0 


i=l 


dhf  0) 
d0 


:  0  *  0 


and 


B  = 


'  B  ' 

.0  . 


This  linear  operator  A  is  generally  unbounded.  In  order  to  use  the 
result  in  subsection  3.2.1,  we  must  show  that  A  satisfies  Eq. 

(3.12). 

Since  sp(A)  lies  entirely  to  the  left  of  the  line  Re  s  =  (3  for 
some  finite  /3  >  0  in  the  complex  plane  [31],  A  -  pi  generates  a 
contraction  semigroup.  Consequently  A  -  pi  is  dissipative  [55], 
i  .e. 

((A  -  pi)x.x )  +  (x.  (A  -  pi)x)  =  0 

Expanding,  this  is  identical  to 

(Ax,x)  +  (x.Ax)  -  2/3 1  |x ||^  =  0 

Therefore  Eq.  (3.12)  is  satisfied.  Hence  Theorem  3.1A  is 
applicable.  We  remark  that  in  fact,  it  is  possible  to  draw  the  same 
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conclusion  for  a  larger  class  of  linear  hereditary  differential 
systems  (Delfour  [76]). 

We  now  show  that  Eq.  (3.14)  implies  the  KFDI  in  Eq.  (2.16)  for 
the  LQHD  systems.  The  optimal  solution  for  the  semigroup 
formulation  gives 

C(ju)  =  H*(ja>I  -  A)-1B 

x 

where  the  operator  H  is  defined  by  [31], 

H*f  =  KQf(t)  +  I  K1(0)f(t+B)d0 
J  _T 

where  KQ  and  K^B)  are  matrices  given  by  Eqs.  (2.5)  -  (2.9). 

It  is  a  simple  algebraic  exercise  to  show  that 

(jul  -  A)  *Be^Wtv  :=  F(ju)e^utv 

Let  u  =  e^utv,  v  e  U,  in  Eq.  (3.14b).  Then 

(I  +  G(  ju)  )u  =  (eJwt  +  G{j(j)eJut)v 

where  F(s)  in  defined  in  Eq.  (2.10).  Also  by  using  the  definition 
of  H*.  we  get 

G(ju)eJutv  =  [K0F(  Ja>)eJ<ot  +  |  K^OJFf  ja>)ejw(  t+6)d0]v 

^  — T 


w 

w' 

S 

sf 


Hence  manipulating  the  above,  we  obtain 


(I  +  G(jo)))eJt,tv  =  (I  +  G(Ja)))eJutv  1 

where  G(s)  is  defined  in  Eq.  (2.14).  Since  u  and  v  are  ) 

l 

arbitrary,  I  +  G(jcj)  and  I  +  G(ju)  are  identical.  Therefore  the  jj 

t 

YFDT  in  this  section  implies  the  KFDI  in  Chapter  2.  j 

| 

3.2.3  Hyperbolic  Systems  with  Inherent  Damping  5 

i 

In  this  subsection,  we  consider  the  class  of  distributed  ] 

control  system  described  by  the  second-order  differential 

equation  [28]  j 

I 

djlil  +  L  dx|tl  +  px(t)  =  BoU(t),  t  >  0  (3.15)  j 

dt  i 

i 

} 

i 

where  x(t)  and  u(t)  belong  to  the  Hilbert  spaces  X  and  U, 
respectively. 

Let  P  be  a  self-adjoint  linear  operator  with  domain  D(P) ,  and 
assume  that  it  satisfies  the  coercive  condition:  there  exists  a  >  0 
that 

(Px,x)  >  aj|x||^  for  all  x  e  D(P)  (3.16a) 

In  addition  assume  that  P  *  is  compact. 

Let  L  be  a  nonnegative,  symmetric  linear  operator  with  domain 
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D(L)  3  D(P),  and  assume  there  exists  nr  >  0  such  that 


|  |Lx|  |  <  i |  |Px |  |  for  all  x  t  D(P) 


(3.16b) 


Rewriting  Eq.  (3.14)  in  the  form  of  the  parabolic  description 
of  Eq.  (3. 13) .  we  get 


x(t)' 

'  o  i 

x(t)' 

0  ' 

d_ 

dt 

dxCt) 

~ 

-P  -L 

dx(Q 

R 

dt  J 

•  • 

dt 

W  J 

"o 

u(t) 


Then  Theorem  3.1A  is  applicable  if  A  = 


satisfies 


l-P  -LJ 

Eq.  (3.12)  when  extended  to  the  space  D(P)  x  X.  from  D(P)  x  D(P). 

<v 

Gibson  [28]  proved  that  A  can  be  extended  to  D(P)  x  X  and  the 
extension  generates  a  contraction  semigroup  of  linear  bounded 
operators  on  this  space;  of  which  the  inner  product  is  defined  as 


>  =  (Px,.x2)x*  (yi.y2)x 


and  A  is  dissipative.  Therefore  [55], 


<  A 


x 

X 

•  • 

X 

PV 

X 

9 

>  +  < 

,  A 

y. 

y. 

y. 

y. 

>  <  0 
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Consequently  for  any  real  a  >  0,  X  =  a,  the  coercive  condition  in 
Eq.  (3.12)  is  satisfied  for  V  =  D(P )  x  D(P)  dense  in  D(P)  x  X. 

3.3  Robustness  Characterization  for  the  Case  Dim  U  is  Finite 
In  this  section  we  discuss  a  version  of  the  multivariable 
Nyquist  theorem  due  to  Desoer  and  Wang  [24]  for  a  distributed  system 
which  has  a  finite  dimensional  control  space.  We  use  this  theorem 
to  extend  the  robustness  characterization  for  LQHD  systems  in 
Chapter  2  to  more  general  distributed  systems.  Then  a  robustness 
theorem  will  be  stated  for  the  LQ  optimal  system. 

3.3.1  Preliminaries 

A 

Callier  and  Desoer  [41]  introduced  an  algebra  35  of  transfer 
functions  that  is  suitable  for  studying  distributed,  linear 

A 

time-invariant  systems.  The  algebra  35  is  formed  as  a  quotient 
algebra  [72]  of  two  convolution  algebras.  It  allows  us  to 
manipulate  distributed  systems  within  the  algebra  in  ways  similar  to 
those  for  the  finite  dimensional  system  within  the  algebra  of  proper 
rational  transfer  functions.  Desoer  and  Wang  [24]  successfully 

A 

applied  the  algebra  28  to  obtain  a  generalized  Nyquist  theorem  for 

A 

distributed  systems.  The  definition  of  28  that  we  adopt  in  this 
chapter  is  similar  to  that  in  [24]. 

Let  a  e  R,  and  let  $ i{a)  denote  the  convolution  algebra  which 
consists  of  elements  of  the  form 
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n  t)  = 


.  t  <  0 


00 


f  (t)  +  1  f i^C t-tj)  .  t  >  0 

1=0 


where 

(2) 

fa< 

(2) 

f 

Jo 

(3) 

co 

(4) 

fi 

00 

(5) 

I 

1=0 

Ua(t)!e  atdt  <  »; 


and 

-at. 


f(*)  Is  said  to  belong  to  rf(a)  if  and  only  If  there  exists  a  R, 
Oj  <  ct  such  that  f(»)  a  rf(cr ^ ) . 

A  A 

Let  f  denote  the  Laplace  transform  of  f(*).  Let  *$(ct)  be  the 

set 


rf(a)  =  {f  | f  £  «*(a)} 


and  (a)  be  the  set 

s4°°(cr)  =  { f  J f  a  d(o)  H  'J'(a)} 


S 

I 

I 


V 


S 
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A  A  A 

where  'J'(a)  is  the  set  consisting  of  all  f  such  that  f  is  bounded 
away  from  zero  at  infinity  in  the  closed  right-half  complex  plane 
Re  s  >  a,  i,e.  there  exist  17  >  0,  p  >  0  such  that  for  all 

A 

s  e  {Re  s  >  a  and  js  -  a\  >  p},  |f(s)|  >  t). 

A 

Finally  let  3(a)  be  the  algebra  with  pointwise  product. 


3(a)  .-s  {f  jf  =  n/d  with  n  t  s4(a)  and  d  e  rf  (a)} 


and  3(a)  be  the  corresponding  convolution  algebra,  i.e. 


3(a)  =  {f (•) |f  e  3(a)) 

a  A  A 

rf( 0),  rf(0),  rf(0),  d  (0),  3(0)  and  3(0)  are  abbreviated  as  d,  d,  rf, 

ACO  A 

d  ,  3  and  3,  respectively. 

3.3.2  Nyquist  Theorem  for  Distributed  Systems 

We  consider  the  feedback  system  in  Figure  3.2,  where 

G(s)  e  (3.17) 

Let  ( N,D )  be  a  right  coprime  representation  of  G(s)  [41],  [24]. 


Equivalently, 


(1) 

G(s)  =  N(s)D-1(s) ; 

(3.18) 

(2) 

H(s)  t  3“".  D(s)  t  2”™-. 

(3) 

there  exist  U(s)  e  simxm, 

V(s)  e  stfmxm  such  that 

U(  s)W(s)  +  V(s)D(s)  = 

I; 

(3.19) 

(4) 

det  D(s)  is  analytic  and 

bounded  away  from  zero  at  » 

in  the 

CRHP. 

Since  D(s)  e  rfc°mxm  c  rimxm,  under  condition  (4)  det  D(s)  has 
only  a  finite  number  of  zeros  in  the  CRHP  [41].  Let  p  denote  this 
number,  counting  multiplicities.  In  addition  we  assume  that 

(a)  G(s)  tends  to  a  constant  matrix  as  |s|  -*  ®  in  the  CRHP;  (3.20) 

(b)  G(s)  has  no  zero  on  the  imaginary  axis.  (3.21) 

Under  this  formulation,  one  can  deduce  a  Nyquist  theorem  based 

on  Desoer  and  Wang  [24]  as  follows. 


S 

! 


5 

i 
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Theorem  3.3  (Generalized  Nvouist  Theorem): 

The  system  of  Figure  3.3  is  closed-loop  stable  if  and  only  if 

(1)  det[I  +  G(ju)]  /  0  for  all  u  e  R+  (i.e.  w  >  0);  and 

(2)  det[I  +  G(j(j)]|wt  ,,  encircles  the  origin  p  times  in  the 

counterclockwise  sense. 
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3.3.3  Robustness  Theorems 

We  call  the  feedback  system  in  Figure  3.2  the  nominal  feedback 
system.  We  define  the  perturbed  feedback  system  to  be  the  system  of 
Figure  3.3,  obtained  by  replacing  G(s)  with  G(s)  in  Figure  3.2. 
Similarly  to  G(s),  we  assume 

(1)  G(s)  e  Smxm;  (3.22) 

/v 

(2)  C(s)  has  a  right  coprime  representation  (N,  D).  (3.23) 

As  in  Chapter  2.  we  define  the  additive  error  matrix  by 


E(s)  =  G(s)  -  G(s)  (3.24) 

For  the  study  of  robustness,  it  is  useful  to  write  down  the 
open-  and  closed-loop  characteristic  functions. 


♦0*(O  =  det(D(s)) 

(3.25) 

♦c4C)  =  det(w(s)  +  °(s)) 

(3.26) 

♦0*(s)  =  det(5(s)) 

(3.27) 

♦c*(s)  =  det(W(s)  +  D(s)) 


(3.28) 
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They  are  related  to  the  return  difference  by 


det(I  +  G( s ) }  = 


♦of(s) 


det(I  +  G(s)  = 


♦c£> 


(3.29) 


(3.30) 


By  the  coprimeness  of  N,  D  (N,  D) ,  s  Is  an  m-th  zero  of 

det(I  +  G(s))  (det(I  +  G(s)))  if  and  only  if  s  is  an  m-th  order  zero 

cf  $c^(s)  (4>ce(s)h  Also  s  is  an  m-th  order  pole  of  det(I  +  G(s)) 

/v 

(det (I  +  G(s)))  if  and  only  if  s  is  an  m-th  order  zero  of  $o^(s)). 
Therefore  the  closed-loop  system  is  stable  if  and  only  if  ^c^(s) 
(*c^(s))  has  no  zero  in  the  CRHP. 

We  can  now  state  the  robustness  theorem  for  additive  modeling 
errors  as  follows. 


Theorem  3.4  (Robustness  Theorem  for  Additive  Error): 

The  perturbed  system  of  Figure  3.3  is  closed-loop  stable  (in 

A/ 

the  sense  that  $c^(s)  has  no  zeros  in  the  CRHP)  if  the  following 
conditions  hold: 

(1)  a.  $0g(s)  and  ♦  g(s)  have  the  same  number  of  zeros  in  the  CRHP; 
b.  $C£(S)  has  no  zeros  in  the  CRHP. 

(2)  CTmin^  +  C(JW)3  >  amax[E(Ju)]  for  a11  w  £  R+- 


m 


nuu  in 
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In  the  above  $ofi(s),  4cg( s).  ^(s).  <f>c£(s)  are  defined  by 
Eqs.  (3.45)-(3.28) .  The  additive  error  matrix  E(s)  is  defined  by 
Eq.  (3.24).  The  notations  a  .  [•]  and  a  [•]  denote  the  minimum 
and  maximum  singular  values  of  the  matrix. 

Proof  of  Theorem  3.4 

This  proof  is  very  similar  to  that  of  Theorem  2.4.  Let  X  e. 
[0,1].  It  is  well  known  (e.g.  [5]  [44])  in  singular  value  theory 
that 


amin^  +  +  >  ^in^1  +  "  ^max^0^ 


(3.31) 


Therefore  condition  (2)  guarantees  that  I  +  G(ju)  +  XE(joj)  is 
nonsingular  for  all  u  e  R+. 

Now  suppose  that  as  X  is  varied  continuously  from  zero  to 
unity,  the  number  of  encirclements  of 


f(X,u)  =  det  (I  +  G(ju)  +  XE(ju)) 


(3.32) 


around  the  origin  changes.  Since  f(X,u)  is  continuous  in  (X.u)  e 
[0,  1]  x  R+,  its  locus  on  the  R+  forms  a  closed  bounded  contour  in 
the  complex  plane  for  any  X  e  [0,  1].  Therefore  the  only  way  for  a 
change  in  the  number  of  encirclements  to  occur  is  for  the  locus  of 
f(Xg,<j)  to  pass  through  the  origin  for  some  Xq  e  [0,  1].  This  is 
equivalent  to  requiring 
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det(I  +  G(uq)  +  XqE(Wq))  =  0  for  some  (Xq.Uq)  e  [0.  1]  x  R+ 

(3.33) 

However  we  have  shown  that  I  +  G(j(i>)  +  X£(ju)  cannot  be  singular  by 

iV 

assuming  condition  (2).  Hence  we  conclude  that  det(I  +  G(jw)), 
corresponding  to  X  =  1,  encircles  the  origin  the  same  number  of 
times  as  det(I  +  G(j<i>)),  corresponding  to  X  =  0,  along  the  path 
u  e  R+.  By  applying  Theorem  3.3  (Nyquist  Theorem),  we  see  that  the 
perturbed  system  is  closed-loop  stable.  (End  of  Proof) 

Another  useful  way  to  describe  the  perturbed  system  with 
respect  to  the  nominal  system  is  to  define  a  multiplicative  factor 
matrix  L(s)  by 

G(s)  =  G(s)  L(s)  (3.34) 

This  representation  generalizes  the  gain  and  phase  margin 
description  of  the  single-input-single-output  case,  and  is 
particularly  useful  to  describe  the  robustness  properties  of  the  LQ 
system.  The  corresponding  robustness  theorem  for  multiplicative 
modeling  error  is  the  following. 
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Theorem  3.5  (Robustness  Theorem  for  Multiplicative  Error): 

The  perturbed  system  of  Figure  3.3  is  closed-loop  stable  (in 

/v 

the  sense  that  $cg(s)  (Eq.  (3.28))  has  no  zero  in  the  CRHP)  if  the 
following  conditions  hold: 

<v 

(1)  a.  $Qg(s)  and  4>0g(s)  have  the  same  number  of  zeros  in  the  CRHP; 
b.  <t>cg(s)  has  no  zeros  in  the  CRHP; 

where  and  are  defined  in  Eqs.  (3.25)  -  (3.27). 

(2)  L(jcj)  has  no  eigenvalue  in  R+1  for  all  w  e  R+. 

(3)  amin[I  +  G(Ju)]  >  CTmax  for  a11  u  fc  V  (335) 

The  proof  of  this  theorem  is  essentially  the  same  as  that  of 
Theorem  2.5,  and  is  therefore  omitted. 

3.3.4  Robustness  Properties  of  Distributed  LQ  Systems 

The  YFDT  presented  in  Section  3.2  for  infinite  dimensional  LQ 
systems  contains  the  KFDI  for  LQHD  systems  as  a  special  case  when 
the  control  space  is  finite  dimensional.  Since  by  using  the  algebra 

A 

S5.  we  were  able  to  characterize  the  robustness  of  distributed  linear 
time- invariant  systems  in  the  same  way  that  we  characterize  the 
robustness  of  LHDF  systems,  the  robustness  results  for  the  LQHD 
optimal  system  in  Chapter  2  can  be  extended  to  include  a  much  wider 
class  of  distributed  systems. 


For  the  purposes  of  investigating  the  robustness  properties  of 
the  LQ  system,  depicted  in  Figure  3.1,  it  is  convenient  to  represent 
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the  model  uncertainties  as  a  multiplicative  perturbation  factor,  as 

shown  in  Figure  3.4.  In  this  configuration,  the  matrix  L(s) 

represents  the  model  uncertainties  of  the  system  reflected  to  the 

point  X  in  Figure  3.1  where  the  loop  is  broken,  as  a  multiplicative 

perturbation  in  the  control  channel.  The  matrix  G(s)  is  the  nominal 

system  loop  transfer  function  matrix,  given  by  Eq.  (3.11). 

Nominally  L(s)  is  the  identity  matrix  and  the  LQ  system  is  closed 

loop  stable.  However  because  of  modeling  errors  and  implementation 

constraints,  the  matrix  L(s)  is  subject  to  changes  from  its  nominal 

value.  Consequently,  the  product  C(s)L(s)  can  be  viewed  as  the 

actual  system  loop  transfer  function  matrix.  This  closed-loop 

system  is  stable  if  and  only  if  the  closed-loop  system  of  Figure  3.5 

1/2 

•is  stable.  In  the  figure,  R  is  a  square  root  of  the  control 
weighting  matrix,  in  the  sense  that 


R  =  R1/2  R1/2 


(3.36) 


With  this  formulation,  we  can  state  some  of  the  robustness 
properties  of  the  LQ  system  as  the  following  theorem. 


Theorem  3.6  (Stability  Margin  of  Distributed  LQ  Systems! 

The  mul tipi icatively  perturbed  LQ  system  in  Figure  3.4  is 

/v 

closed-loop  stable,  i.e.  the  characteristic  function  $c^(s)  (Eq. 
(3.28))  has  no  CRHP  zeros,  provided  that  the  following  conditions 
are  satisfied: 
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I  Perturbed  system  G($)  | 


Figure  3.4  Mul tipi icatively  Perturbed  Feedback  System 


the  Derivation  of  LQ  Robustness  Properties 
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(1)  G(s)  is  specified  by  Eqs.  (3.3),  (3.4),  (3.11),  as  the 
closed-loop  solution  to  the  LQ  control  problem. 

<v 

(2)  ^0^(s)  and  have  the  same  number  of  CRHP  zeros, 

<v 

where  ^^(s),  4>0g(s),  are  defined  in  Eqs.  (3.25),  (3.27), 

a/ 

respectively,  and  both  $0g(s)  and  ^0^(s)  have  no  zeros  on  the 
imaginary  axis. 

(3)  L(ju)  has  no  eigenvalue  in  R+. 

(4)  ct  (R1/2L_1( ju)R_1/2  -  I)  <  1  for  all  u  e  R+. 

Proof  of  Theorem  3.6 

By  theorem  3.2,  $c^( s)  (Eq.  (3.26))  has  no  zero  in  the  CRHP. 
Note  that 

det[I  +  (R1/2G(s)R“1/2)(R1/2L(s)R'1/2)]  =  det[I  +  G(s)L(s)] 

(3.37) 


Hence 


det[I  +  (R1/2G(s)R~1/2)(R1/2L(s)R'1/2)] 


*Qg(S) 


(3.38) 


Similarly, 


det[I  +  R1/2G(s)R_1/2]  = 


(s) 

JJ) 


(3.39) 


94 


Theorem  (3.1)  YFDT  Implies  that 


a  [I  +  R1/2C(Jw)R"1/2]  >1  for  all  w  e  R+ 


(3.40) 


Therefore  in  conjunction  with  Theorem  3.5,  condition  (4)  guarantees 


that  $c^(s)  has  no  zeros  in  the  CRHP.  (End  of  Proof) 


The  following  are  corollaries  of  Theorem  3.6.  Their  proofs  are 


omitted  since  they  are  the  same  as  the  standard  proofs  (e.g.  [5])  in 


the  finite  dimensional  case. 


Corollary  3.7 


The  multiplicatively  perturbed  LQ  feedback  system  is  stable  if 


amax(L  Wl)  <  [cond(R)]  *^2  ,  for  all  u  e  R+  (3.41) 


where  cond(R)  denotes  the  condition  number  of  R 


cond(R)  = 


a  (R) 
max'  ' 


<Wn<R) 


(3.42) 


defined  as  the  ratio  of  its  maximum  singular  value  to  its  minimum 


singular  value. 
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Corollary  3.8 

If  R  >  0  is  diagonal,  then  simultaneously  in  each  feedback 
loop,  the  LQHD  system  has  the  following  guaranteed  stability 
margins, 

1.  [1/2,  »]  gain  margin 

o 

2.  +60  phase  margin 

Remark :  Tsitsiklis  and  Athans  [25]  have  shown  that  the  robustness 
properties  in  this  corollary  actually  hold  for  nonlinear  systems, 
under  certain  assumptions.  Their  proof  is  basically  an  extension  of 
the  one  found  in  Yakubovich  [22].  Our  formulation  in  this  thesis, 
however,  provides  a  more  general  framework  for  characterizing 
robustness. 

Corollary  3.9 

If  R  is  block  diagonal  of  the  form 


and  L(s)  consists  of  crossfeed  perturbation  of  the  form 


L(  s)  = 


I  X(s) 

P  I 


(3.44) 


nBnmnilBnDiRQn  JWteseiwiaiHrassnHraMaaiaju^reuwwww'.;  nvuma 
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Then  the  LQ  system  Is  stable  If 


. . ,<'mln*V,l/2 

°max'X<J"»  <  [5^T] 


,  for  all  w  t  R 


(3.45) 


Similarly,  If  L( s)  Is  of  the  form 


L(s)  = 


I  0 
X(s)  I 


(3.46) 


then  the  LQ  system  Is  stable  if 


<WX(J"))  <  [Hoi1"  •  ^  “I*  "  ‘  R+ 

max'  2 1 


(3.47) 


In  particular,  if  R  =  pi,  where  p  is  a  positive  scalar,  then  the  LQ 
system  is  stable  if 


a  (X(Ju))  <  1  ,  for  all  u  t  R 
max  t 


(3.48) 


for  both  types  of  crossfeed  perturbations  (Eqs.  (3.44)  and  (3.46)), 


3.4  Lyapunov  Control  Systems 

Lyapunov  control  systems  form  a  subclass  of  LQ  control  systems, 
in  which  the  control  gains  is  computed  via  the  solution  of  a 
Lyapunov  equation  instead  of  the  Riccati  equation.  A  nice 


! 

I 

! 

I 


Introduction  to  the  finite  dimensional  Lyapunov  control  system  is 
presented  in  Wong  [12].  In  this  section,  we  discuss  the  semigroup 
description  of  Lyapunov  control  systems,  and  their  robustness 
properties. 

Let  A  be  the  Infinitesimal  generator  of  a  strongly  continuous 
semigroup  of  operators  over  X,  and  let  B  be  a  linear  operator. 
Suppose  A  is  Hurwitz.  Then  the  Lyapunov  controller  for  the  dynamics 

=  Ax(t)  +  Bu(t);  x(0)  given  (3.49) 

is  given  in  feedback  form  by 

u(t)  =  -  H*x( t)  (3.50) 

where  H  is  the  solution  to  the  following  set  of  equations, 

KA  +  A*K  +  Q  =  0  (3.51) 


KB  -  HR  =  0 


(3.52) 


The  self-adjoint  operators  Q  and 


(3.2)  as  in  the  LQ  problem  (Eqs. 

.*  ,  nl/2* 

A  and  Q 


R  are  defined  in 
(3.8)  -  (3.10)). 


Eqs.  (3.1)  and 
Also  assume  that 


are  stablizable. 


Since  A  is  assumed  to  be  Hurwitz,  the  solution  K  to  Eq.  (3.51) 
Is  at  least  positive  semi-definite.  By  direct  manipulation  of  Eqs. 
(3.51)  -  (3.52),  we  get 

K(A  -  BH*)  +  (A  -  BH*)*K  =  -Q  -  2 HRH*  (3.53) 

Thus  A  -  BH**  is  Hurwitz.  In  addition,  Eq.  (3.51)  can  be  written  as 

KA  +  A*K  +  (Q  +  HRH**)  -  HR H*  =  0  (3.54) 

This  shows  that  the  Lyapunov  control  system  belongs  to  the  class  of 
LQ  systems. 

Yakubovich  [23]  has  a  frequency  domain  theorem  which  is  useful 
for  studying  the  robustness  of  Lyapunov  control  systems. 

Theorem  3.10= 

For  the  existence  of  K  and  H  satisfying  Eq.  (3.51)  -  (3.52),  it 
is  necessary  and  sufficient  that 

RG(ju)  +  C*(ju)R  >  0  for  all  u  >  0  (3.55) 

where  G(ju)  is  the  loop  function  defined  by 
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Recall  that  a  perturbed  system  can  be  represented 
multiplicatively  as 

G( jw)  =  G( ju)L( Jw)  (3.57) 

The  Lyapunov  control  system  will  remain  stable  as  long  as  the 
multiplicative  perturbation  L(jtj)  is  "passive"  [40]. 

For  Theorem  3.11  and  its  Corollary  3.12,  we  assume  that  the 
control  space  U  is  finite  dimensional. 

Theorem  3.11: 

The  perturbed  system  will  remain  closed-loop  stable  as  long  as 
RL( j<*>)  +  L*(ju)R  >  0  for  all  u  >  0  (3.58) 

Proof  of  Theorem  3.11: 

This  proof  uses  an  imbedding  argument  similar  to  that  for 
Theorem  3.4.  In  this  case,  we  consider  the  continuous  function 

f(X.u)  =  det(I  +  G( J<d)(XL(j<i>)  +  (1  -  X)I))  (3.59) 

where  X  a  [0,  1],  and  u  >  0.  We  want  to  show  that  f(X,u>)  ^  0  for 
all  X  and  u.  Hence  there  is  no  change  in  encirclements  when  the 
graph  f(0,u)  is  continuously  collapsed  to  that  of  f(l,u). 


We  prove  this  by  contradiction.  Let 

P(X.u)  *  XL(Jw)  +  (1  -  X)I  (3.60) 

and  suppose  that  ^Xq.Oq)  =  0  for  some  Xq,  Wq.  Then  there  exists  a 
nonzero  vector  u  e  U  such  that 

[I  +  C(JWo)P(X0.u0)]u  =  0  (3.61) 

l.e. 

us-  G(Jw0)P(X0.<j0)u  (3.62) 

Define 

z  =  P(XQ.w0) u  (3.63) 

Then  Eq.  (3.62)  becomes 

z  =  -  P(X0,w0)G(jw0)z  (3.64) 


Thus 

z*RG(J<Jq)z  +  z*G*(j<i>0)Rz 

=  -  zV(jw0)[RP(X0.u0)  +  P**(X0.Uo)R]G(jW())z  >  0 


(3.65) 
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I 

I 

I 

V 

I 

Since  G(JWq)z  s  u  is  nonzero, 

RP(X0,<j0)  +  P**(X0,(J0) R  <  0  (3.67)  I 

but  this  contradicts  Eq.  (3.5S)  because  of  Eq.  (3.60). 

\ 

(End  of  Proof) 

I 

1*  ' 
“  2 

As  in  the  LQ  case  (see  Theorem  3.6),  by  considering  R 

1  1*  _  I  ! 

”  2  ~  2  9 

C(s)R  and  R  L(s)R  in  place  of  G(s)  and  L(s),  one  can  easily 

obtain  the  following  Corollary.  j 

Corollary  3. 12 

If  R  is  diagonal,  then  the  Lyapunov  control  system  has. 

simultaneously  in  all  control  channels,  at  least 
o 

(1)  +90  phase  margins, 

(2)  100%  reduction  gain  margins,  and 

(3)  00  gain  margins. 

In  Chapter  4,  we  study  a  Lyapunov  control  design  for  the 
vibration  suppression  of  a  flexible  beam. 

In  case  U  is  infinite  dimensional,  we  can  state  the  following 
theorem.  Part  (1)  of  this  theorem  is  an  extension  of  the  necessary 
part  of  Theorem  3.10.  Part  (2)  provides  a  test  for  stability, 
similar  to  the  condition  in  Eq.  (3.58),  except  here  one  must  go 
through  the  test  for  the  whole  CRHP. 
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Theorem  3.13: 

For  Lyapunov  control  systems, 

(1)  RC(s)  +  G*(s)R  >  0  for  all  s  t  CRHP 

(2)  The  perturbed  system  will  remain  stable  as  long  as 

RL(s)  +  L(s)R  >  0  for  all  s  t  CRHP 

Outline  of  Proof: 

The  result  (1)  can  be  easily  obtained  by  manipulating  Eqs. 
(3.51)  -  (3.52).  For  (2),  we  want  to  show  that  I  +  G(s)L(s)  is 
invertible  for  all  s  e  CRHP.  This  can  be  done  by  using  a 
contradiction  argument  similar  to  that  in  Eqs.  (3.60)  -  (3.67),  i.e. 
by  supposing  I  +  G(sQ)L(s0)  is  singular  for  some  sQ  e  CRHP. 

(End  of  Proof) 


CHAPTER  4 


SYNTHESIS  OF  A  ROBUST  CONTROLLER  FOR  1 

i 


THE  SUPPRESSION  OF  VIBRATION  OF  A  FLEXIBLE  BEAM 


4.1  Introduction 

The  robustness  characterization  that  we  discussed  in  Chapters  2 
and  3  are  useful  for  the  robustness  analysis  of  distributed  feedback 
systems,  regardless  of  what  the  control  design  method  is.  However, 
the  guaranteed  stability  margins  for  LQ  optimal  systems  that  we 
derived  has  one  additional  Important  impact.  They  allow  the  control 
system  designer  to  synthesize  closed-loop  stable,  realizable  (finite 
dimensional)  feedback  controllers  for  infinite  dimensional  systems 
without  going  through  any  trial  and  error  iteration  in  design.  As 
discussed  in  Chapter  1.  the  trial  and  error  iterative  process  is  a 
notable  shortcoming  of  design  methods  which  depend  on  the  use  of 
reduced-order  models.  In  this  chapter,  we  use  a  flexible  beam, 
simply  supported  at  both  ends,  to  illustrate  a  design  synthesis 
method  which  leads  to  a  stable  closed-loop  system  when  the  inherent 
damping  is  sufficiently  large. 

Also  we  use  this  opportunity  to  study  the  importance  of 
inherent  damping  in  large  flexible  structure.  Stein  and  Greene  [1] 
showed,  by  example,  that  sufficient  damping  must  be  present  in  order 
to  design  stable  finite  dimensional  compensators.  Our  results 
support  their  findings.  We  show  that  if  either  all  the  temporal 
modes  of  the  flexible  structure  have  the  same  damping  ratio  or  the 


same  damping  coefficient,  then  a  stable  finite-dimensional 
compensator  can  be  Implemented  by  using  spatial  sampling.  In  the 
latter  case,  a  dynamical  compensator  is  required  to  attenuate  the 
high  frequency  components  in  addition  to  the  sampled  measurements. 

We  organize  this  chapter  as  follows.  In  Section  4.2  we 
describe  a  uniform  slender  beam,  simply  supported  at  both  ends.  In 
Section  4.3  we  describe  the  system  we  want  to  control.  A  point 
actuator  is  used  to  control  the  vibration  of  the  beam  and  the 
damping  ratio  is  assumed  to  be  a  constant  for  all  modes.  Also  we 
present  an  open-loop  analysis  of  the  beam.  In  Section  4.4  we 
formulate  a  Lyapunov  control  problem  for  the  flexible  beam  using  the 
semigroup  framework.  The  closed- loop  behavior  of  the  resulting 
optimal  feedback  system  is  discussed.  In  Section  4.5  we  show  that 
the  continuous  distributed  measurements  required  by  the  Lyapunov 
controller  can  be  adequately  approximated  with  finitely  many 
samples.  In  Section  4.6  we  take  a  closer  look  at  the  relationship 
between  the  size  of  the  structure’s  inherent  damping  and  the  density 
of  spatial  sampling  required  to  maintain  a  stable  approximation. 

In  Sections  4.3-4. 6,  a  constant  damping  ratio  is  assumed  for 
all  modal  deflections.  In  Section  4.7,  we  consider  instead  the  case 
in  which  the  inherent  damping  corresponds  to  a  constant  damping 
coefficient  for  all  modes.  In  this  case,  since  the  open-loop  system 
frequency  response  has  an  envelope  of  peaks  that  does  not  roll  off. 


a  high  frequency  attenuator  is  used,  in  addition  to  Increasing  the 
number  of  spatial  samples. 


4.2  A  Simply  Supported  Flexible  Beam 

The  dynamics  of  a  flexible  beam  freely  hinged  at  both  ends  with 
constant  density  per  unit  length,  denoted  by  p,  and  constant 
flexural  rigidity  El,  can  be  represented  by  the  following  system  of 
partial  differential  equations  [47] 

a2  a  a4 

P  — 9  y(x.t)  +  L  T-  y(x,  t)  +  EI  —T  y(x,t)  =  w(x)u(x,t) 

atz  01  ax 

(4.1) 


along  with  the  boundary  conditions 

2  2 

y(0,  t)  ,  y(a.t)  =0;  ^-s-y(O.t)  ,  ^-s-y(a.t)  =0  .  (4.2) 

dx  dx 

The  independent  variable  x  is  the  longitudinal  coordinate  ranging 
from  0  to  a,  where  a  is  the  length  of  the  beam,  and  the  variable 
y ( • . •) .  which  depends  on  x  and  time  t,  denotes  the  vertical 
displacement  of  the  beam  from  its  equilibrium  position.  The 
function  u(x,t)  and  the  product  w(x)u(x,t)  represent  the  control 
input  and  the  force,  respectively,  acting  on  the  beam  at  x  and  t. 
Figure  4.1  illustrates  these  arrangements.  The  control  input 
function  u(*,»)  is  an  element  of  ^([O.a^tO,00) ;  R).  The  function 
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w(x)  Is  real  and  may  vanish  for  some  x  and  is  integrable  i.e. 

d 

w(*)  fc  Lj.  The  damping  term  L  ^-y(x.t)  represents  the  inherent 
energy  dissipation  [47]  [48]. 

The  assumption  of  an  Inherent  damping  term  is  extremely  crucial 
to  vibration  control  problems,  as  evident  from  the  following  in 
literature. 

(1)  Gibson  [28]  showed  that  positive  definiteness  of  L  is  a 
sufficient  condition  for  his  approximate  method,  based  on 
optimal  control  of  a  sequence  of  reduced-order  modal  models,  to 
converge,  and  for  the  limiting  control  law  to  be  closed-loop 
stable. 

(2)  Gibson  [27]  showed  that  open-loop  uniformly  exponential 
stability  is  a  necessary  condition  for  the  existence  of  a 
compact  feedback  compensator  that  is  closed-loop  stable. 

(3)  Stein  and  Greene  [1]  used  classical  control  arguments  to 
demonstrate  that  a  practical  implementation  of  a  closed-loop 
stable  finite  dimensional  feedback  compensator  is  impossible 
unless  sufficiently  large  inherent  damping  is  assumed. 

In  addition,  we  have  discussed  in  Chapter  3  that  it  is  possible  to 
use  the  infinite  dimensional  robustness  theorems  when  a  sufficient 
amount  of  inherent  damping  is  present. 


WY1 
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4.3  Open-Loop  Analysis 

Beginning  in  this  section  we  study  a  simplified  model  of  the 
flexible  beam.  We  consider  the  single  actuator  case  in  which  a 
point  actuator  is  assumed  to  be  situated  at  the  midpoint  of  the 
beam,  exerting  a  force  on  the  beam  for  all  time  t  >  0.  Also  we 
assume  that  the  operator  L  which  accounts  for  the  inherent  damping 
in  the  beam  dynamics  is  given  by 


L  = 


.  C  >  0 


(4.3) 


where  f  is  a  small  positive  number  representing  a  small  constant 
damping  ratio  for  the  beam’s  modal  responses.  In  a  moment,  we  will 
clarify  this  in  terms  of  a  series  of  second-order  transfer 
functions. 

In  addition,  referring  to  Eqs.  (4.1)  -  (4.2).  we  let,  without 
loss  of  generality, 


a  =  1 


P  =  1 


El  =  1 


(4.4) 

(4.5) 

(4.6) 


These  assumptions  simplify  the  ensuing  analysis  considerably. 


Rewriting  Eqs.  (4.1)  and  (4.2)  under  the  assumptions  of 
Eqs.  (4.3)  -  (4.6),  we  obtain  the  following  set  of  equations  for  the 
dynamics  of  the  uniform  flexible  beam  simply  supported  at  both  ends. 


4  2  2 

^-4  y(x. t)  ■  c  ^2  ft  y(x- 0  +  ^“2  y(x- 0  =  «(t)*s(x  -  A)  (4.7) 

dx  dx  dt 

where  6(x  -  ^ )  is  the  spatial  Dirac  delta  function  at  x  =  The 
boundary  conditions  that  govern  the  dynamics  are 

y(O.t)  =  y(l.t)  =  0  (4.8) 

2  2 

y(O.t)  *  2— y(l.t)  =  0  (4.9) 

ax'  ex¬ 

it  is  easy  to  show  (Appendix  C)  that  the  solution  to  Eqs. 

(4.7)  -  (4.9)  is  representable  in  modal  form, 

00 

y(x,  t)  =  ^  sin  (ifrxJT^t)  (1.10) 

i=l 

where  the  functions  sin(iirx)  and  T^(t)  are  the  so-called  mode  shapes 
and  modal  deflection  of  the  beam.  The  modal  deflection  T^t)  is 
given  by  the  equation 


I 


T,(t)  +  CiV  |j-Tt(t)  +  i4ir\(t) 


2u(t)sin  (3j)  ;  (4.11) 


i  *  0,  1 , . . . 


Observe  chat  because  of  the  boundary  conditions  in  Eqs.  (4.8)  - 

(4.9),  the  mode  shape  functions  v£  sin  (Ittx),  i  =  1,2,3 .  form  a 

complete  orthonormal  basis  for  the  beam  shape  between  x=0  and  x=l . 

A  A 

Let  u(s),  T^s)  denote  the  Laplace  transforms  of  u(t)  and 
Ti(t),  respectively.  It  is  useful  to  write  down  the  transfer 

A 

function  from  the  control  of  u(s)  to  each  individual  modal 

A 

deflection  T^s).  This  is  obtained  in  a  straightforward  fashion 
from  Eq.  (4.11). 


T^s)  2sin(^|) 

~  .  =  2  r .2  2  .44 

u(s)  s+j;iTrs+i7r 


i  =  1,  2.... 


(4.12) 


By  examining  this  modal  deflection  transfer  function,  we  see 

2  2 

that  the  temporal  frequency  is  i  w  ,  and  the  damping  ratio  is  f. 

This  explains  our  earlier  claim  in  Eq.  (4.3)  that  f  corresponds  to  a 
constant  damping  ratio  for  all  deflection  modes. 

Let  Tj(t)  denote  the  rate  of  change  of  T^t),  i.e. 


ti(o  =  ^Tl(t) 


(4.13) 


Ill 


Then 


Tt(s) 

u(s) 


2s  sin  (^|) 


2  A  -.2  2  ,  .44 
s  +fiirs+iir 


1  =  1 ,  2, . . . 


(4.14) 


where  T(s)  is  the  Laplace  transform  of  T^(t).  This  last  transfer 
function  is  important  here  because: 

(1)  we  assume  that  the  available  sensors  are  of  the  linear-velocity 
type; 

O 

(2)  it  introduces  a  90  phase  lead  which  simplifies  the  control 
design  task. 

Let  y(x,t)  denote  the  time-derivative  of  y(x,t), 

y(x.t)  *  ft  y(x<t>  (4.15) 


Then 


y(x.t)  =  )  sin(iirx)T,  (t)  (4.16) 

i=l  1 

Observe  that  in  Eqs.(4.12)  and  (4.14),  since 

sin  (^)  =  0  for  i  even  (4.17) 


the  even  modes  of  the  beam  are  not  controllable  by  u(t).  However 
all  modes  are  stabillzable  because  of  the  inherent  damping 


mechanism.  Any  vibration  suppression  control  problem  for  the  system 
of  Eqs.  (4.7)  -  (4.9)  is  in  effect  a  control  problem  for  the  odd 
numbered  modes  only.  This  is  a  consequence  of  our  choice  of  placing 
the  actuator  at  the  mid-point  of  the  beam.  By  placing  the  actuator 
elsewhere,  for  example  at  an  irrational  point,  complete 
controllability,  albeit  small  for  some  modes,  can  be  obtained. 

4.4  Optimal  Lyapunov  Control 

In  this  section  we  derive  the  optimal  Lyapunov  control  for  the 
system  in  Eqs.  (4.7)  -  (4.9).  A  semigroup  formulation  is  used. 

Recall  from  Chapter  3  that  the  Lyapunov  control  system  is  a 
subclass  of  the  LQ  control  system.  The  Lyapunov  control  system  has 
the  additional  properties  that  it  is  in  general  easier  to  compute 
and  has  superior  stability  margins. 

Schaechter  [49]  considered  the  LQ  formulation  for  the  control 

of  the  flexible  beam,  and  also  for  the  estimation  of  distributed 

parameter  systems  in  [50].  In  both  cases,  the  optimal  solution 

calls  for  solving  a  matrix  Riccati  equation  over  a  ring  of 

polynomials  [72],  i.e.,  the  matrices  in  the  equation  as  well  as  the 

solution  contain  elements  which  are  polynomials  in  D,  where  D 

d 

denotes  the  operator  The  solution  of  this  Riccati  equation, 
however,  is  not  easy  to  obtain  [51].  In  addition,  the  exact 
solution  is  generally  Irrational,  i.e.  it  is  not  an  element  of  the 
quotient  field  [72]  of  the  ring  of  polynomials,  and  hence  is  often 
difficult  to  characterize.  Brockett  and  Willems  [52]  considered  a 
wider  class  of  LQ  problems  for  partial  differential  equations,  in 


their  generalized  path  integral  method.  Their  method  also  required 
the  control  system  designer  to  deal  with  irrational  functions  of  D. 

We  use  instead  an  approach  by  solving  a  Lyapunov  equation.  The 
linearity  of  the  Lyapunov  equation  eliminates  the  aforementioned 
pitfalls  associated  with  solving  the  Riccati  equation.  Exact 
solutions  [73]  can  be  obtained  for  this  Lyapunov  formulation.  In 
addition,  the  Lyapunov  control  sets  up  a  framework  in  which  one  can 
study  the  dependence  of  stability  on  the  Inherent  damping. 

We  rewrite  Eq.  (4.12)  in  the  following  semigroup 
representation. 

0 

l  u(t) 

6(x  -  \) 

(4.18) 


The  spectrum  of  the  infinitesimal  generator  of  the  semigroup. 


(4.19) 


lies  entirely  in  the  open  left-half  plane,  thus  permitting  the  use 
of  the  Lyapunov  method. 
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For  the  Lyapunov  design,  we  choose 


Q  = 


0 

0 


0 

2f. 


(4.20) 


and 


R  =  1  (4.21) 

This  choice  of  Q  will  lead  to  a  control  gain  that  Is  Independent  of 
the  damping  ratio.  The  operator  A  maps  the  space  S(0,1)  x  S(0,1) 
into  itself,  where  S(0,1)  is  the  set  of  smooth  functions 
f:  [0.1]  -»  R 

S(0,1)  =  { f ( •)  |  H  exists  for  all  n;  f(0)  =  f(l)  =  0 
dx11 

2  2 

and  11  (0)  =  ^4  (U  -  °)  H-22) 

dx*  dx* 


It  is  trivial  to  verify  that  S(0,1)  is  a  separable  Hilbert  space  to 
which  we  can  associate  the  inner  product  <•,•> 


and  with  an  abuse  of  notation,  the  Inner  product  for  the  space 
S(0,1)  x  S(0,1)  Is 

fl  S1  rl 

<  .  >  =  J  C^OOgjM  +  f2(s)g2(x)]dx; 

f 2j  g2j  0 

frf2,gl,g2  £  S(0,1)  (4.24) 

Before  turning  to  the  derviation  of  the  optimal  controller,  we 
would  like  to  show  that  the  flexible  begun  control  system  belongs  to 
the  class  of  systems  discussed  in  Section  3.2;  hence  we  can  use  the 
robustness  results  in  Chapter  3  for  this  application. 

Referring  to  Eqs.  (3.15),  for  the  flexible  beam, 


Clearly  for  all  mode  shapes  sin(i7rx), 


sin(iirx) ,  sin(iirx)>  >  sin^(iTrx) 
dx 
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and  for  small  f. 


2  4  4 

<~C  ^-s- sin(lirx).  sin(i7rx)>  <  sin(iTrx),  sin(iirx)> 

dx  "  dx  dx 


Therefore  Eqs.  (3.16a)  and  (3.16b)  are  satisifed,  and  it  Is 
legitimate  for  us  to  use  the  LQ  robustness  properties  to  design  the 
vibration  suppression  controller. 

Let  A  be  the  adjoint  operator  of  A  with  respect  to  the  inner 
product  defined  by  Eq.  (4.24).  i.e.  for  all  fj.fg.gj.gg  £■  S(0.1). 


V 

< 

.  A 

f2. 

g2 

>  =  <  A 


M 

• 

[f2j 

M 

Then  by  direct  comparison  of  the  two  sides  in  Eq.  (4.25), 


dx 

a2 

r  - — 

»  2 
dx 


We  form  the  following  Lyapunov  equation. 


A  K  +  KA  =  -  Q 


(4.25) 


(4.26) 


(4.27) 


& 
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where  K:  S(0,1)  x  S(0,1)  -»  S(0,1)  x  S(0,1)  is  a  self-adjoint  and 
positive  definite  linear  operator.  The  solution  K  to  this  equation 
is  unique  [23]  and  is  given  by  Eqs.  (4.28)  -  (4.30)  below. 


K  = 


K, 


0 


(4.28) 


K 


1 


(4.29) 


K2^~2  =  -  1  (4-30) 

3x 


where  and  K g  are  self-adjoint  and  positive  definite  linear 
operators  on  S(0,1).  Note  that  conceptually.  Kg  can  be  intepreted 
as  the  double  integral  operator,  which  satisfies  the  appropriate 
boundary  values. 


MHiNi  PH  JWJUH  mi  mm 
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The  Lyapunov  control  law  constructed  by  using  the  weights  Q  and 
R  specified  by  Eqs.  (4.20)  and  (4.21)  is  given  in  the  feedback  form 
as 


u(t)5(x  -  ~) 


y(x. t) 

ft  (Xlt) 


(4.31) 


where  G:  S(0,1)  x  S(0,1)  «*  is  the  map 


G  =  R-1  [0  6(x  -  |)]K 


i  .e. 


G  =  [0  fi(x  -  |)K2] 


(4.32) 


Or  simply 


The  operator  Kg  can  be  represented  in  terms  of  a  Green's  function 
[69]  [70]  h(x.f), 


Kg f  =  J1  h(x,f )f (f )df  :  f  e  S(0.1) 


(4.34) 


where 


h(x,f)  =0  at  f  =  0  and  f  =  1 


(4.35) 


2-5-  h(x.f)  =  0  at  f  =  0  and  f  =  1 

ar 


(4.36) 


The  Green's  function  is  determined  in  Appendix  D.  It  is  given  by 


h(x.f)  =  2  ~22  sin(iirx)  sin(iirf) 
i=l  i  v 


(4.37) 


Combining  this  result  with  Eqs.  (4.31)  -  (4.34)  we  can  write  down 


the  optimal  Lyapunov  control  law  as 


u(t)6(x  -  i)  =  [  [  -sto  sln(i™)  sin(iirf)] 

J0  i=l  in 


*  6<x  '  b  Ik  Cf  -  c)df 


(4.38) 
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Integrating  both  sides  with  respect  to  x  to  get  rid  of  the  delta 
functions,  we  obtain 

u(t)  =  f  C  5  sin(5I)  sln(iirf)]  (f.t)df  (4.39) 

J0  i=l  iV  z  ot 

To  derive  the  loop  transfer  function  for  this  optimal  system, 
breaking  the  loop  at  the  control,  we  note  from  Eq.  (4.16)  that 

§£(f.t)s]|  sin(krrf)  T.(t)  (4.40) 

k=l 


Substituting  this  into  Eq.  (4.39)  we  get 


[  )  -5%  sin(|£)  sin(iirf)]  [  )  sin(kmf )T.  (t)]df 
i=l  i  ir  *  k=l  k 

(4.41) 


Because  of  orthogonality,  this  is  equivalent  to 

u(t)  =  \  ~Y2  sin(|ZL)  (4.42) 

k=l  k  t r 

Combining  this  with  Eq.  (4.14)  we  obtain  the  loop  transfer  function 
for  this  optimal  feedback  system. 


g(s)  = 


i  2  2 

k=l  k  ir 


.  2,k7T. 
sin  (s-) 


2s 


r,2  2  , 

fk  ir  s+ 


.4  4 
k  ir 


(4.43) 


s-r.  r,  iv  -f *.< 4 


This  loop  transfer  function  has  several  properties  which  are 
significant  for  the  control  design. 

According  to  Chapter  3.  g(s)  is  passive  [40]  [53],  i.e. 

g(jw)  +  g*(ju)  >  0  for  all  w  t  R+  (4.44) 

This  can  be  readily  verified  by  inspecting  Eq.  (4.43),  since 
2  lCTT 

sin  {—)  is  either  0  or  1 .  The  passivity  guarantees  that  the 
optimal  system  is  closed-loop  stable  (Willems  [53]).  In  fact  one 
can  characterize  the  additive  stability  margin  of  this  Lyapunov 
control  system  as  follows. 

Lemma  4.1  (Guaranteed  Additive  Stability  Margin  for  Lvapunov 
Control! • 

Let  g(s)  be  the  loop  transfer  function  of  a  unity  feedback 
system,  and  g(s)  has  no  open-loop  poles  in  the  CRHP.  Then  this 
system  is  closed- loop  stable  if 


|g(Ju)  ~  g(ju)l  <  (1  +  |g(ju)|2)1/2  for  all  w  e  R+  (4.45) 


where  g(ju)  is  the  loop  transfer  function  of  the  optimal  Lyapunov 
control  system. 


Proof  of  Lemma  4.1 
|| 

Adding  1  +  g(Ju)g  (Ju)  to  both  sides  of  Eq.  (4.44),  we  get 


(1  +  g(J«))(l  +  g*(ju))  >  1  +  g(ju)g*(J«) 


for  all  w  e  R+ 
(4.46) 


or  simply, 


|l  +  g(jw)|2  >  1  +  |g(Jw) |2  for  all  u  t  R 


(4.47) 


Thus  if  the  condition  in  Eq.  (4.45)  holds,  then  by  the  triangle 
Inequality,  i.e. 

|l  +  g|  =  |l  +  g  +  g  -  g|  >.  |l  +  g|  -  |g  -  s| 


we  get 

|l  +  g(ju>)|  >  0  for  all  w  e  R+  (4.48) 

Hence  the  system  is  closed  loop  stable  by  the  Nyquist  theorem 
(Willems  [53]).  (End  of  Proof) 

Now  let  us  examine  the  shape  of  g(s).  We  show  in  Appendix  E 
that  at  low  frequencies,  |g(ju)|  consists  of  a  series  of  peaks  at 
the  natural  frequencies 


K 


% 

mmM 


1 

' 

I 

In  addition,  |g(Jw)|  has  a  finite  cutoff  frequency  u  ,  i.e.  |g(jw)| 

c 

is  less  than  unity  for  all  u  >  u  .  At  high  frequencies,  |g(ju)| 

c 

rolls  off  at  the  rate  of  20  db/decade.  This  behavior  is  shown  in 
the  frequency  plot  in  Figure  4.2. 

i 

I 

4.5  Implementation  of  a  Closed-Loop  Stable  Suboptlmal  Design 

The  optimal  feedback  control  law  that  we  derived  in  Section  4.4 
cannot  be  realized  exactly  in  actual  practice.  It  requires 
continuous  measurements  along  the  entire  length  of  the  beam  as 
indicated  by  Eq.  (4.33)  or  complete  knowledge  of  all  modes  as 
Indicated  by  Eq.  (4.42).  Neither  of  these  is  achievable  physically. 

We  consider  using  a  suboptlmal  control  of  the  form 

M 

u(t)  =  1  ~T2 sin  (r}  ^k(t)  <4-51> 

k=l  k  ir 

where  T^( t )  is  an  estimate  of  T^(t)  from  the  measurements  (x.t) 
at  N  discrete  locations,  described  below,  along  the  beam.  The 


estimate  T^(t)  is  obtained  by  correlating  the  measurements  with  the 


log  freq 


6 

E+00 


Figure  4.2  Optimal  Loop  Function 
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k-th  mode  shape.  In  practice,  the  series  summation  over  k  cannot  be 
computed  for  infinitely  many  terms.  Therefore,  we  assume  that  the 
series  summation  has  been  truncated  to  M  terms  in  actual 
implementation. 

We  assume  that  the  beam  is  divided  equally  into  N  parts  as 
follows. 


x 


=  0 


i  i 

N  N .  N 


A  point  sensor  that  measures  the  deflection  rate  is  situated  at  each 
of  the  marked  point.  In  other  words, 


z«(t)  =  (f-  0  :  *  =  0 . N 


(4.52) 


are  known  at  all  times. 

We  construct  the  correlation  estimate 

iftt)  I  sin  (f^)  z{(t) 

i=0 


Note  that  this  estimate  uses  an  approximation  of 


1  N_1 

J  sin(kirf)  §*■(£. t)df  :i  l  sin(^L)z^t) 

0  e=o 


(4.53) 


(4.54) 


In  the  sense  of  Riemann  integration.  A  comparison  of  this  with 


Eq.  (4.39)  motivates  our  attempt  of  using  the  suboptimal  control  in 


Eqs.  (4.51)  -  (4.53). 


We  summarize  the  suboptimal  control  scheme  in  Figure  4.3. 


First  we  compute  the  loop  transfer  function 


g(s)  = 

u(s) 


(4.55) 


We  show  in  Appendix  F  that 


g(s)  =  )  g,(s) 
i=l  1 


(4.56) 


where 


K^s) 


2s  sin  (j-) 


2  A  r,2  2  x  ,441 
s  +  fi  TT  s  +  i  IT 


(4.57) 


and  the  quantity  is  evaluated  in  Appendix  F.  For  M  sufficiently 


large,  it  is  given  for  i  odd  by 


,  ,  >  2  1  2  .iv, 

(-1)  — o  cosec  (ov)  +  O(jrj) 


(4.58) 


where  O(^)  is  a  number  of  the  order 
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Figure  4.3  Suboptimal  Control  Implementation. 
Constant  Damping  Ratio  Case 


Similarly  we  can  decompose  the  optimal  loop  transfer  function 


(Eq.  (4.43))  as 


g(s) 


(4.59) 


where 


g^s) 


„  .  2,iiTv 

2s  sin  (^-) 

2  _2  2  ^ 
s  +  fi  7T  s  + 


.4  4 
i  it 


.2  2 

i  ir 


(4.60) 


By  Lemma  4.1,  the  additive  stability  margin  of  the  optimal 
system  is  at  least  1.  Therefore  the  suboptimal  system  is 
closed-loop  stable  as  long  as 


|g(Ju)  -  g(ju) I  <  1 


for  u  e  R+ 


(4.61) 


In  Appendix  G,  we  show  that  there  exists  an  N  large  enough  such 
that  the  stability  condition  in  Eq.  (4.61)  is  satisfied.  In 
particular,  we  show: 

(1)  The  quantity  Ig^Jw)  -  g^jw)  |  is  arbitrarily  small  for  all  i 
and  for  all  w  >  0. 

2  2 

(2)  Let  (j  =  f  7 r  ,  and  k  is  the  nearest  odd  integer  to  f.  Then  the 


summation 
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l  |gt(J«)  -  gj(J«)|  <  0.05  (4.62) 

1=1 

i*k 

Remark:  The  bound  in  (2)  is  rather  conservative.  However,  it  is 

sufficient  for  our  design.  We  have  not  pursued  removing  this 

conservatism  in  this  thesis,  but  we  think  this  bound  can  be 

estimated  to  within  the  order  of  i. 

N 

Since 

|g(Ju)  -  g(jw)|  <  |gk(j<*>)  -  gk(Ju)| 

+  ^  IgjUu)  -  gj(J»)l  (4.63) 

1=1 

i/k 

Eq.  (4.61)  is  satisfied.  Consequently  the  suboptimal  system  is 
closed-loop  stable. 

4.6  Relationship  between  Sample  Density  and  Damping  Ratio 

Our  design  procedure  in  Section  4.5  calls  for  the  determination 
of  N  such  that  for  all  i. 

.  2  2  l]Z2  cosec  ^2N^  "  ~2~2  +  I  <  1  (4.64) 

fi  i  4N  i  7T 


v.y. 

-.  -r.K.  •  < 

£ L A:  ^  ..L  i*. 


We  discuss  in  Appendix  G  how  this  can  be  achieved.  A  finite  number 
of  constraints  are  used  to  determine  N  that  satisfies  Eq.  (4.64) 
which  guarantees  that  the  closed-loop  system  is  stable.  In 
particular,  we  require 

7"2  N  <  2  (4.65) 

where  0.95  >  6  >  0  is  chosen  by  the  control  designer  for  bounding 
the  error  term  0(-^).  This  constraint  Indicates  that  for  our  design 
procedure,  N  is  at  least  inversely  proportional  to  f.  Thus  when  f 
gets  smaller,  the  number  N  that  is  computed  by  the  procedure  gets 
larger.  Figure  4.4  summarizes  our  observation,  for  6  = 

We  remark  that  the  sample  measurements  are  meaningful  only  for 
N  >  2.  In  other  words,  even  if  the  damping  ratio  increases,  there 
must  still  be  at  least  one  sensor  at  the  mid-point  (the  other  is  at 
the  end-point  x  =  0  which  does  not  move)  of  the  beam  for  feedback 
measurements.  However,  when  this  stage  is  reached,  the  size  of  the 
damping  ratio  is  probably  so  large  that  it  would  have  deviated  from 
the  flexible  structure  assumption,  anyway. 

Remark:  numerical  evaluations  of  the  expression  inside  the  absolute 

value  signs  in  Eq.  (4.64)  for  some  i  and  N  indicate  that  the  error 

term  0 .  (rr)  dominates  over  the  difference  — ^  cosec^  (~) - 

4N  i  7T 

This  has  not  been  verified  rigorously,  however. 


where  T^(t)  is  related  to  u(t)  by  the  modal  transfer  function 


Tk(s)  2sin(jp) 

s  ~2~  A  .  4  4 
u(s)  S  +T7S  +  K7T 


(4.70) 


As  can  be  inferred  from  this  transfer  function,  the  natural 
2  2 

frequency  is  k  ?  for  the  k-th  deflection  mode,  and  the  damping 
coefficient  is  rj  for  all  modes.  Except  for  the  difference  in  the 
damping  coefficient  (in  the  denominator)  of  the  transfer  functions, 
this  open- loop  solution  is  the  same  as  the  previous  case,  where  the 


2  2 

damping  coefficients  are  fk  ir  . 

In  what  follows,  we  describe: 

(1)  how  one  can  derive  an  optimal  control  law  which  is  identical  in 
form  to  the  one  used  in  the  previous  case; 

(2)  an  attempt  to  design  a  suboptimal  controller  using  exactly  the 
same  strategy  as  described  in  Section  4.5; 

(3)  why  this  suboptimal  controller  cannot  be  guaranteed  closed-loop 
stability  no  matter  how  dense  one  samples  the  measurements, 
unless  a  dynamical  compensator  is  also  used; 

(4)  an  interpretation  of  the  above  results. 


4.7.1  Optimal  and  Suboptimal  Control  Laws 

As  in  the  constant  damping  ratio  case,  a  Lyapunov  control 
problem  can  be  formulated  and  solved.  The  details  of  these 
derivations,  however,  are  not  presented  because  they  are  almost 
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identical  to  those  in  Sections  4.4  -  4.5. 

Referring  back  to  Section  4.4,  the  infinitesimal  generator  of 
the  semigroup  (Eq.  4.19)  is  now 


A  = 


~V 


(4.71) 


and  we  use,  in  order  to  achieve  rolloff,  an  integral  operator  in  the 
control  weighting  operator 


Q  = 


0 

P 


0  ' 

2qK2. 


(4.72) 


where  the  operator  is  defined  in  Eqs.  (4.34)  and  (4.37). 

After  some  manipulations,  one  finds  that  the  optimal  Lyapunov 
feedback  control  is  given  in  the  same  manner  as  the  constant  damping 
ratio  case,  i.e.  Eqs.  (4.39),  (4.42), 

1  ^ 

u(t)  =  JQ  l  -jh  sinty  sin(i7rf )3  (f.t)df  (4.73) 


or  equivalently, 


u(t)  =  l  sin(ii)  T  (t) 

1=1  1  ’ 


(4.74) 
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The  optimal  loop  transfer  function  is 


g(s)  =  1  2  2  sin2^2)  2  2S  4  4 

. i  TT  S  +T)S+i  7T 


(4.75) 


Eqs.  (4.73)  -  (4.75)  motivate  us  to  try  to  construct  the 
suboptimal  control  exactly  as  before  (Eqs.  (4.51)  -  (4.54)), 


«<*)  *  l 

i=l  k  W 


(4.76) 


ft‘>  -  I  I  »‘"<TT>  £  <1  *> 


(4.77) 


The  loop  transfer  function  of  this  suboptimal  system  is  given  by 


r 

g(s)  =  2  g}(s) 


(4.78) 


~  2s  sin(-^J)  N 

gi(s)  =  “2 - 1"4  ci 

s  +77S+1  TT 


(4.79) 


„N  1  V  1 _  .  ,kir.  T  ,£kir x  ,lln\ 

“i  "  N  2  ,2  2  ir^2  >  2  ^  N  ^  sin(  N  * 


(4.80) 


This  coefficient  c^  is  identical  to  the  one  computed  for  the 
constant  damping  ratio  case,  in  Appendix  F.  Indeed  except  for  the 
change  of  the  damping  coefficient,  the  representation  of  g(s)  here 
is  the  same  as  the  one  formulated  for  the  constant  damping  ratio 
case. 
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N 

Using  the  computed  from  Appendix  F,  we  obtain 

2  r  1  2,1*.  1 _  ,  -,ln  . 

rj  ^  2  cosec  '22  °W'  1  odd 

4N  i  it 

0  ,  i  eve  (4.83) 

1  1  2  2 
where  O(^)  is  a  number  of  the  order  j^.  At  low  u  =  i  ir  ,  this  is 

approximately  the  peak  value  of  g(Ju)  -  g(ju),  at  low  frequencies. 

Let  us  assume  that  N  is  very  large.  Then  for  all  positive 

integers  m, 


~  ,  .  .2  2.  t  ..2  2. 
fiEi( Ji  »  )  "  g^Ji  ir  )  =  ■ 


cosec2  (fi™!!,  ~  ^ 


4N 


(4.84) 


O(^)  ~  0 


(4.85) 


(2mN+l)27r2 


(4.86) 


Thus  from  Eq.  (4.83), 


*2ntN+l(J(2mNt1)2’'2)  -  g9mN+!  (Jf2mN+l)2F2)  2  \  (4.87) 


2mN+l 


rjir 


In  general  q  is  an  extremely  small  number,  and 


>  1 


W 


(4.88) 
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2 

unless  tj  >  -g”.  or  approximately  tj  >  0.2.  Therefore  generally 

IT 

speaking  the  unity  guaranteed  stability  margin  of  the  LQ  system  is 
not  sufficient  to  establish  an  integer  N  for  sampling  such  that  the 
suboptimal  control  (Eq.  (4.76))  is  closed-loop  stable. 

The  above  conclusion  is  also  true  when  one  considers  the 
Lyapunov  guaranteed  stability  margin  instead  of  the  unity  LQ 
stability  margin,  even  though  the  Lyapunov  guaranteed  stability 
margin  (Eq.  (4.45)) 


(1  + 


2.4  4' 

TJ  1  IT 


1/2 


is  larger.  This  is  not  enough  because  when  i  is  sufficiently  large, 
the  term  -rp^-^will  no  longer  be  significant.  At  these  high 

T)  1  IT 

frequencies,  the  guaranteed  stability  margin  is  infinitesimally 
close  to  1. 

Figure  4.5  depicts  the  envelope  of  the  peaks  of  the  suboptimal 

loop  function  g(ju).  As  can  be  seen  from  the  figure,  the  loop 

2  2 

frequency  response  rolls  off  initially  until  the  frequency  N  ir  is 

reached,  and  since  the  sample  density  is  finite,  it  rises  back  to 

2  2 

the  original  peed:  value  at  (2N+1)  it  ,  which  is  shown  in  Eqs. 

(4.83)  -  (4.87).  This  phenomenon  repeats  indefinitely  as  frequency 
moves  higher. 

The  above  phenomenon  can  be  explained  by  looking  at  the 
open-loop  frequency  response  plot  of  any  actuator-sensor  pair.  For 
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Figure  4.5  Envelope  of  Peaks  for  Suboptimal  Loop 

Function.  Constant  Damping  Coefficient  Case 


example,  Figure  4.6  shows  the  loop  function  for  the  sensor  colocated 
with  the  actuator.  This  loop  function  has  a  constant  envelope  of 
the  peaks.  Thus  its  high  (spatial  and  temporal)  frequency 
components  are  aliased  to  lower  spatial  frequencies  by  the 
sampling-correlation  process,  and  are  then  multiplied  by  the  high 
gains  appropriate  for  the  lower  modes.  Consequently  we  have  a 
problem  with  the  stability  of  the  closed-loop  system. 

The  reason  why  this  does  not  create  a  problem  for  the  constant 
deunping  ratio  case  can  be  seen  by  observing  Figure  4.7.  It  shows 
that  for  the  constant  damping  ratio  case,  the  same  actuator-sensor 
pair  gives  rise  to  a  loop  function  whose  envelope  of  peaks  rolls  off 
indefinitely.  Therefore,  spatially  aliased  components  have  low 
amplitudes  and  do  not  cause  a  problem  in  that  case. 

The  above  difficulty  with  the  constant  damping  coefficient  case 
can  be  overcome  by  including  a  low  pass  temporal  filter  to  attenuate 
high  frequency  components.  This  corresponds  to  the  suboptlmal 
implementation  configured  in  Figure  4.8. 

Recall  that  a  Lyapunov  control  design  guarantees  that  the 
0 

system  has  at  least  90  phase  margins.  Thus  by  using  a  single-pole 

filter,  the  system  is  guaranteed  to  be  marginally  stable.  By 

placing  the  pole  of  the  filter  to  the  right  of  <j  ,  the  cutoff 

*  c 

frequency  of  the  inital  rolloff  of  |g(Jw)  |  in  Figure  4.5,  additional 
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stability  margin  is  obtained.  Also  one  should  place  sufficiently 
2  2 

to  the  right  of  N  rr  if  one  desires  to  retain  u  as  the  cutoff 

c 

frequency  of  the  overall  system.  Summarizing,  the  dynamical  filter 

^f  2  2 

has  a  transfer  function  of  — ■ — ,  where  u>  <  u.  <  N  ir  .  This 

S+OJj.  c  f 

subop timal  implementation  guarantees  that  the  closed-loop  system  is 
stable  when  N  is  large  enough.  Also  at  low  frequencies,  the  system 
behaves  approximately  as  the  optimal  system.  The  loop  function  of 
this  implementation  is  shown  in  Figure  4.5  for  a  comparison  with  the 
one  before  the  filter  is  inserted  in  the  loop. 


CHAPTER  5 


CONCLUSIONS  AND  SUGGESTIONS 
FOR  FUTURE  RESEARCH 

5. 1  Conclusions 

This  thesis  has  addressed  the  following  Issue.  Given  a 
distributed  linear-time-invariant  system,  there  are  two  sources  of 
errors  that  can  lead  to  a  difference  between  the  actual  implemented 
control  system  and  the  computed  nominal  control  system.  The  sources 
of  errors  are: 

(a)  modeling  errors  of  the  plant  dynamics;  and 

(b)  inexact  implementation  of  the  derived,  usually  infinite 
dimensional,  controller. 

It  is  imperative  that  the  control  system  designer  knows  how  much 
error  the  system  can  tolerate  without  becoming  unstable.  This 
knowledge  will  allow  him  to  synthesize  feedback  controllers  that  are 
implementable  and  closed-loop  stable. 

We  have  studied  this  problem  by  looking  at  the  robustness  of 
stability  for  distributed,  linear-time-invariant,  feedback  systems. 
We  have  been  successful  in  developing  the  following  results: 

(1)  A  robustness  characterization  for  linear  hereditary 

differential  systems  using  singular  values  (Ch.  2).  A  Nyquist 
theorem  for  this  class  of  systems  was  also  established. 


(2)  The  Kalman  frequency  domain  inequality  (and  equality)  for 
linear  quadratic  hereditary  differential  systems  (Ch.  2). 

(3)  A  set  of  guaranteed  stability  margins  for  the  optimal  linear 
quadratic  hereditary  differential  systems  (Ch.  2). 

(4)  A  robustness  characterization  for  distributed, 
linear-time-invariant  systems  (Ch.  3).  The  number  of  CRHP 
open-loop  poles  and  the  dimension  of  the  control  space  are  both 
assumed  to  be  finite. 

(5)  A  set  of  robustness  properties  for  linear  quadratic, 
distributed,  linear- time-invariant  systems  (Ch.  3). 

(6)  An  analysis  of  a  Lyapunov  controller  for  the  vibration 
suppression  problem  of  a  flexible  beam  (Ch.  4). 

(7)  An  analysis  of  the  role  played  by  the  inherent  damping  in  the 
synthesis  of  an  implementable,  closed-loop  stable  controller 
(Ch.  4).  Two  cases  are  studied  in  detail,  the  first  of  which 
corresponds  to  a  small  constant  damping  ratio  for  all  modes, 
and  the  second,  a  constant  damping  coefficient. 

(8)  Synthesis  using  sampled  measurements  along  the  beam  of  a 
suboptimal  closed-loop  stable  controller  when  either  a  constant 
inherent  damping  ratio  or  damping  coefficient  is  present 

(Ch.  4).  The  robustness  properties  of  the  optimal  system  are 
used  in  the  design  procedure. 

The  finite  dimensional  version  of  the  results  in  (1)  -  (5)  have 
been  known  for  some  time.  This  thesis  has  extended  the  finite 
dimensional  results,  while  at  the  same  time  maintaining  their 
practical  uses  in  an  infinite  dimensional  context. 


With  the  advent  of  techniques  for  computing  the  singular  value 
decomposition  of  a  matrix,  e.g.  Laub  [65],  Klema  and  Laub  [66], 
singular  values  have  proven  themselves  to  be  an  extremely  valuable 
tool  for  analyzing  multi-input-multi-output  feedback  control 
systems.  By  plotting  the  singular  values  of  the  return  difference 
over  all  frequencies,  the  class  of  modeling  errors  that  would  not 
destabilize  the  control  system  are  exposed. 

We  exploited  the  above  concepts  for  distributed,  linear¬ 
time-invariant  systems,  and  have  come  up  with  results  that  allow  the 
control  system  designer  to  use  the  singular  value  technique  almost 
in  the  same  way  as  in  the  case  of  finite  dimensional  systems.  The 
modification  here  is  that  the  plant  dynamics  are  no  longer 
represented  by  a  rational  transfer  function,  but  in  most  cases, 
their  representations  are  programmable  into  the  computer. 

Therefore,  there  is  practically  no  loss  in  utility  of  the  singular 
value  analysis  in  the  distributed  system  case. 

Invariably  literature  on  the  robustness  of  finite  dimensional 
systems  assumes  that  since  the  computed  controller  is  finite 
dimensional,  it  can  always  be  implemented  exactly.  However,  this  is 
not  true  for  distributed  systems. 

Generally  the  controller  design  approach  for  the  distributed 
system  fall  into  one  of  the  following  two  categories: 

(a)  Reduced-order  model  method,  and 

(P)  Infinite  dimensional  model  method,  followed  by  a  finite 

dimensional  approximation  of  the  resulting  infinite  dimensional 


controller. 


Referring  back  to  the  two  types  of  errors  introduced  at  the 
beginning  of  this  section,  the  approach  (a)  induces  the  (a)  type 


error,  i.e.  the  modeling  error  in  the  plant  dynamics,  and  the 
approach  (P)  induces  the  (b)  type  error,  i.e.  the  implementation 
error.  The  methods  developed  in  this  thesis  are  suitable  for  design 
analysis  for  either  of  these  approaches.  The  latter  approach  has 
been  examined  in  depth  in  this  thesis  using  a  vibration  suppression 
problem  for  an  idealized  model  of  a  flexible  beam.  Our  analysis 
reveals  that  the  inherent  damping  of  the  flexible  structure 
determines  the  degree  of  approximation  that  is  required  of  the 
suboptimal  Implementation  in  order  to  maintain  closed-loop 
stability. 

In  what  follows,  we  give  an  account  of  the  results  in  Chapters 
2-4,  and  then  a  few  suggestions  for  future  research. 

5.2  Account  of  Results  in  Chapters  2-4 

In  Chapter  2.  we  develop  a  robustness  characterization  of 
linear  hereditary  differential  systems  (LHDS),  for  both  additive  and 
multiplicative  perturbations.  By  showing  that  LHDS  have  finitely 
many  open-loop  poles  in  the  CRHP,  we  are  able  to  state  a  Nyquist 
theorem  for  LHDS.  which  provides  a  basis  for  robustness 
characterization.  As  in  the  finite  dimensional  case,  the  minimum 
singular  value  of  the  return  difference  is  used  to  characterize  the 
robustness  of  the  system  against  unstructured  error  in  the  system 
loop.  Because  of  this  formulation,  the  guaranteed  stability  margins 


for  the  linear  quadratic  optimal  system,  termed  the  linear  quadratic 
hereditary  differential  system,  are  derived  directly  from  the  Kalman 
frequency  domain  inequality.  We  have  shown,  by  directly  working 
with  complex-valued  matrix  functions  in  the  frequency  domain,  how 
one  can  derive  the  Kalman  Inequality  for  LQHD  systems. 

In  Chapter  3  we  generalized  the  results  in  Chapter  2  to 
distributed,  linear-time-invariant  systems  by  considering  a 
formulation  using  linear  operators.  This  is  achieved  for  the  case 
of  finite-dimensional  control  space  by  combinining  the  Yakubovich 
frequency  domain  theorem  and  the  generalized  Nyquist  theorem  by 
Desoer  and  Wang. 

The  results  that  we  obtain  in  Chapters  2-3  contain  the  finite 
dimensional  system  case  as  a  special  case.  In  particular,  we  showed 
that  the  guaranteed  stability  margins  of  the  distributed,  LQ  optimal 
systems  are,  generally  speaking,  the  same  as  those  of  the  finite 
dimensional  LQ  systems. 

In  Chapter  4  we  employ  a  flexible  beam,  simply  supported  at 
both  ends,  to  illustrate  the  use  of  the  above  robustness  results  to 
synthesize  a  closed-loop  stable  control  design  for  the  vibration 
suppression  problem.  Lyapunov  control  is  used,  which  offers  the 
advantages  of  being  computationally  less  complex  than  the  general  LQ 
method  with  better  stability  margins.  By  using  a  suboptimal  scheme 
which  calls  for  discrete  spatial  sampling  of  the  deflection  rate 
along  the  beam,  we  are  able  to  obtain  a  closed-loop  stable  finite 


dimensional  controller  when  there  is  a  sufficient  amount  of  inherent 
damping.  Specifically  we  show  that  we  can  design  such  a  controller 
if  the  modes  of  the  structure  has  a  constant,  albeit  small,  damping 
ratio.  The  fact  that  the  number  of  spatial  samples  required  to 
approximate  the  controller  increases  as  this  damping  ratio  decreases 
is  exposed.  The  case  of  constant  damping  coefficient  bears  similar 
results,  except  that  the  suboptlmal  controller  contains  spatial 
sampling  as  well  as  a  dynamical  compensator  to  attenuate  the  high 
frequency  components  of  the  loop  response. 

5.3  Suggestions  for  Future  Research 

We  believe  that  the  linear  quadratic  gaussian  (LOG)  approach  is 
a  practical  method  for  designing  feedback  controllers  for 
distributed  systems,  provided  that  the  control  system  designer 
understands  the  fundamental  characteristics  of  the  compensated 
system.  This  thesis  has  laid  down  the  basic  groundwork  for  the 
designer  to  do  a  credible  job.  However,  there  is  a  whole  realm  of 
other  design  techniques,  which  have  proven  to  be  successful  in 
finite  dimensional  system  design,  that  this  thesis  has  not  explored. 

We  have  studied  the  implication  of  full-state  LQ  feedback  for 
distributed  system  in  detail,  and  have  analyzed  the  use  of  spatial 
sampling  to  approximate  the  full-state  measurements. 

Depending  on  the  application,  one  or  more  of  the  following 
techniques  may  be  useful  for  control  designs  for  distributed 
systems.  When  combined  with  the  results  in  this  thesis,  these 


techniques  are  likely  to  produce  control  designs  which  have 
desirable  properties.  Theoretical  treatments  as  well  as 
applications  of  these  techniques  for  distributed  systems  will  allow 
us  to  obtain  a  unified  framework  for  dealing  with  infinite 
dimensional  system  design. 

(11  loop-transfer-recoverv  (LTR)  technique 

The  LQG/LTR  technique  [3],  [6]  is  one  of  the  most  popular 
methods  for  designing  controllers  for  multivariable  systems.  In 
Chapter  2,  we  mentioned  that  this  technique  seems  to  be  applicable 
to  linear  hereditary  differential  systems  also.  We  think  by 
properly  choosing  the  stochastic  filter  (observer),  the  same  can  be 
said  of  distributed  systems  in  general.  Although  in  theory  the 
LQG/LTR  method  is  limited  in  applicability  to  minimum  phase  systems, 
Athans  [6]  has  reported  that  desirable  results  are  obtained  in 
actual  applications  even  when  the  system  is  non-minimum  phase,  and 
has  recommended  the  use  of  this  method  under  all  circumstances. 
Therefore,  we  suggest  developing  this  method  for  infinite 
dimensional  systems,  and  carrying  out  a  feasiblity  study  of  this 
method  on  some  concrete  applications. 

(21  Design  Synthesis  Technique  that  Utilizes  Crossfeed  Structure 
In  multi-input-multi-output  systems,  it  is  possible  [13]  to 


adjust  the  control  weighting  matrix  to  trade  off  the  stability 
margins  against  the  various  crossfeed  perturbations  among  the 


control  channels.  Lee  and  Gully  [67]  have  applied  this  technique  to 
the  control  of  lateral  dynamics  of  a  high  performance  missile.  This 
technique  is  applicable  to  all  distributed,  linear-time-invariant 
systems  without  modifications. 

(31  Design  Analysis  technique  that  Explores  Error  Structure 

Lehtomaki ,  et  al .  [63]  have  shown  how  one  can  relax  the 
conservatism  when  evaluating  the  feedback  design  by  taking  into 
account  the  structure  of  the  perturbation  error.  As  a  result,  the 
control  system  designer  has  more  room  to  trade  off  the  various 
system  characteristics  in  his  design.  This  technique  is  applicable 
to  distributed  systems  without  modifications. 

A  notable  limitation  of  our  results  for  distributed  system  is 
that  the  control  space  is  finite  dimensional.  This  is  because  all 
the  known  Nyquist-type  of  criteria  have  the  same  limitation. 

However,  the  Yakubovich  theorem  is  valid  for  infinite  dimensional 
control  spaces  also.  Perhaps  by  exploring  a  function  other  than 
det(»),  the  robustness  results  in  Chapter  3  (which  requires 
verification  on  the  Nyquist  path  only)  can  be  extended  to  this  case 
also.  Until  this  suggested  research  direction  is  undertaken,  one 
must  contend  with  results  such  as  that  in  Appendix  B  which  requires 
the  control  designer  to  check  the  inver tibi lity  condition  over  the 
whole  CRHP,  instead  of  the  Nyquist  path  alone. 

There  are  two  important  types  of  modeling  errors  that  we  have 
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not  considered  in  this  thesis.  Quite  often  the  delay  time  in  a 


time-delay  systems  is  not  known  exactly.  It  is  useful  to  know  how 


much  error  in  this  delay  time  can  a  time-delay  system  tolerate 


without  becoming  unstable.  Also  we  have  not  considered  nonlinearity 


in  the  system  loop.  We  think  both  of  these  error  types  are 


important  and  more  work  should  be  done  to  expose  their  impact  on  the 


system  stability. 


In  this  research,  we  have  not  fully  analyzed  the  impact  of 


classifying  perturbation  errors  into  modeling  errors  and 


implementation  errors.  For  infinite  dimensional  control  designs,  it 


may  be  advantageous  to  distinguish  between  these  two  classes  of 


errors,  and  analyze  each  of  them  separately,  as  in  the  localization 


theorem  in  Lee,  et  al.  [64]. 


Lastly  we  would  like  to  suggest  that  the  control  problem 


formulation  in  Chapter  4  for  the  flexible  beam  can  be  extended  to 


control  problems  for  flexible  structures  with  multiple  components. 


The  exact  Lyapunov  solution  method  by  Djaferis  and  Mitter  [73]  seems 


to  be  suitable  for  handling  these  problems. 
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APPENDIX  A 


Derivation  of  Identitites  1-5  for  LQHD  Systems 

f 


Derivation  of  Identity  1: 


Kj(s) 


K1(0)eS0d0 


(A.l) 


Integrating  by  parts,  we  get 


Kj(s)  =  i  (K1(0)-KQBeST)  -  ±  J°  eS0  Kj(0)d0  (A.2) 


Derivation  of  Identity  2- 

Substituting  Eq.  (2.7)  into  the  integral,  we  readily  see  that 


J 


N-l 


-T 


eS0  Kx(0)d0 


s0. 


(A '  “K0BR-1B  *  ]Kj  (s)  +  l  K^e  i+K0A0(s) 

i=l 


i0 


-T 


eS0K2(O,0)d0 


(A. 3) 


where  K2(O,0)  is  given  in  Eq.  (2.9),  i.e. 


S0 


K2(O,e)d0  =  J°  es9Kj(-0-r)ANd0 


-  J0  eS0  K^-O+aJBR'^’KCaJdodO 


,0  N-1 

J!,  •“  i 

1=1 


{ 


AjKjfO+Oj)  ,  -t  <0+0i 


,  otherwise 


K)  N_1 

J!,  i 

i=l 


{ 


Ki(-0+0i)Ai  .  ei  <  0 


,  otherwise 


1°  es0  r°  { Ao«>Ki 
J-T  ” — T  0 


(OK,  (C+9)  ,  f>-0-r 


,  otherwise 


}df  d0 


{  ‘S9J°  <Kl 

J-T  J-T  0 


(C-0)Ao(O  .  C>0-T 


,  otherwise 


}dC  d0 


(A.4) 


On  the  right-hand-side  of  Eq.  (A.4),  it  is  easy  to  see  that  the  1st 
term 


eS0K  j  (-0-t  )Ajjd0  =  e"STKj(-s)AN 


(A.5) 


by  a  change  of  variable.  The  2nd  term  is  simply  ^(s)  as  defined  in 


Eq.  (2.21).  r^s),  r2(s),  Ej(s),  E2(s)  are  obtained  from  the  3rd  - 

6th  terms,  respectively,  by  applying  the  appropriate  upper  and  lower 
limits  to  the  integrals. 

Derivation  of  Identity  3: 


EJ 


*(s)  =  1  r  K^-e+oJBlrt'K^CTjdCTde 


Changing  the  order  of  integration  gives 


tf(s 


(A .  6) 


es0  Kj(-0+o)  BR^B'K^oJdedo  (A. 7) 


# 

V  <*“ 


By  a  change  of  variable  n  =  6+a,  we  get 


#(s 


)  =  J°  e_S”rK j  (”r )BR- !B '  Kj  (a) 


eSCTdtda 


=  t  J°  e'S7Kj(nr)d7]BR_1BIK1(a) 


sa, 
e  da 


(A.8) 
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and 


dV  =  BR~  *B ' Kj (a) esada 


we  obtain 


*(s 


)  =  [j°  e"S7Kj(nr)d7]  [J^  BR'^’K^Oe^dC] 


a  -  0 


ct  =  -  T 


j°  e-SCTKj(a)  J7  BR"1B,K1(f)esfdCda 


(A. 10) 


Evaluating  the  1st  term  on  the  right-hand-side,  we  see  that  it  is 

i  ”  1  1  i 


the  same  as  K^(-s)BR  B'K^(s).  Comparing  the  2nd  term  with  Eq. 
(A. 8),  we  see  that  it  is  actually  -^’(-s).  Therefore, 


""1dm 


^(s)+^*  (-s)  =  Kj  (-s)BR  B'K^s) 


(A. 11) 


as  claimed. 


Derivation  of  Identity  4: 

From  Eqs.  (2.22)  and  (2.23),  we  know  that 


** 


m 


s 
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iy  s) 


Ay^e+e^e^de 


N-i 

r2(-S)  =  ]>  J  AjK1(-0+0i)eS0d0 


1=1  '8i 


By  a  change  of  variable,  Eq.  (A.  12)  can  be  written  as 


ri(s)  =  l  J^1  AjKjtOe^dCe  ^ 

1=1  ~T 


Similarly,  Eq.  (A.  13)  can  be  written  as 


r2(-s)  =  l  ]°  AjKjOe^dfe  *8) 


1.1  8i 


Adding  Eqs.  (A. 14)  and  (A. 15)  produces  the  desired  result. 

Derivation  of  Identity  5: 

From  Eq.  (2.23)  we  know  that 


Hj(s)  =  £  A^fJKjCf+OJdfe^de 


(A. 12) 


(A. 13) 


(A. 14) 


(A.  15) 


(A. 16) 


WMIUINI  li.  PU/VUm^IUVU^  -SEW u 
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Changing  the  order  of  Integration,  we  get 


AQ(f )Kj (f +0)eS^d0df 


(A. 17) 


By  a  change  of  variable,  we  obtain 


AoCfJKjfOe^e'^dfdf 


(A. 18) 


Similarly  from  Eq.  (2.24),  we  get 


S2(_s)  =  £  A^CfJKjCOe^^dfdf 


(A. 19) 


Adding  Eqs.  (A. 18)  and  (A. 19)  gives 


EjCsJ+Sy-s)  =  J°  J°  A0(f)K1(f)esfe"sfdCdf  (A.20) 


and  the  claim  is  immediate. 


APPENDIX  B 


Robustness  Theory  for  Infinite  Dimensional  Control  Space 

The  purpose  of  this  appendix  is  two-fold.  First  we  give  a 

2 

tutorial  on  singular  values  of  linear  bounded  operators,  which  are 
useful  for  characterizing  stability  margins  of  infinite  dimensional 
systems.  Second  we  use  this  opportunity  to  show  how  one  can 
characterize  robustness  of  unity  feedback  system  in  general  by 
exploring  the  invertlbility  condition  within  the  Banach  algebra  of 
linear  bounded  operators.  This  result  differs  from  that  in  Chapter 
3  because  the  control  space  may  be  infinite  dimensional. 

B.l  A  Tutorial  on  Singular  Values  of  Linear  Bounded  Operators 

The  definitions  and  results  in  this  tutorial  are  largely  due  to 
or  are  easy  consequences  of  Gohberg  and  Krein  [68].  Also  we  have 
drawn  on  Ben-Israel  and  Greville  [44]  for  a  lot  of  the  motivation 
behind  the  definitions  and  interpretations  of  singular  values.  We 
have  limited  this  tutorial  to  contain  only  results  that  are  relevant 
to  the  study  of  robustness  in  this  appendix. 

For  the  development  of  singular  values,  we  list  the  following 
definitions  and  notations. 


Singular  values  are  referred  to  as  s-numbers  in  [68]. 


■swnwanooBownrai 
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H  =  a  separable  Hilbert  space  with  inner  product  (•.•) 

B(H)  =  set  of  all  linear  bounded  operators  on  H 
C(H)  =  set  of  all  linear  compact  operators  on  H 
I  =  identity  operator  H 
For  A  t  B(H) . 

|  | A  |  |  =  induced  operator  norm  of  A 
A  =  adjoint  operator  of  A  in  the  sense  that 
(x,Ax)  =  (A*x,  x)  for  all  x  e  H 
R(A)  =  range  of  A 
r(A)  =  dimension  of  R(A) 

Xj (A)  =  j-th  eigenvalue  (in  nonincreasing  order)  of  A,  where  A 
is  compact 

(A  complex  number  X  is  called  an  eigenvalue  of  A  if  (M-A)  $  »  0  has 

at  least  one  nonzero  solution  $  e  H,  i.e.  eigenvector) 

*  1/2 
Let  A  =  A  and  (x,  Ax)  >  0  for  all  x  e.  H.  Then  A  =  square 

root  of  A  in  the  sense  of 

*  =  A 1/2V/2 

B.1.1  Singular  Values  of  Linear  Compact  Operators 

Roughly  speaking  the  set  of  compact  operators  is  the  closest 
family  to  matrices.  Since  a  compact  operator  on  a  separable  Hilbert 
space  admits  a  Schmidt  expansion,  defining  the  singular  values  of  a 
compact  operator  is  almost  as  natural  as  defining  its  eigenvalues. 


Definition  (Singular  Values  of  a  Compact  Operator): 


Let  A  e  C(H) .  Then  the  eigenvalues  of  the  operator 

D  -  (A*A)1/2  (B.l) 

are  called  the  singular  values  of  the  operator  A. 

#  1/2 

Note  that  (A  A)  is  linear  compact  because  A  is.  Hence  its 
eigenvalues  are  denumerable.  We  shall  enumerate  the  nonzero 
singular  values,  which  are  all  nonnegative,  in  nonincreasing  order, 
taking  their  multiplicities  into  account.  Specifically,  let 


<?j(A)  =  Xj(D)  ;  J=1 .2 . r(A)  =  r(D)  (B.2) 

Oj(A)  =0  ;  j  =  r(A)  +  1,...  (B.3) 


Much  of  the  significance  of  singular  values  is  due  to  the  fact 
that  one  can  decompose  any  linear  compact  operator  into  a  series  of 
orthogonal  parts.  The  following  theorem  is  due  to  the  well-known 
Schmidt  expansion  of  compact  operators. 

Theorem  B.l  (Singular  Value  Decomposition): 

Let  A  e  C(H) .  Then  A  admits  a  Schmidt  expansion,  i.e.  there 
exist  {♦j}.  orthonormal  bases  for  II,  such  that  the  operator  A 

can  be  represented  by 


•wmwwwnu*  mo  5  «juu.uu 
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r(A) 

Ah  =  J  CTj(A)  (h.*  )*j;  *>  £  H  (B.4) 

j=l 


This  representation  of  A  is  called  the  singular  value  decomposition 
(SVD)  of  A.  The  series  summation  converges  uniformly  in  the  induced 
operator  norm.  Taking  multiplicities  into  account,  the  basis 
elements  ^  are  unique  for  j  =  l,...,r(A),  up  to  a  complex 


scalar  multiple. 


where 


=  set  of  all  finite  dimensional  operators 

of  dimension  smaller  than  k  (B.7) 

Moreover  the  solution  to  Eq.  (B.6)  is  given  by  a  partial  sum  in 
the  SVD  of  A, 


k-1 

*k  ■  l 

J=1 


Corollary  B.2 
(a)  a1(A)  =  ||A|| 


(B.8) 


(B.9) 


(b)  Let  r(A)  =  n  <  “.  Then  ^n(A)  is  nonzero,  but  on+j(A)  is  zero, 

for  j  =  1.2 .  Suppose  r(A  +  AA)  <  n-1,  where  AA  e  B(H) . 

Then 


1 1 AA 1 1  >  on(A) 


(B. 10) 


Moveover,  equality  in  Eq.  (B.10)  occurs  if 


AA  =  -  on(A)  (•.♦nHn  (B.  11) 

(c)  Let  A,  B  e  C(H) .  Then 


KKmmcnrwm  mm  w.u  wmvwwmv.w^jrtm^wm.mmm^irA  wv:*v rj  itj  i ru  *viiru  irv  nvm 
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ai+J-1(A  +  B)  <  ai(A)  +  aj(B)’*  1'Js1'2 .  (B.12) 

The  proof  of  Theorem  B.2  can  be  found  in  [68,  pp.  28-29]. 
Corollary  B.2  contains  easy  consequences  of  this  theorem.  Corollary 
B.2(a)  states  that  the  largest  singular  value  of  a  compact  operator 
is  actually  equal  to  the  norm  of  the  operator.  Since  singular 
values  are  similar  to  eigenvalues,  they  are  often  computationally 
more  tractable  than  using  notions  such  as  the  induced  operator  norm. 
It  seems,  therefore,  that  singular  values  are  useful  as  a  measure  of 
robustness,  provided  that  we  can  extend  Theorems  B.l  -  B.2  to 
include  linear  bounded  operators. 

Conceptually  the  appeal  of  Theorems  B.2  lies  in  the  fact  that 
it  relates  the  singular  values  of  a  compact  operator,  which  is  an 
algebraic  property,  to  a  geometrical  propoerty  in  the  space  of 
linear  bounded  operators.  Loosely  speaking,  the  SVD  gives  the 
shortest  distance  and  directions  to  the  "more  singular”  operators. 

If  H  is  finite  dimensional,  the  C(H)  =  B(H)  =  space  of  nxn 
matrices,  where  n  =  dim  H.  Let  A  be  an  nxn  matrix.  Then  A  is 
invertible  if  and  only  if  on(A)  >  0.  By  letting  i  =  n,  J  =  1  in  Eq. 

!  (B. 12)  and 

A  =  A  +  AA  (B.  13) 

B  =  -  AA  (B. 14) 


where  AA  is  an  nxn  matrix,  one  gets 
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cn(A)  <  aj(-AA)  +  an(A  +  AA)  (B.15) 

But  by  the  definition  of  singular  values 

a^-AA)  =  Oj(AA)  (B.16) 

Hence  if 


CTj(AA)  <  an( A) 


(B. 17) 


then 


an(A  +  AA)  >  0  (B. 18) 

i.e.  A  +  AA  is  invertible  also. 

We  shall  return  to  discuss  this  concept  when  we  discuss  the 
robustness  theorems. 

B.1.2  Singular  Values  of  Linear  Bounded  Operators 

In  general  linear  compact  operators  are  too  restrictive  for 
modeling  distributed  dynamical  systems.  In  this  section  we  extend 
the  notion  of  singular  values  of  linear  operators  to  include  linear 
bounded  operators.  The  main  difficulty  here  is  that  the 
"eigenvalues"  (more  precisely,  spectrum)  of  a  bounded  operator  do 
not  form  a  countable  set.  However  this  problem  can  be  resolved  by 


85 


5* 


decomposing  the  linear  bounded  operator  into  three  parts,  the  first 
of  which  is  a  linear  compact  operator. 

To  give  meaning  to  the  singular  values  of  linear  bounded 


multiplicity  m.  Define 


RSUUIIUWiJAPJ 
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Xj(D)  =  m  :  J  =  1.2 . m.  (B.22) 

and  the  remaining  Xj , 

Xm+J(D)  =  XjJDj)  ;  j  =  1.2....  (B.23) 

where 

Dj  =  D  -  pP  (B.24) 

and  P  is  the  orthoprojector  onto  the  eigenspace 
of  the  operator  D  corresonding  to  the  eigenvalue  p. 

Definition  (c)  gives  a  recursive  procedure  to  arrange  the 
largest  eigenvalues  of  D,  which  are  nonnegative,  in  a  nonincreasing 
order. 

By  this  construction  the  sequence  X^  has  a  limit,  which  we 
denote  by  Xro(D) .  Clearly, 

\o(D)  =  sup  spJD)  (B.25) 

Now  we  can  formally  define  the  singular  values  of  a  linear 
bounded  operator. 


£ 


Ar'  *»'  ■*!  ,  


Definition  (Singular  Value  of  Linear  Bounded  Operators): 

Let  A  t  B(H).  Then  the  singular  values  of  A  are  defined  in 
nonincreasing  order  as 


aj(A)  *  V,(D)  ;  j  =  1.2 . » 


(B.26) 


where  X  (D)  is  defined  in  Definition  (c)  above. 

J 

Figure  B.l  illustrates  the  relationship  between  the  spectrum  of 
D  and  the  singular  values  of  A. 


cJqM  )  (orthoprojector  representation) 


a.  (JJ )  (compact  operator 
J  representation) 


The  points  x  and  the  intervals  (/////) 

*  1/2 

are  contained  in  sp[(A  A)  ] 


Figure  B.l  Illustration  of  Singular  Values  of  Linear 
Bounded  Operator  A 


In  the  figure,  we  see  that  in  defining  the  singular  values  for 
the  operator  A.  we  have  broken  A  into  three  parts,  the  first  of 
which  is  a  compact  operator.  Hence  the  singular  values  can  be 
enumerated.  The  second  part  is  an  orthoprojector  on  a  plane  which 
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separates  the  eigenspace  of  the  compact  part  from  the  ''non-compact" 
part.  For  the  readers  who  desire  to  know  this  geometrical 
interpretation  more  precisely,  we  refer  them  to  [68,  pp.  60-61], 

The  quantity  ainf(A)  in  the  figure  will  be  defined  later. 

The  following  theorem  provides  the  motivation  behind  the 
construction  of  the  singular  values. 

Theorem  B.3  (Approximation  Property  of  Singular  Values'): 

Let  A  e  B(H).  Then  for  any  k  =  1,2, . . .  ,r(A), .  .  . 


a.(A)  =  min  |  |A  -  K.|  | 

K  K.  eF,  K 

k  k 


(B.27) 


where  F^  is  defined  in  Eq.  (B.7).  Moreover 


°«>(A)  =  mln  I  |A  -  K| 
KeC(H) 


(B.28) 


The  proof  of  this  theorem  is  given  in  [68,  pp.  61-62].  This 
theorem  has  an  important  corollary. 

Corollary  B.3  (Fundamental  Relationship  between  Singular  Values  and 
Norm) : 


a 


I 


(B.29) 


aj(A)  =  | | A  |  | 

Proof :  =  {0}  C  B(H).  Therefore  letting  k  =  1  in  Eq.  (B.27) 

gives  Eq.  (B.29).  (End  of  Proof) 

Referring  back  to  Figure  B.l.  we  define  the  leftmost  point  of 
the  spectrum  as 

Oinf(A)  =  inf  sp(D)  (B.30) 

We  will  show  that 

^(A'1)  =  [CTinf  (A)]-1  (B.31) 

First  we  establish  a  lemma. 

Lemma  B.4: 

Let  A  e  B(H)  be  invertible.  Then 

(a)  if  A  a  sp(A),  then  X  *  e  sp(A  *}  ;  (B.32) 

(b)  if  X  is  a  nonzero  limit  point  in  sp(A),  then  X  *  is  a  limit 

point  in  sp(A  *).  (B.33) 

Proof  of  Lemma  B.4: 

First  we  note  that  the  invertibi  1  i ty  of  A  implies  0  €  sp(A). 
Hence  XI  -  A  is  invertible  if  and  only  ifA*-X*Iis  invertible. 
Therefore  X  *  e  sp(A  *).  (b)  follows  from  (a)  because  the  function 

X  *  is  analytic  in  the  region  X  /  0. 


.'httu 


We  state  Eq.  (B.31)  formally  In  a  theorem. 


Theorem  B.5: 

Let  A  e  B(H)  be  invertible.  Then 

®lnf(A)  ■=  OjCA"1)]'1  (B.34) 


Proof  of  Theorem  B.5= 

*  1/2 

Since  sp[(A  A)  ]  contains  only  positive  numbers,  by  Lemma 

B.4, 


{inf  sp[(A*A) 1  =  sup  sp[(A*A) 


(B.35) 


But 


-1  M  -1/9 

o^A  )  =  sup  sp[(A  A)  ]  (B.36) 

because  A  ■*  A  is  an  involution  [46]:  i.e. 

(A^A)"1  =  A-1^"1)*  (B.37) 

and 

a1(A~1)  =  Oj[( A_1)M]  (B.38) 


.A  y.  v.  <./.<•  <  •-  r.  /.  a.  /, 
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Thus  by  the  definition  of  ajnf(A)  in  Eq.  (B.30) 

[®,  nf(A)D_I  =  ajlA'1)  (B.39) 

(End  of  Proof) 

This  completes  our  discussion  on  singular  values  of  linear 
bounded  operators. 

B.3  Generalized  Robustness  Theorem 

In  this  section  we  present  a  small  gain  theorem  for  additive 
perturbations,  and  use  Theorems  B.3  and  B.5  to  develop  a  generalized 
robustness  theorem  for  feedback  systems.  We  adopt  the  definition  of 
feedback  systems  in  Willems  [53]. 

B.3.1  Definitions  and  Notation 

We  shall  use  the  following  notation  in  this  section: 

H  =  a  separable  Hilbert  space  with  an  inner  product  (•,•) 

F  =  a  separable  Hilbert  space  of  functions 
f:  [0 ,«)  -»  H,  with  inner  product  <•«•> 

P^,  =  truncation  operator  on  F  defined  by 

’  f(t)  ;  t  <  t 

(PTf)(t)  =  ;  f(.)  e  F 

.0  ;  t  >  T 
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Fg  =  extended  space  of  F 

=  {£(•) |PTf  £  F  for  all  0  <  t  <  »} 

I  =  Identity  operator  on  Fe(or  F) 

B(F  )  =  set  of  all  linear  bounded  operators  on  F 

B(F)  =  set  of  all  linear  bounded  operators  on  F 

Let  G  e  B(Fe)  (B(F)).  Then 

G*  =  adjoint  of  G  with  respect  to  <•«•> 

|  |C|  |  =  Induced  operator  norm  by  <•,•> 

Let  G  be  an  operator  on  Ffi(F).  Then  G  is  called  causal  if 

P  GP  =  P  G  for  all  t  >  0,  where  P  has  been  extended  to  F  .  Let 
r  t  r  -  t  e 

B+(F  )  =  {G  e  B(F  )  |G  is  causal} 

6  6 

and 

B+(F)  =  {G  e  B(F) |G  is  causal} 

With  an  abuse  of  notation,  given  G:F  -»  F  ,  the  restriction  of 

e  e 

G  to  F.  G|F.  is  also  denoted  by  G  when  no  confusion  arises. 

Definition  (Nominal  Unity  Feedback  Control  System!: 

In  Figure  B.2 

y  =  Ge  (B.40) 


e  =  u-y 


(B.41) 
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G:  F  -*  F  is  causal 
e  e 


(B.42) 


The  system  in  Figure  B.2  is  called  the  nominal  unity  feedback 


system.  It  is  called  linear  if  G  is  linear. 


Definition  (Uni tv  Feedback  System  subject  to  Additive  Perturbation 


In  Figure  B.3 


y  =  (G  +  AG)e 


(B.43) 


e  =  u-y 


(B.44) 


G:  F  -*  F  is  causal 
e  e 


(B.45) 


AG:  F  -*  F  is  causal 
e  e 


(B.46) 


The  system  in  Figure  B.3  is  called  a  unity  feedback  system  subject 


to  additive  perturbation. 


Definition  (Wei  1-Posedness  of  Feedback  Systems'): 


Referring  to  the  unity  feedback  system  in  Figure  B.4,  let  F'-  F 


-♦  Fg  be  causal.  Then  the  basic  feedback  equation  is 


(I  +  F)e  =  u 


(B.47) 
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Figure  B.4  Unity  Feedback  System 
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(a)  We  say  that  the  unity  feedback  system  is  well-posed  if 

(I  +  F)  F“e  -»  Fe  exists  and  is  causal  (B.48) 

(b)  Given  that  the  unity  feedback  system  is  well-posed;  we  say  that 
the  unity  feedback  system  is  (finite  gain)  stable  if  the 
inverse  in  Eq.  (B.48)  is  bounded  when  restricted  to  the  space 
F,  i.e. 

(I  +  F)  * |F  is  bounded  (B.49) 


Remarks: 

(1)  Referring  to  the  definitions  of  the  nominal  unity  feedback 
system,  and  system  subject  to  additive  perturbation,  two  cases 
are  of  interests: 

a)  F  represents  the  nominal  plant  operator,  i.e.  G:  F  ■*  F  ; 

c  © 

b)  F  represents  the  nominal  plant  operator  plus  the 
perturbation  operator,  i.e. 

F  =  G  +  AG:  Fe  -»  Fe  (B.50) 

(2)  A  unity  feedback  system  not  satisfying  the  well-posedness  is 
ill  conditioned.  Therefore  in  the  rest  of  this  appendix,  when 
we  speak  of  a  unity  feedback  system,  well-posedness  is  always 
assumed  implicitly. 


‘'■■■'fa 


MM 


(3)  Suppose  that  G,  AG  e  B+(Fe),  and  assuming  well-posedness,  then 
the  nominal  unity  feedback  system  is  stable  if  and  only  if 

(I  +  C)'1  e  B+(F)  (B.51) 

and  the  unity  feedbak  system  subject  to  additive  perturbation 
AG  is  stable  if  and  only  if 

(I  +  G  +  AG)"1  e  B+(F)  (B.52) 

This  is  true  since  a  causal  operator  on  F  can  be  extended 
uniquely  to  a  causal  operator  on  F^. 

(4)  The  notion  of  stability  here  says  that  a  bounded  input  produces 
bounded  output. 

Given  that  the  operators  G  and  AG  are  bounded  and  the  nominal 
unity  feedback  system  is  stable,  a  basic  and  legitimate  question  to 
ask  is: 

"What  is  the  region  of  robustness  of  stability  of  the  unity 
feedback  system?"  In  other  words,  "what  is  the  range  of  AG  for 
which  the  perturbed  system  is  (guaranteed  to  be)  stable?" 

One  answer  to  this  question  is  the  following  theorem. 


Theorem  B.6  (Small  Gain  Theorem  for  Additive  Perturbation) : 
Consider  the  feedback  system  in  Figure  B.3  (assuming 
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.  n  .  »»-'■ .  Kii*-  . 


.  iT  ,4«  ■  >«  ii  .fr»  .>J  .<■  ■  >«  -  , 


»'n#«rl^«r»inr,jrw»A/1.i^jwv?^-4..^:i;.irii.»rvf  r-J'.rv  nuru  av«vi  ia 


well-posedness) .  Suppose  G:  F  ->  F  and  AG:  F  -*  F  are  both  linear 

©  ©  ©  © 

and  causal  and  the  nominal  unity  feedback  system  is  stable,  i.e. 
there  exists 


(I  +  G)-1  t  B+(F) 


(B.52) 


Then  a  sufficient  condition  for  the  perturbed  unity  feedback  system 
in  Figure  B.3  to  be  stable,  i.e.  there  exists 


(I  +  G  +  AG)"1  e  B+(F) 


(B.53) 


||(I  +  G)_1 AG | |  <  1 


(B.54) 


Corollary  B.6 


Another  sufficient  condition  is 


| |(I  +  G)-1 1 1  |  |AG| |  <  1 


(B.55) 


-1 i i-l 


I  I AG 1 1  <  ||(I  +  G)-i 1 1 


(B.56) 


Proof : 


Since  F  is  a  Hilbert  space,  under  the  induced  operator  norm 
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B(F)  forms  a  Banach  algebra  (see  e.g.  Rudin  [46]  or  Willems  [53]). 
Moverover  B+(F)  C  B(F)  is  a  subalgebra.  Therefore  if 


||(I  +  G)“1AG| |  <  1  ( B . 57 ) 

then  I  +  (I  +  G)  *AG  is  invertible  in  B+(F).  Hence  I  +  G  +  AG  is 
invertible  in  B+(F),  since  the  invertible  elements  of  B+(F)  forms  a 
group.  Therefore  the  unity  feedback  system  in  Figure  B.3  is  stable. 
The  corollary  follows  since  under  this  norm 


|  |(I  +  G)-1AG|  |  <  ||(I  +  G)"1  1 1  1 1 AG |  | 


(B.58) 


(End  of  Proof) 


Remarks : 

(1)  This  theorem  shows  that  the  region  of  robustness  of  stability 
can  be  characterized  by  the  gains,  i.e.  induced  operator  norm 
of  the  inverse  of  I  +  C  and  AG. 

1 1 AG 1 1  <  ||(I  +  G)_1 1 1"1  (B.59) 

This  equation  will  form  the  basis  for  the  development  of  the 
robustness  theorem  in  terms  of  singular  values  below. 

(2)  The  requirement  that  F  be  a  separable  Hilbert  space  can  be 
relaxed  to  a  Banach  space,  but  the  former  is  needed  for  the 
development  of  singular  values. 


The  following  robustness  characterization  in  terms  of  singular 
values  is  a  direct  consequence  of  Theorems  B.6,  B.3  and  B.5. 


Theorem  B.7  (Small  Singular  Value  Theorem  for  Additive 


Per  turbat i on) : 


Under  the  conditions  of  Theorem  B.6,  a  sufficient  condition  for 


the  perturbed  unity  feedback  system  in  Figure  B.3  to  be  stable,  is 


OjfAG)  <  ainf(I  +  G)  (B.60) 

By  using  the  Parseval's  theorem,  one  can  easily  show  that 
Theorem  B.7  implies  the  usual  finite  dimensional  robustness 
characterization  results  (Sandell  [74],  Lehtomaki  [5],  Doyle  [10]). 

This  completes  our  discussion  on  robustness  and  the  singular 
values  of  linear  bounded  operators. 


APPENDIX  C 


Nodal  Expansion  of  Flexible  Beam 


First  we  use  the  principle  of  separation  of  variables  to  solve 
Eqs.  (4.7)  -  (4.9).  Let 


y(x,t)  =  )  X.(x)T.  ( t) 

k=0 


Then  Eq.  (4.7)  can  be  written  as 


(C.l) 


w  4  2  2 

l  £7?  -  r  jt  Vx>  arV*> +  V«>  Jr 

k=0  dx  dx  dt 


=  u(t)  6(x  -  i) 


(C.2) 


Now  try 


”2  Xk(x)  +  ^k(x)  =  0 
dx 


(C.3) 


Then  Eq.  (C.2)  becomes 


00 


l 


at 


=  U(t)  6(x  -  j) 


(C.4) 


A.  L&..1 


t..  t .  -  ».  l-f.  ->rft  *W  ■»>.  «  «  W  1*  »  A  t  A,  •»«%  «  !■»■  «  li  **!'» 


ranwm  avkwkj^kkk  w« »»  *»  '*  » 


Using  the  boundary  conditions,  Eqs.  (4.8)  -  (4.9),  in  Eq.  (C.3),  we 


<d,  =  kTT 

k 


;  k  =  0.  1.  2. 


(C.5) 


X^(x)  =  sin(kTrx)  ;  k  =  0,  1,  2,... 


(C .  6) 


Since  v5  sin(krrx),  k  =  1,  2,...  forms  a  complete  orthonormal  basis 


for  the  beam  shape  between  x  =  0  and  x  =  1,  satisfying  the  boundary 


conditions,  we  conclude  that  the  solution  to  Eq.  (4.7)  -  (4.9)  is 


given  by 


y(x.t)  =  )  sin(kirx)T,  ( t) 
k=l  K 


(C.7) 


where  ( t )  is  given  by  Eq.  (C.4).  This  same  result  can  also  be 


obtained  by  using  the  Fourier  series  expansion. 


Now  multiplying  both  sides  of  Eq.  (C.4)  by  sin(kirx),  and 


integrating  over  x,  from  0  to  1,  we  obtain 


“JWO  ♦  fujj  ^  Tk(t)  +  2-5-  yt)  =  I  u(t)  sin(£L)  (C 


2  ~  2 


r  * 


Vif. 
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APPENDIX  D 

Green's  Function  Characterization  of  Kg 

The  Green’s  function  in  Eq.  (4.34)  is  determined  as  follows. 
Try 


h(x.f)  =  )  a  sin(iTrx)  sin(iTrf)  (D.  1 ) 

i=0 

where  a^'s  are  to  be  computed.  Multiplying  both  sides  by  sin(mrf) 
and  integrating  over  the  interval  between  f  =  0  and  £  =  1,  we  get 

[  )  a.sin(iTrx)sin(iirf )]sin(mrf )df 
i=0  1 

(D.2) 

By  using  Eq.  (4.34)  and  the  orthonormality  of  the  sine  functions, 
this  is  reduced  to 


f  h(x,f )sin(nirf )df  =  [ 
J0  J0 


a 

Kg  sin(mrx)  =  tj—  sin(mrx) 


(D.3) 


But 


— k  sin(mrx)  =  -  n  w  sin(mrx) 
dx 


(D.4) 


Therefore, 


o  2  2 

a  ann  n 

Kg  — g  sin(rnrx)  = - 2 —  sin(nTrx) 


(D.5) 


But  from  Eq.  (4.30), 


K„ 


— s-  sin(nirx)  =  -  sin(mrx) 
chi 


(D.6) 


This  implies  that  the  value  of  a  must  be 


a 

n 


2 

2  2  ’ 


n  ir 


n  =  1 ,  2,  3 


(D.7) 


Since  sin(nTrx)  forms  a  complete  basis  for  S(0,1)  and  all  the 
boundary  conditions  are  satisfied,  we  conclude  that  the  operator  Kg 
can  be  characterized  in  the  form  of  Eq.  (4.34),  by  using 


h(x.f) 


.2  2 
i  t r 


sin(  Ittx) 


sin(imf ) 


(D.8) 


It  is  easy  to  verify  that  K0  thus  defined  is  positive  definite 


U  Tywn1  TBimTf -f  ^  WML  W-HBnTi^*' 
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APPENDIX  E 

Optimal  Loop  Function  Analysis  < 

I 

i 

One  can  write 


g(s)  =  J  gk(s) 
k=l 


(E.l) 


where 


Ku(s)  = 


.  2Ak7r. 
sin  (s-) 


2s 


2  _.  2  2it4  4 
s  +£k  it  s+k  ir 


(E.2) 


2  2 

First  we  examine  |g(s)  |  at  s  =  jk  it  .  k  =  1,2,...  Separating  gk(s) 
from  the  series,  we  get 


lg(s)|  <  |gk(s)|  +  |hk(s)  | 


(E.3) 


where 


Vs>  -  1 


1  .  2  ,imr. 

sin  (^-) 


2s 


m=l 

m*k 


2  2 
m  7 r 


2  ’  2. r,  2  2  ,,  4  4 

s  +fk  it  s+k  IT 


(E.4) 


Claims* 


(1)  For  k  odd, 


max|gk(jw)l  = 


0) 


|gk(  jk2'2)  I  =  —5-4 

fk  1 r 


(2)  For  k  even,  gk(jco)  is  identically  zero, 


(3)  |hk(jk27T2)|  <  -|"2  (1  +  — i-o) 

k  ir  4k  ir 


Proof  of  Claims: 


(1)  For  k  odd,  since  the  function, 

n  ■ 4  4  2.2-2.  4  4  2  .  . 

1  ,,  ^  i-2  (k  ir  -<o  }  +C  K  ir  co  „4  4, 

|gk(J«)  I  =  * - -  (k  "  ) 

4co 

kV  rr2,  4  4  1  ,.4  4  2.2, 

=  -4—  tt  +  “2  (k  tt  -o  )  ] 

<0 

achieves  a  minimum  when  the  second  term  in  the  sum  is  zero. 

2  2 

co  =  k  ir  ,  |gk(  jco)  |  achieves  a  maximum  there  also,  and 


(E.5) 


(E.6) 


i.e.  at 
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kTT 

(2)  This  is  trivial  since  sin  {^— )  =  0  for  k  even. 

(3)  Let  us  split  h^s)  into  two  parts, 


k-1 


Pk<*>  *  I 


m=l 


1  _,_2  ,imr,  2s 

2  2  Sln  '2*2  22  44 

m  ir  s  +fm  ir  s+m  ir 


f  \  X  1  .  2  ,imr, 

qk(8)  =  2  “T2 8ln 


2s 


m=k+l 


m  ir 


2,r  2  2.44 
s  +|,m  ir  s+m  ir 


2  2 

At  s  :  jk  r  .  taking  the  magnitude  of  each  term,  and  using 
that  sin2(rp)  <  1, 


|pk(Jk2ir2)|  <  J 


2  4  A  A  2  A  A 
m=l  [ ( m  — k  )  k  j 


1/2 


Droping  fm^k^  from  the  demonlnator. 


k-1 


|pk(Jk2ir2)|  <  ^ 


m=l 


2  4 
m  ir 


.4  4 
k  -m 


Now  observe  that  by  partial  fraction  expansion. 


2_  k2  _  1_  ,2__  1 

2,4  4  “  ,  2  '  2  ,2  I 

m  k  -m  km  k  -m 


iA  2 
k  +m 


) 


(E.8) 


(E.9) 

the  fact 


(E.10) 


(E. 11) 


(E- 12) 
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Therefore, 


|pk(Jk2ir2)  |  <  -y-4 
k  ir 


,2  2‘ 

k  -m 


(E. 13) 


By  the  use  of  table  for  series  [71], 


k-1 

l 

m=l 


.2  2 
k  -m 


X  _JL_ 
A.  2  Z  2 .2 

414  m=k+l  m  “k 


/  2 
2ki+r 


00 


< 


(E.14) 


Also, 


(E. 15) 
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Thus 


|pk(  Jk2Tr2)  |  <  -i-j 

k  2k V 


(E. 16) 


Similarly, 


lqk( Jk  *  ^  <  „.  2  2  +  .  4  4 
2k  ir  4k  v 


(E.17) 


Now  since 


l\(s)l  <  |pk(x)|  +  |qk(s)| 


(E. 18) 


we  obtain, 


IhjJjk  ir  )|  <  2  2  +  4  4 

k  ir  4k  ir 


(E. 19) 


establishing  claim  (3). 

Using  the  above  claims, 


2.2, ,  ,  1 _  „  ^  3_a  a  I  fkV 


|gk(JkV)|  <i—(i  +  -fr>  + 
k  ir  4k  ir 


,  k  odd 


0  k  even 


(E.20) 


For  k  relatively  small,  and  odd,  since  f  is  small. 


|g(jk2ir2)|  ~ 

Ck  ir 


(E.21) 
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But  for  k  sufficiently  large,  the  term  — y-r  no  longer 

CkV 

1  3 

dominates  over  — »  a  (1  +  — s-jr).  ConsetIuently>  the  values 
K  7T  4k  IT 

2  2 

|g(jk  »  )!  are  small  at  high  frequencies  and  rolls  off  at  a  rate  of 
no  less  than  20  dB/decade. 

It  Is  not  hard  to  show,  by  modifying  the  above  analysis,  that 
2  2 

for  u  b  Jf  ir  ,  where  f  >  0  is  not  necessarily  an  Integer, 


|g(Jw)l  <  1 


,2  2 

f  IT 


(1  ♦  "%)  + 
4ir 


2  r 


o  >  ^  .  2  A  A  1/2 

k2w4((k4-f4)  +Ck4f4) 


(E.22) 


where  k  is  the  nearest  odd  integer  to  f.  Therefore,  on  the  whole. 

|g(J«)  I  rolls  off  at  20  dB/decade  at  high  frequencies,  and  there 

exists  a  cutoff  frequency  <j  . 

c 


APPENDIX  F 


Derviatlon  of  Suboptlmal  Loop  Function 

Recall  from  the  open-loop  analysis  that 

§£(*,t)=  ^  sln(lirx)  T  (t)  (F.l) 

Ol  ife1  1 

Hence  z^(t),  defined  in  Eq.  (4.52),  can  be  written  in  terms  of  T^t) 
as 

ze(t)  =  J  sin(^)  TjC t) ;  £  =  0.  1 . N  (F.2) 

Substituting  this  into  Eq.  (4.53)  gives 


5in<iFL)  Tj(t) 


(F.3) 


Combining  this  result  with  Eq.  (4.51),  and  translating  the  result  in 
the  frequency  domain,  we  get  the  suboptlmal  feedback 
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where  for  the  moment  we  assume  M  =  «  in  Eq.  (4.51).  We  shall  remove 
this  assumption  later. 


Using  the  transfer  function  T1(s)/u(s)  from  Eq.  (4.14),  we 
obtain 


v(s)  ?  2 _  .  ,k ir, 

*77  *  1  72  2  sln<r> 


u(s) 


k=l 


k  w 


N-l  _  .  ,iir, 

rl  V  ^  .  ,ikir.  .  1  far «  2s  8ln^  2^ _ , 

X  2n  1  sin*  N  *  8ln^  N  ^  2  .  _,2  2  .  .4  4^ 


«S0  151 


A  mm  ,  a*  A#  .  .  A  J 

s  +  fi  w  s  ♦  i  w 


(F.5) 


Rearranging  the  above,  we  get  (replacing  the  right-hand-side  with 

g(s)). 


g(s)  = 


2s  sin(^) 


kir. 


t-t  mm  j  mm  tm  J  * 

S  +  fl  IT  S  +  1  IT 


2  ,44^2  ,22  sin^2  ^ 


k=l 


k  TT 


N-l 


r  1  V  .  ,lkir,  .  ,if7r,,. 

x  tjr  f20  sl"(— >  sln(— >]> 


(F.6) 


where  the  order  of  the  summations  over  1,  k  and  t  have  been  changed. 

This  is  possible  because  the  summation  over  l  is  finite,  and  those 

over  i  and  k  are  easily  seen  to  be  uniformly  convergent  because  of 
4  4  2  2 

the  i  ir  and  k  ir  terms  in  the  denominators. 
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00 

1  X2  C08[k,r(i  +  J)3 

.  k  t 


k=l 


1  1  ,1  .  «.  .  1  rl  .  £.2 

6  ‘  2  (2  *  if5  *  4  (S  + 


(F.ll) 


00 


I 

k=l 


7T2  “■CMj  -  f)3 

k  it 


=  < 


I-i(i 

6  2  l2 


£)+!(!_  £)  *  <  N 

N;  4  l2  N'  •  e  -  9 


1+1(1 
fi  O  vo 


«,  + 1  (i .  i)2  a  (i 

nJ  a  VO  n2  •  9  N  e 


6  2  v2  N'  4  v2  N 


Therefore, 


N 


c,  = 


N-l 

n  2  sln<— > 
e=o 


e 

,  <  N 

2N 

•  e  <  2 

i  e 

N  < 

2  "  2N 

•  2  <  * 

Let  N  be  even  and  r  is  the  integer  N/2.  Then 


N-l 


N 

ci ■ 


-  hj  I  {sln(*r>  *  3  I  {sl"(*r> 


2N 


e=o 


e=o 


N-l 


1  V  .  ITT,  1  V  .  /^iTTx 

2N  Z  sln(  N  *  "  2N  Z  Sln^  N  3 


(F.12) 


(F. 13) 


(F. 14) 


I 


vWvfl 
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Using  Cables  [71] 


N-l 

^  ,  ,£i  ir%  ,  ,iir\  ,  /(N-Diir^  ,iu\ 

2  sin(-j-p)  =  sin(^— )  sin(H^ — )  cosec(^) 

€=0 

This  can  be  reduced  to 

N-l 

I  8in(T>  =  COt(^ 

e=o 


Also  from  tables  [71] 


l 


e=o 


.  .  ,H V\ 
€sin(-^) 


sln(  It t) 

A  .  2,i tr. 
4Sin  (w) 


m  /2N-1  .  . 
Ncos(-^j-  Itt) 

„  .  2,iir. 

2sin 


This  can  be  reduced  to 


N-l 


£=0 


N  ,  ,,i  .,iiF\ 
=  2  (-1)  cot(2N> 


For  the  other  terms  in  Eq.(F.M), 


I 


e=o 


.  . £iw x  .  ,r+l  .  ,  .  riTr  ,i7r, 

sln(— )  =  sin(~  lir)  sin~  cosec(ow) 


'2N 


2N 


l2NJ 


(F.15) 


(F  - 16) 


(F.17) 


(F* 18) 


(F- 19) 


(F.20) 


sin  “  lir  (r+1)  cos  "2n“  i7r 

4sln2(^)  2s,„2(ii) 


But  by  construction, 

r  -  |  (F.21) 

Hence  Eqs.  (F.19)  -  (F.20)  can  be  simplified  to,  for  i  odd, 

r 

l  aln(i^)  .  I  [cot  (g)  ♦  alnfil-)]  (F.22) 

e=o 

and  for  i  odd, 

1-1 

r  “2“ 

l  «sin(^L)  =  (-1)  cosec2(^)  +  |  (N+l)]  (F.23) 

e=o 

Therefore,  for  i  odd 


,  ,  A  r  1  2  f  Itt »  A  N+l 

(-1)  [— o  cosec  (d  +  — s- 


(F.24) 


Now  observe  that  Cj  above  is  obtained  by  summing  over  k  for 
infinitely  many  terms  in  Eqs.  (F.8)  -  (F.12).  In  actual 


Implementation,  this  is  not  possible,  and  we  can  only  compute  the 

series  from  k  =  1  to  k  =  M  for  some  finite  M.  However,  we  can 

always  choose  M  sufficiently  large  such  that  the  error  introduced  by 

this  truncation  is  less  than  — ^r,  for  each  B  in  Eq.  (F.13). 

4JT 

Therefore,  since  B  raises  from  0  to  N-l,  the  total  truncation  error 

N  1 

between  the  actual  and  nominal  Cj  in  Eq.  (F.24)  is  less  than 

This  gives  us  Eq.  (4.58),  in  which  the  symbol  O(j^)  is  used  to  denote 
N+l 


It  is  important  to  note  that  in  the  design  procedure,  H  is 
computed  only  after  N  has  been  chosen  (Appendix  G) . 
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APPENDIX  G 


Analysis  of  Implementation  Error 


Consider  the  difference  between  the  corresponding  terms 


g^s)  -  gjs). 


g^s)  -  g^s)  = 


,  2/  iTT x 

2s  r _ j _  / llTv  _N  sin  (2  K 

2  2  2  4  4  '■*  n'  2  Ci~  2  2  J 


s  +  fi  7T  S  +  i  TT 


(G.l) 


Clearly, 


max  Ig^Jw)  -  g1(J<j)  |  =  |gt( Jcj)  -  gjCjw)!  2  2  (G.2) 


U=i  V 


This  is  identically  equal  to  zero  for  i  even. 


By  performing  an  analysis  similar  to  the  one  in  Appendix  E,  we 


/v 

see  that  at  low  frequencies,  for  f  very  small,  |gj(ju)  -  gj(j<j)|  is 


much  larger  than  )  |g  (ju>)  -  g  (ju)|  around  u  =  i  ir  ,  i.e. 
m  i  n  n 

n=l 


|g(ju)  "  g(Jw)l  ~  |g± ( j^)  -  gjCJw)!  around  u  =  12tt2  (G.3) 


but  this  is  not  the  case  at  high  frequencies. 
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Let  us  examine  this  more  closely.  For  all  i  odd, 


|gj(J")  -  g,(J“)|  22  =  rT2lcl 

cj=i  ir  fi  7T 


1-1 

(-D  2  I 
2  2  1 
if 


(G.4) 


Let  6  >  0  and  let  Nq  be  chosen  such  that  for  N  >  Nq 


-2_<  ££ 

A  2 


(G.5) 


Then  since 


N+l  ,  1 

2  *  N 
2TT  " 


(G.6) 


we  have 


_2_  ,N±U  <  £6 
i  ir  2ir  4 


(G.7) 


So 


a  sufficient  condition  for  |gj(ji2ir2)  -  gj( Ji'V2)  |  <  6  is,  by 


,2  2, 


N 


using  Cj  from  Eq.  (F.24), 


cosec2(— ) - —I  <  ££  for  all  i  odd 

if  4NZ  2N  i  V  4 


(G.8) 


In  what  follows,  we  show  that  N  can  be  chosen  such  that  Eq.  (G.8)  is 


satisf led. 
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It  is  easy  to  see  that  as  N  tends  to  ®,  for  fixed  1, 

I  2  /lTr\  1  I  sy  f\  \ 

cosec  (gj-)  - •  -5-2  (0.9) 

4N  i  7T 

In  particular,  we  can  choose  a  small  number  e  >  0  and  Nj  such  that 
for  N  >  Nj 

l“V  cosec2  (^)  -  K?\  <  (G.  10) 

4N  7 Tir 


2  lir 

Since  the  maximum  of  the  discrete  sequence  cosec  (^-)  occurs  at 
1=1  and 


1  2  .  i7T»  1  2, IT  v 

max  — *■  cosec  (5^)  =  — s-  cosec  (srr) 

i  4N  4N 


(G- 11) 


We  deduce  from  Eqs.  (G.10)  and  (G.U)  that 


1  .  l+e6  .  . 

— o  cosec  (^jrr)  <  — sr-  for  all  i 


4N 


k2N'  2 

IT 


(G- 12) 


Therefore  for  all  i, 


1  1  2,iw\ 

l—o  cosec  (sjt) 

4N 


1  _ 1  *  fl+e.6  1  x 

2  2 1  ^  iriax  (  2  '  2  2 


1  7T 


iir 


(G.13) 


because  both  terms  on  the  left-hand  side  are  positive.  Clearly  for 
all  i  >  1,  the  quantity  on  the  right-hand  side  is  — 

IT 


^  Wi  1  ^  r_W  1 
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If  1  is  large  enough  such  that 


2ff2  (i|^)1/2<  1 


(G. 14) 


then 


l+e6  .  [6 
,2  4  M 

i  IT 


(G. 15) 


Thus  according  to  Eq.  (G.13),  for  these  i's, 


1  |  cosec2(i|)  -  -T2  I  < 

i  ir 


,2  2  1  Afl2 

i  ir  4N 


(G- 16) 


Let 


I  =  the  first  integer  larger  than  2 ir2  (^£-)^2 


(G- 17) 


Then  what  remains  for  us  to  do  is  to  make  sure  that  our  spatial 
sampling  is  dense  enough  that  Eq.  (G.16)  is  satisfied  for  i  <  I 
also.  Specifically,  we  want  to  choose  Ng  such  that  for  N  >  Ng 


1  I  1  2,  iir.  1 

.2  2  '  AN2  COSCC  W  '  .22 

i  ir  4N  i  ir 


< 

'  2 


(G- 18) 


for  all  i  <  I. 
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Finally  let 

N  =  max(N0.  Nr  N2)  (G.19) 

where  Nq,  N^,  Ng  are  defined  along  with  the  inequalities  in  Eqs. 
(G.7),  (G.16)  and  (G.18),  respectively.  This  N  guarantees  that 

|g(Jfa>)  “  g(Jw)|  <  6  for  all  t>  £  R+  (G.20) 

and  the  suboptlmal  system  is  closed-loop  stable. 

Next  we  establish  that  for  small  f,  t.  and  6 

l  IgjCJu)  “  gjfJuJl  <  0.05  (G.21) 

1=1 

i*k 

where  k  is  the  nearest  odd  integer  to  f  >  0,  and  f  is  related  to  w 


u 


.2  2 
f  t r 


(G.22) 


Note  that  by  Eqs.  (4.56)  -  (4.60), 


Ig^ju)  -  gA( J<^)  I  = 


2f 


A  A  2  A  A  1/2  O  ' 

[(f4-i4)  i4f4]  IT2 


gf- 


.4  .4,  2  fci 

f  -i  |7T 


(G.23) 
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where 


4N 


eosee  2  2 

i  TT 


°«(R)I 


By  Eqs.  (G.13).  (G.5).  (G.6), 


(G .  24) 


for  all  i 


(G.25) 


Since  t,  6,  and  f  are  usually  small  numbers,  a  is  very  close  to 

TT 

Also, 


o 

2  r  11 

f4-i4  "  f2-!2  f2+i2 


Clearly 


I 

i=l 

i*k 


f2  +  i2 


s  1 

i=l 


(G.25) 


(G.26) 


and  using  similar  analysis  to  that  in  Appendix  E,  we  obtain 


l 

1=1 

i*k 


ir2-i2i 


3 

4 


(G.27) 
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Therefore,  summing  Eq.  (G.23)  over  i,  and  using  Eqs.  (G.24)  - 
(G.27),  we  get 


1 

1=1 

i=k 


IgjUw)!  <  (g*  +  “\)  a 

4ir 


(G.28) 


~  1 


For  a  Z  which  is  usually  the  case,  the  last  quantity  in  Eq. 

IT 

(G.28)  is  less  than  0.05. 
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